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The  image  shows  sediment  scattering  from  a  mine  hunting  SONAR .  The  position  of  the  SONAR  is  indicated  by  the  red  star.  The  four 
panels  show  the  scattered  field  0. 15,  0.30,  0.45  and  0.60  seconds  after  the  SONAR  has  transmitted  its  ping.  The  forward  scattered  field 
is  shown  to  the  right  of  the  SONAR  position,  while  the  backscattered  field  is  shown  to  the  left,  (see  pages  87  to  94). 


See  pages  53  to  57  for  a  description  of  the  figures  on  the  front  cover 


Annual  Progress 
Report  1999 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


20010515  044 


CONTENTS 


Foreword  by  SACLANTCEN  Director  1 


Project  03-A:  Detection  and  classification  of  buried  mines  27 
Project  03-B:  Remote  large  area  search  in  coastal  waters  35 
Project  03-C:  Mine-ship  interaction  41 

Project  03-D:  Environmental  impact  of  MCM  on  sonar  design  and  performance  45 

Project  03-E:  High  fidelity,  broadband  structural  acoustic  scattering  models 
for  buried  and  partially  buried  mines  53 

Project  03-F:  Minehunting  sonar  performance  model  59 

Thrust  04  Tactical  Active  Sonar  (TAS) 

Project  04-A:  Advanced  active  sonar  61 

Project  04-B:  Deployable  underwater  surveillance  systems  (DUSS)  75 

Project  04-C:  Low  frequency  shallow  water  reverberation  and  propagation: 
adaptation  to  large  bandwidths  81 

Project  04- D:  Broadband  modelling  of  ASW  monostatic  and  multistatic 
propagation  and  reverberation  87 


Thrust  05  Command  Support  (CDS) 


Project  05-A:  Command  and  NATO  support  95 
Project  05-B:  Planning  support  for  underwater  warfare  operations  101 


Thrust  06  Exploratory  Research  (EXR) 


Project  06-B:  Focused  acoustic  field  (FAF)  studies  103 
Project  06-C:  Sound,  oceanography  and  living  marine  resources  (SOLMR)  109 


Annexes 


A  Ship  Management  Office,  NRV  Alliance  and  T-Boat  Mannings 

B  Annual  bibliography  of  SACLANTCEN  reports  with  abstracts  118 

C  Organization  and  staff  members  with  effect  from  1st  January  2000  129 

D  Visitors  and  meetings  133 

E  Scientific  Committee  of  National  Representatives 
and  National  Liaison  Officers  137 

F  Personnel  by  category  and  nationality  141 


Plus  ga  change ...  1959-1999 


The  year  1999 saw  the  40th  anniversary  of  the  Centre  and  the  50th  of  NATO.  The  wisdom 
and  vision  of  the  founding  fathers”  manifested  by  the  creation  of  an  undersea  research 
centre  for  the  NATO  Alliance  "...  pooling  the  technical  data  ana  scientific  personnel  of 
NATO  countries,  in  an  effort  to  solve  some  of  the  problems  that  none  of  these  countries  had 
been  able  to  solve  independently”  have  allowed  scientists,  engineers,  technicians  and  naval 

fcers  from  national  research  organizations,  universities,  naval  establishments  and  other 
inguished  and  respected  institutes  to  work  at  the  Centre.  All  brought  their  wide  national 
experience  and  knowledge,  a  deep  understanding  of  the  challenge,  enthusiasm  for  scientific 
research  and  most  importantly,  a  shared  belief  in  its  relevance.  Several  hundreds  of  NATO's 
best  brains  in  underwater  acoustics,  oceanography,  signal  processing,  mathematics,  physics 
and  ocean  engineering  have  inspired  each  other  to  excel,  supported  by  the  expertise  and 
ingenuity  of  technical  support  staff  and  the  dedication  of  the  small  but  effective  administrative 
support  department.  They  have  made  the  Centres  relevance  to  NATO  outstanding  and 
unique,  under  the  policy  guidance  ofSACLANT,  augmented  by  invaluable  counsel from  the 
Scientific  Committee  of  National  Representatives. 

Although  the  Kosovo  crisis  prevented  the  Secretary  General,  SACLANT  and  the  Military 
Committee  from  being  present  as  planned,  the  Fortieth  Anniversary  of  the  Centre  was 
appropriately  marked  by  a  number  of  events  attended  by  staff,  VADM  Manlio  Galliccia, 
RADM  Giuseppe  Celeste  and  RADM  Dino  Nascetti,  ( Chairman  of  the  SACLANTCEN 
Scientific  Committee  of  National  Representatives),  representing  the  host  nation  which  makes 
an  inestimable  contribution  to  the  success  of  the  Centre  and  six  former  directors  and  deputy 
directors.  Congratulatory  messages  were  received  from  the  Secretary  General,  Dr.  Javier 
Solana;  SACLANT,  ADM  Harold  Gehman  and  the  Chairman  of  the  NATO  Military 
Committee,  ADM  Guido  Venturoni. 

The  continuing  quest  for  adequate funding  remained  a  dominant  theme.  Numerous  briefings 
to  NATO  HQ  committees  resulted  in  an  increasing  awareness,  that  although  the  problem 
has  not  yet  been  resolved,  the  Centre’s  predicament  is  firmly  on  the  agenda.  There  were  some 
successes,  such  as  funding  for  replacement  of  the  46-year  old  poor  relation  of  the  research 
fleet,  the  coastal  research  vessel  Manning  and  for  an  autonomous  underwater  vehicle.  The 
new  vessel,  provisionally  named  Leonardo,  is  scheduled  to  be  in  service  in  2002,  while  the 
AUV  will  be  delivered  in  2000,  inaugurating  a  new  era  in  the  Centre’s  programme  of  work. 
We  continue  to  plan  for further  updates  in  our  equipment  and  infrastructure,  while  making 
steady,  recognized progress  in  improving  efficiency  commensurate  with  the  nature  of  scientific 
research. 

The  reader  will  find  in  this  report  how  the  Centre  discharges  its  mission  under  the  terms  of 
its  Charter.  The  raison  d’etre  of  the  Centre  as  quoted  earlier  is,  after forty  years,  mirrored  in 
the  recently  approved  new  NATO  research  and  technology  strategy  "  . .  to  ensure  that  the 
Alliance  has  at  its  disposal  the  best  scientific  and  technical  capability  that  member  nations 
are  prepared  to  make  commonly  available”.  In  order  to  ensure  that  eminent  scientists  will 
continue  to  contribute  to  fulfilling  NATO’s  requirements,  it  is  of  utmost  importance  that 
this  strategy  is  translated  into  a  stable  funding  regime.  Only  then  can  NATO,  the  most 
successful  alliance  in  history,  be  guaranteed  that  the  quality  of  our  work  as  exemplified  in 
this  annual  report,  can  continue. 


Jan  L  Spoelstra 
Director 


See  photographs  on  pages  2  and  1 1 7 


The  celebrations  to  commemo¬ 
rate  the  fortieth  anniversary  of 
the  Centre  in  May  1999  were  at¬ 
tended  by,  back  row,  left  to  right: 
Dr  Richard  Nagelhout,  (US), 
Deputy  Director,  1991-1 994;  Dr 
John  H.  Foxwell,  (UK),  Director, 
1990- 1993;  Dr  David  L.  Bradley, 
(US),  Director,  1993-1996; 
RADM  (retd)  Jan  L  Spoelstra, 
(NL),  Director  since  October 
1996;  Prof.  Dr.  Peter  C.  Wille, 
(GE),  Director,  1987-1990;  Ing 
Giancarlo  Vettori,  (IT),  Deputy 
Director,  1994-1998;  Dr  William 
J.  Jobst,  (US),  Deputy  Director 
since  October  1998;  front  row 
seated:  Ms  Heather  Barrett, 
Secretary  to  SACLANTCEN  Di¬ 
rectors,  1959-1999  and  Mr 
Melvin  Arsove,  (US),  Deputy 
Scientific  Director,  1959-1961. 
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RADM  Dino  Nascetti,  Direttore,  Arsenale  Mifitare 
Marittimo  and  Chairman,  SACLANTCEN  Scientific 
Committee  of  National  Representatives,  presenting 
certificates  commemorating  40  years  of  service, 
signed  by  the  Secretary  General,  to  Mrs  MariaPia 
Vergassola  and  Mrs  Gabriella  Parmigiani. 


A  group  photograph  taken  during 
the  annual  maintenance  of  the 
Centre  laboratory  on  the  island  of 
Formica  Grande.  Standing  left  to 
right:  Angelo  Spairani,  Sauro 
Giusti,  Emilio  De  Cola.  Seated  left 
to  right:  Umberto  Varlese,  the 
Director,  Mauro  Lombardi,  Mauro 
Pini,  Umberto  Fabiani. 


Thrust  01  Rapid  Environmental 
Assessment  (REA) 

Project  01 -A:  Rapid  assessment  of  ocean  parameters 
(Alliance  days  -  49,  Manning  days  -  7) 


Operational  Relevance 

Rapid  Response,  the  third  in  the  series  of  demonstrations  of  skills  in  rapid  environmen¬ 
tal  assessment  (REA)  ended  in  1998.  A  comparable  demonstration  bringing  together 
research  organizations  and  military  command  structures  with  common  purpose  will  not 
reoccur  until  2002  in  support  of  Strong  Resolve.  REA  methodologies  continue  to  be 
developed  in  parallel  with  validation  of  the  SACLANT  concept  of  REA  operations  by  the 
NATO  MILOC  community. 


MILOC  support ,  data  fusion  and  transmission 

In  preparation  for  operational  REA  surveys,  data  transfer  between  ships  was  effected  for 
the  first  time  with  spread  spectrum  radios.  Two-dimensional  sound  velocity  fields  on  the 
track  between  an  acoustic  source  and  a  receiving  array,  continuously  acquired  on  a  tow 
ship,  were  accessible  via  the  computer  network  of  NRV  Alliance  by  standard  file  transfer, 
the  high  bandwidth  radio  link  between  the  ships  being  fully  transparent.  Spread  spectrum 
technology  is  easy  to  install  and  to  use.  It  does  not  suffer  from  multipath  signal  cancellations 
when  ships  are  close  together.  After  having  gained  positive  experience  in  another  sea  trial, 
spread  spectrum  technology  was  proposed  by  the  Centre  and  accepted  by  SACLANT  for 
operational  REA  data  communication  between  survey  ships.  It  will  be  in  place  for  Linked 
Seas  2000,  planning  for  which  commenced 
in  January  1999.  During  the  REA  survey, 

SACLANTCEN  personnel  will  embark  on 
one  of  the  survey  ships  and  work  on  bottom 
assessment  for  mine  warfare.  In  addition, 
the  Centre  will  be  engaged  in  ocean 
modelling  and  data  fusion. 


Figure  01-A.1.  Positions  of  CTD  (conductivity, 
temperature,  depth)  stations  (white  circles)  in  the  survey 

regions  of  I  ON  I AN99  and  temperature  distribution  (°C) 
at  10  m  depth  obtained  by  objective  analysis.  Note  the 
cold  water  on  the  Tunisian  Plateau  and  around  the  heel 
of  the  Italian  boot.  Objective  analysis  smooths  fields  and 
occasionally  creates  unrealistic  extrapolations  such  as 
the  discontinued  Adriatic  coastal  current. 
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Oceanographic  field  work 

Of  the  three  research  cruises  with  NRV  Alliance,  only  the  IONIAN99  cruise  in  February  was 
exclusively  dedicated  to  oceanography,  investigating  the  origin  and  generating  mechanism 
of  patches  of  exceptionally  cold  water  observed  at  intermediate  depth  in  the  central  Ionian 
Sea  in  1 998.  Subsurface  lenses  of  deviating  temperature-salinity  characteristics  affect  sonar 
performance  because  of  the  changes  in  the  sound  velocity  profile.  Volumes  of  cold, 
subsurface  water  that  might  be  detected  at  any  season,  can  only  be  generated  in  winter, 
when  heat  loss  of  the  ocean  increases  surface  water  density  to  a  level  equal  to  or  exceeding 
that  of  intermediate  water. 

Winter  surface  temperatures  in  the  vicinity  of  the  Ionian  Sea  observed  by  satellites  are 
lowest  in  the  outflow  from  the  Adriatic  Sea  and  on  the  Tunisian  shelf.  Both  areas  were 
visited  during  IONIAN99.  Positions  of  conductivity,  temperature,  depth  (CTD)  profile 
measurements  are  given  in  Fig.  01 -A.  1  on  surface  temperature  maps  smoothed  by  objective 
analysis. 

The  most  notable  observation  of  the  IONIAN99  cruise  was  the  outflow  of  dense  Adriatic 
Surface  Water  on  the  100  m  deep  Italian  shelf.  In  February  the  density  of  the  coastal 
current  can  even  exceed  that  of  Adriatic  Deep  Water  in  the  Strait  of  Otranto.  The  dense 
water  is  prevented  from  immediately  gliding  down  the  shelf  slope  only  by  Coriolis  force. 
After  streaming  around  the  southern  cape,  it  leaves  the  shelf  and  under  entrainment  of 
surrounding  less  dense  water,  proceeds  to  the  depth  of  neutral  buoyancy.  Figure  01-A.2 
shows  profiles  south  of  Capo  S  Maria  di  Leuca.  Close  to  the  coast,  the  water  is 
homogeneously  cold,  fresh  and  dense.  Further  out,  the  surface  water  continuously 
approaches  Ionian  Surface  Water  conditions.  At  the  bottom,  water  density  decreases.  At 
the  second  position  from  the  south,  the  sea  bottom  is  deeper  than  the  density  horizon  of 
the  modified  coastal  water,  which  in  160  m  depth  exists  as  a  subsurface  patch  replacing 
the  high  salinity,  Levantine  Intermediate  Water  of  the  southernmost  profile.  When  a 
subsurface  patch  spreads  horizontally,  it  gains  anticyclonic  vorticity  by  which  it  is  enabled 
to  survive  over  a  long  period  and  travel  through  the  Ionian  Sea. 


Figure  01-A.2,  Temperature,  salinity,  density  anomaly  and  sound  velocity  profiles  truncated  at 300  m  depth.  The  cold 
and  least  saline  coastal  water  is  denser  than  intermediate  water  further  out.  In  the  temperature-salinity  diagram,  it 
appears  as  a  dot  on  the  same  isopycnal  as  the  Eastern  Mediterranean  Deep  Wa  ter. 
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Figure  01 -A. 3.  Sample  high  resolution  sound  velocity  field  between  the  acoustic  source  and  the  receiving  vertical 
hydrophone  array  Spatial  variability  was  measured  with  the  towed  CTD  chain  during  all  fixed  range  experiments. 


The  data  set  obtained  during  IONIAN99  south  of  35  0  N  covers  a  “white  spot”  in  winter 
climatology.  The  dense  water  produced  on  the  very  shallow  Tunisian  Plateau  by  cooling  and 
evaporation  flows  as  a  bottom  current  towards  the  deeper  Ionian  Sea.  It  is  identified  far 
from  its  origin  on  the  African  continental  slope  before,  by  dilution  with  adjacent  water,  it  is 
totally  merged. 

The  oceanographic  component  of  the  second  sea  trial  ADVENT99  in  May,  supported  the 
acoustic  measurements  by  providing  high  resolution  ocean  stratification.  An  example  of 
the  sound  velocity  field  on  the  main  acoustic  track  is  shown  in  Fig.  01-A.3. 

During  CORFU99,  the  third  cruise  of  the  project,  oceanographic  measurements  contributed 
to  mine  warfare  related  rapid  environmental  assessment.  Prior  to  and  after  the  investigations 
in  the  vicinity  of  the  island,  CTD  stations  were  taken  in  an  area  of  1 15  x  175  km,  which  in 
addition  to  the  Greek  shelf,  contains  deep  water  of  the  northern  Ionian  Sea  and  the  Strait 
of  Otranto.  The  data  were  used  for  initialization  and  verification  of  numerical  ocean  model 
runs  (Fig.  01-A.4). 

The  SKERKI96  data  set,  which  was  collected  in  the  Strait  of  Sardinia  and  eastward  to 
Sicily,  has  been  subjected  to  detailed  investigation  and  prepared  for  publication.  Skerki 

Bank  is  a  barrier  between  the  Sardinian 
channel  and  the  Strait  of  Sicily,  which 
impedes  the  flow  of  intermediate  water. 
The  choice  of  the  area  was  motivated  by 
the  fact  that  knowledge  of  this  region  is 
essential  for  understanding  the  exchange 
between  western  and  eastern 
Mediterranean.  The  high  quality  of  data 
is  ideal  for  use  in  numerical  ocean  models 
and  model  validation.  A  sample  of  the 
analysis  is  displayed  in  Fig.  01 -A. 5, 
showing  the  colour  coded  distribution  of 


Figure  01-A.4.  Forecast  of  surface  temperature  for 
1 1  September  1999, 12  days  after  model  setup  during 
CORFU 99.  The  temperature  ranges  from  16  to  28 

0  C.  Upwelling  of  cold  water  off  the  coasts  of  Albania 
and  Corfu  is  a  consequence  of  atmospheric  forcing. 
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Figure  01-A.5.  Salinity 
(colour)  at  the  salinity 
maximum  and  geostrophic 
currents  (vectors)  vertically 
averaged  between  100  m  and 
700  m  depth.  Currents  were 
calculated  from  horizontal 
pressure  differences  which 
are  functions  of  the 
temperature  and  salinity 
distribution.  Areas  left  white 
indicate  that  no  maximum 
was  found  or  that  the 
statistical  error  of  the  analysis 
exceeds  a  certain  threshold. 


salinity  at  the  depth  of  its  maximum  and  the  currents  at  the  same  level.  The  salinity  maximum 
is  associated  with  Levantine  Intermediate  Water  that  is  formed  in  the  eastern  Mediterranean 
in  winter  and  mostly  occupies  the  depth  range  from  200  to  500  m  of  the  entire  Mediterranean. 


Numerical  modelling 

The  Harvard  Ocean  Prediction  System  (HOPS)  is  the  primary  tool  for  ocean  modelling  at 
SACLANTCEN.  Its  forecast  kernel  is  based  on  the  system  of  hydrodynamic  “primitive” 
equations  with  bottom  following  (sigma)  coordinates.  Packages  for  objective  analysis,  data 
assimilation  and  additional  capabilities  such  as  a  coupled  model  for  biological  activity  are 
components  of  the  system.  HOPS  is  applied  operationally  for  rapid  environmental 
assessment  and  used  for  oceanographic  research  in  order  to  analyze  an  observed  situation 
and  to  forecast  its  evolution.  Figure  01-A.6  is  a  10-day  forecast  of  surface  currents  and 
salinity  based  on  IONIAN99  data  from  the  Tunisian  shelf  and  the  Sicilian  Channel.  The 
jet-like  inflow  of  low  salinity  Atlantic  Water  through  the  Strait  of  Sicily  from  the  west  is  clearly 
visible.  During  CORFU99,  HOPS  was  applied  to  a  shelf  break  and  large  coastal  area 
containing  islands.  The  HOPS  prediction  of  surface  temperatures  in  Fig.  01 -A. 4  shows 
upwelling  areas  off  the  Albanian  and  Corfu  coast  that  are  known  from  satellite  images. 

Numerical  ocean  models,  which  integrate  forward  in  time  the  equations  of  ocean  dynamics 
are  the  most  intuitive  tools  for  ocean  prediction.  They  provide  a  complete  physical  view  of 
the  complexities  of  ocean  evolution.  On  the  other  hand,  numerical  forecasts  of  the  ocean 
require  profound  understanding  of  the  dynamical  laws  controlling  ocean  processes. 
Consequential  to  the  wide  range  of  space-time  scales  arising  from  the  turbulent  nature  of 
ocean  dynamics,  small  scale  oceanic  processes  cannot  be  resolved  in  coarse  schemes, 
but  have  to  be  parameterized  according  to  numerical  mesh  width. 

A  numerical  and  theoretical  study  of  the  influence  of  small  scale  variability  on  large  scale 
currents  has  been  carried  out  with  oceanographers  from  the  Universidad  de  le  Islas  Balears. 
It  has  been  found  that  small  scale  eddies  may  reverse  the  direction  of  large  scale  currents 
when  they  interact  with  the  underlying  bottom  topography. 

An  alternative  approach  to  forecasting  ocean  evolution  consists  of  extracting  dynamical 
information  from  empirical  data  without  imposing  an  explicit  dynamical  model. The  extracted 
information  about  the  past  of  the  system  is  used  to  predict  future  evolution. The  problem  of 
ocean  forecasting  is  thus  reduced  to  finding  the  dynamical  model  which  best  explains  the 
features  of  a  given  time  series  of  observations. 
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Satellite  imagery  is  appropriate  to  continuous  monitoring  of  the  space-time  variability  of  the 
ocean.  Although  satellites  provide  only  a  limited  view  of  the  state  of  the  ocean,  the  information 
may  be  sufficient  for  operational  purposes.  Promising  attempts  have  been  made  to  develop 
an  ocean  forecasting  system  using  satellite  imagery.  The  analysis  of  ocean  dynamical 
tendencies  in  time  series  of  satellite  data  and  forecasting  of  future  states  of  the  ocean  is 
based  on  genetic  algorithms.  Figure  01-A.7  shows  the  mean  SST  pattern  observed  in  the 
Alboran  Sea  in  November  1998  and  the  SST  pattern  forecast  one  month  in  advance  by  this 
approach. 

Satellite  remote  sensing 

Images  of  sea  surface  temperature 
provided  by  NOAA  satellites  have  been 
processed  for  many  years  and  used  in 
many  projects  with  great  benefit.  Since 
spring  1 999,  when  a  new  antenna  (1 .5  0) 
was  installed  and  a  contract  signed  with 
NASA  as  an  authorized  SeaWiFS 
research  station,  SACLANTCEN  is  also 
receiving  colour  images  of  the  sea  surface 
that  can  be  converted  into  quantitative 
information  about  transparency  and 
chlorophyll.  Raw  data  are  preprocessed 
some  hours  after  arrival  and  delivered  to 
NASA.  A  key  required  for  the  conversion 
of  preprocessed  data  into  physical  units 
is  given  to  the  Centre  a  few  weeks  later, 
in  exceptional  cases  for  real-time 
conversion.  Ocean  colour  images 
processed  at  SACLANTCEN  have  been 
used  in  Project  06-C. 

Instruments  for  rapid  oceanographic  assessment 

The  prototype  of  a  Shallow  Water  Expendable  Environmental  Profiler  (SWEEP)  was 
completed  in  1998.  The  final  design  refinement  prior  to  volume  production  was  delayed  in 
order  to  allow  development  of  a  Shallow-water  Environmental  Profiler  in  Trawl-safe,  Real¬ 
time  configuration  (SEPTR).  The  SEPTR  system  combines  the  advantages  of  the  trawl- 
safe  bottom  instrument  platform  and  the  SWEEP  water  column  profiler  developments.  It  is 
intended  for  extended  duration  deployments  with  real-time  data  return  and  control  via  two- 
way  cellular  or  satellite  communication  in  areas  where  water  column  instruments  are  at  risk 
from  fishing  trawlers  (Fig.  01-A.8). 

The  SEPTR  design  consists  of  a  trawl-safe  bottom  platform,  which  houses  an  Acoustic 
Doppler  Current  Profiler  (ADCP),  wave/tide  gauge,  ambient  noise  sensor  array  and  buoyant 
water  column  profiler. The  profiler  buoy  and  associated  winch  are  designed  for  autonomous 
vertical  profiling  of  CTD  and  optical  properties  of  the  water  column  down  to  1 00  m.  Features 
include  DGPS  navigation,  two-way  communication  when  surfaced,  acceleration  and  magnetic 
field  sensors.  The  bottom  platform  includes  an  extended  duration  battery  package  in  the 
recoverable  barnacle-shaped  housing.  For  recovery,  a  messenger  buoy  is  released  by  a 
command  through  the  normal  communication  channel  to  the  profiler  or,  acoustically,  to  the 
bottom  mooring.  If  the  platform  becomes  inverted,  the  system,  which  after  release  of  the 
ballast  is  buoyant,  can  still  be  recovered  by  command  to  a  backup  acoustic  transponder 
system  in  the  bottom  platform. 


37° 
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35° 


34° 


33° 
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Figure  01-A.6.  Forecast  of  surface  currents  and  salinity 
(colours)  for  26  February  1999.  The  forecast  is  based  on 
data  collected  from  10  to  16  February  during  IONIAN99. 
The  scale  arrow  in  the  upper  left  corner  corresponds  to  a 
speed  of  20  cm/s.  Note  the  highly  variable  currents  in  the 
Strait  of  Sicily  and  the  quiet  areas  on  the  Tunisian  shelf. 
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Figure  01-A.7.  a)  Surface  temperature  of  theAlboran  Sea  (western  Mediterranean)  obtained  as  an  average  of  satellite 
observations  in  November  1998.  b)  Sea  surface  temperature  one  month  later  predicted  from  the  time  behaviour  of 
orthogonal  functions  determined  by  a  genetic  algorithm. 


The  system  is  designed  for  360  profiles  during  a  3-month  deployment  in  100  m  water 
depth,  with  current,  wave/tide  and  noise  measurements  every  hour. Two  way  communication 
of  data,  position  and  control  allows  profile  results  to  be  returned  in  real  time  and  operational 
commands,  profile  schedules  and  DGPS  correctors  to  be  sent  to  multiple  profiler  instruments. 
Complete  profiles  will  also  be  stored  on  board  the  recoverable  units.  During  sea  trials  in 
November  1999,  critical  subsystems  were  successfully  tested. 

The  real-time  data  transmission  capability  of  SEPTR  is  well-suited  to  provision  of  water 
property  and  ocean  current  data  for  assimilation  into  numerical  ocean  models.  Two  units 
will  be  operationally  tested  during  the  sea  trial  GOATS  2000,  which  in  addition  to  bottom 
acoustics  and  mine  warfare  related  tasks,  encompasses  an  oceanographic  observation 
and  modelling  programme. 


Figure  01-A.8.  The  Environmental 
Profiler  in  Trawl-safe,  Real-time 
configuration  (SEPTR)  is  a  high 
value  multi-sensor  platform  for  REA 
in  accessible  areas.  After 
deployment,  it  does  not  interfere  with 
fishing  activity.  At  preprogrammed 
times,  a  float  rises  for  measurements 
in  the  water  column  and  data 
delivery. 


Comprehensive  experiments  and  REA  methodology  validation 

REA  system  components  such  as  oceanographic  and  acoustic  instruments,  remotely  sensed 
data,  real-time  ocean  modelling,  inverse  acoustic  methods,  communications  and  data  fusion 
are  developed  as  the  responsibility  of  the  respective  task  leaders.  When  several  components 
or  the  entire  REA  system  is  tested  or  a  specific  REA  survey  is  requested,  trials  organization 
is  transformed  into  a  complex  task,  as  during  the  ADVENT99  cruise  in  May  1999,  which 
combined  REA  tasks  with  the  acoustic  REA  Project  01 -B  and  a  multinational  joint  research 
programme.  Figures  01-A.9,  01 -A.  10  and  01-A.3  illustrate  data  communications  and  high 
resolution  oceanography. 
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Figure  01 -A. 9.  Signals  are 
transmitted  and  received  over  a 
fixed  range.  Environmental 
conditions  on  the  acoustic  range 
are  measured  at  the  same  time. 


Following  a  Greek  proposal  for  a  limited  REA  experiment  in  their  territorial  waters  in  order 
to  demonstrate  skills  acquired  during  Rapid  Response,  the  mine  warfare  related  CORFU99 
sea  trial  was  conducted  in  September  1999  with  the  Hydrographic  Service  of  the  Hellenic 
Navy  and  HNS  Pytheas  in  conjunction  with  Project  03-D,  which  specializes  in  environmental 
characteristics  for  mine  warfare.  Figure  01-A.11,  shows  bottom  mapping  from  a  shallow 
water,  multi-beam  echo  sounder,  installed  especially  for  this  cruise. 


Geophysical,  oceanographic  and  acoustic  results  of  the  area  assessment  were  processed 
on  board,  assembled  into  a  data  server  structure  and  transferred  to  a  CD-ROM,  which  was 
sent  to  the  Hellenic  Navy  immediately  after  the  cruise. 


1 


Figure  01-A.10.  CTD  chain  towed  by  the 
Italian  navy  tug  Ciclope  along  the  acoustic 
transmission  track.  The  two-dimensional 
data  of  conductivity,  temperature  and  depth 
were  transmitted  to  NRV  Alliance  in  real¬ 
time  via  spread-spectrum  radio. 


Figure  01-A.11.  Very  high  resolution 
bathymetry  in  a  test  area  northwest  of  Corfu 
obtained  from  a  Simrad  EM-3000  echo 
sounder  deployed  from  Alliance.  The  size  of 
each  box  is  713  by  926  m.  For  structure 
enhancement,  a  hypothetical  light  source  was 
applied. 
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Grab  sampling  and  the  Greek  vessel 
Pytheas  during  CORFU99 
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Project  01 -B:  Rapid  environmental  assessment  of 
operational  acoustic  parameters 
(Alliance  days  -  25,  Manning  days  - 12) 


Operational  relevance 

Remote  and  rapid  assessment  of  unknown  littoral  environments  provides  the  necessary 
acoustic  input  parameters  to  enhance  performance  prediction  of  sonar  platforms. 


Seabed  geoacoustic  inversion  techniques 

Seabed  geoacoustic  inversion  techniques  have  been  developed  and  validated  to  rapidly 
and  accurately  predict  the  range  dependent  geoacoustic  properties  of  littoral  areas.  These 
techniques  will  eliminate  or  reduce  the  need  for  seabed  coring. 

The  strong  dependence  of  shallow  water  sound  propagation  on  sea-bed  type  has  led  to 
the  development  of  inversion  methods,  which  use  measured  acoustic  transmissions  to 
determine  bottom  properties  such  as  sound  speed,  density  and  attenuation  constant.  Using 
acoustics  for  seabed  characterization  has  several  advantages  over  other  methods  (e.g. 
coring). The  measurement  equipment  is  easily  and  quickly  deployed  and  because  acoustic 
signals  are  used,  the  method  is  well  suited  to  characterize  the  bottom  properties  which 
have  the  greatest  influence  on  propagation.  The  inversion  method  is  based  on  Matched 
Field  Processing  (MFP),  the  starting  point  for  which,  is  to  parameterize  the  environment 
with  an  assumed  geoacoustic  model  consisting  of  several  unknown  parameters  (Fig.  01- 
B.1).  Computer  simulation  models  the  acoustic  response  to  different  seabed  types  and 
efficient  search  algorithms  identify  the  environment  which  matches  modelled  and  measured 
data. 

To  illustrate  (and  validate)  the  utility  of  MFP  inversion,  two  sites  were  chosen  which  have 
very  different  bottom  types.  The  first  data  set  was  collected  during  the  EnVerse97  experiments 
on  the  Adventure  Bank  off  Sicily,  where  the  seabed  is  of  sand  over  a  rock  sub-bottom. The 
second  data  set  was  taken  from  the  PROSIM97  experiments  near  Elba,  where  the  seabed 


Figure  01  .B.1  Left  panel  Is  a  representation  of  the  true  environment  and  the  right  panel  is  the  geoacoustic  model 
consisting  of  several  seabed  parameters.  There  is  a  sediment  layer  of  thickness  “h"  over  an  infinite  (homogeneous) 
sub-bottom.  The  sediment  layer  has  a  sound  speed  gradient  denoted  as  “csed ,  ’’and  “csed?  "at  the  upper  and  lower  points 
in  the  sediment,  the  sub-bottom  has  sound  speed  “cMOT”.  Attenuation  “a”  and  density  “p”  are  also  included  in  the 
geoacoustic  model. 
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consists  of  an  acoustically  slow  clay  layer  over  silt  sub-bottom.  At  both  sites,  towed  sources 
transmitted  broadband  signals  from  200-800  Hz  which  were  received  on  a  vertical 
hydrophone  array.The  data  sets  were  inverted  using  identical  procedures  and  environmental 
parameterization  (Fig.  01-B.1).  Figure  01-B.2  shows  the  best  fit  between  the  computer 
simulated  and  measured  data  for  the  EnVerse97  site.  The  inverted  values  for  the  bottom 
are  summarized  in  Table  01-B.1.  The  slow  clay  layer  over  silt  (PROSIM97  site)  and  the 
faster  sand  bottom  over  rock  (EnVerse97  site)  were  correctly  identified  from  the  inversions. 
The  significance  of  this  work  is  to  illustrate  correct  seabed  characterization  without  assuming 
prior  environmental  knowledge  using  an  unbiased,  automatic  inversion  procedure. 


Figure  01  -B.2  Modelled  (a)  and  measured  (b)  acoustic  data  taken  along  a  62  m  vertical  array  from  the  EnVerse97 
experiments.  Measured  and  simulated  pressure  magnitudes  are  shown  on  a  log  scale  as  a  function  of  depth  normalized 
at  each  frequency  in  the  band  between  200-800  Hz. 


The  importance  of  correctly  identifying  the  seabed  type  can  be  illustrated  by  example.  For 
very  different  bottom  types  such  as  for  the  sites  in  Table  01-B.  1,  large  differences  in 
transmission  loss  (TL)  are  expected.  But  how  doTL  predicitions  differ  when  the  environmental 
properties  are  taken  from  an  archive  such  as  the  Allied  Environmental  Support  system 
(AESS)  compared  to  sea-bed  properties  from  MFP  geoacoustic  inversion?  Broadband 
acoustic  measurements  with  2  km  source-receiver  separation  were  taken  during  the 
Advent99  experiments  which  took  place  at  a  shallow  water  site  (80  m  depth)  in  the  Strait  of 
Sicily.  These  data  were  used  for  geoacoustic  inversion  to  characterize  the  seabed.  The  TL 
predictions  as  a  function  of  range,  using  MFP  inverted  values  for  the  seabed  and  the 
archived  values  are  shown  in  Figure  01-B.3.  Also  shown  are  TL  measurements  (which 
were  not  used  in  the  inversion)  taken  at  2,  5  and  10  km,  which  are  in  agreement  with  the 
predictions  using  MFP  inverted  seabed  values.  The  differences  in  the  TL  predictions  are 
large.  Using  the  archived  seabed  properties  would  lead  to  over  predicting  the  range  of 
target  detection  by  15-20  km.  In  Fig.  01-B.4,  measured  TL  at  10  km  range  is  compared  with 
predicted  TL  using  the  MFP  inverted  values  for  the  seabed  with  those  taken  from  the 
archive.  The  TL  level  and  depth  structure  is  matched  well  using  the  MFP  inverted  seabed 
properties,  but  not  when  archived  values  are  applied. 


Site 

Csed,1 

(m/s) 

^sod,2 

(m/s) 

h 

(m) 

(dB  /) 

(g/cm3) 

£ 

^bottom 

(m/s) 

PR097 

1464 

1496 

4.5 

0.03 

1.6 

1549 

ENV97 

1545 

1559 

9.4 

0.73 

1.5 

1739 

Table  01-B.1  Inversion  results  for  EnVerse97  (ENV97)  and  PROSIM97  (PR097)  experimental  data.  Both  data  sets 
were  collected  on  a  vertical  array  using  a  towed  source  near  mid-water  depth  in  the  frequency  band  200-800  Hz. 
Using  the  same  data  type  and  inversion  procedure,  the  two  different  bottom  types  were  correctly  identified.  The 
PR097  inverted  sediment  and  sub-bottom  are  in  agreement  with  the  expected  seabed  type  of  clay  over  silt  and  the 
ENV97  inverted  values  are  consistent  with  sand  over  rock. 
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In  Fig.  01-B.5,  matched  field  beamforming  is  used  to  localize  the  acoustic  source  (located 
near  2  and  10  km  from  the  receiver  array). The  figure  shows  good  localization  of  the  source 
in  range  and  depth.  Rather  than  using  plane  waves  as  with  conventional  beamforming,  the 
matched-field  beamformed  image  uses  the  computer  model  to  generate  wavefronts  more 
appropriate  for  the  environment.  Here,  the  computer  model  uses  the  environment  taken 
from  the  MFP  inversion. 

The  research  demonstrates  that  the  geoacoustic  inversion  procedure  can  be  successfully 
applied  without  prior  knowledge  of  the  seabed.  The  Advent99  data  and  analysis  showed 
that  data  acquired  at  relatively  short  source-receiver  range  (2  km)  give  a  stable  prediction 
of  seabed  properties  over  many  hours.  These  seabed  properties  were  used  to  accurately 
predict  TL  measured  on  a  different  day  with  10  km  source-receiver  range  separation.  The 
Advent99  data  and  analysis  also  show  the  value  of  using  correct  seabed  properties  for  TL 
prediction  and  target  localization  in  range  and  depth. 


Figure  01 -B. 3  Green  curve  shows 
predicted  TL  using  archived  seabed 
parameters,  black  curve  uses  seabed 
parameters  determined  by  MFP 
geoacoustic  inversion.  The  red  diamonds 
show  measured  TL  (The  data  shown  as  red 
diamonds  was  not  used  in  the  geoacoustic 
inversion).  The  TL  is  depth  and  frequency 
averaged  using  tones  spaced  at  100  Hz 
from  200-700  Hz. 
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Figure  01 -B. 4  Red  curves  show  measured  TL  as  a  function  of  depth  at  10  km  source  receiver  separation  at  the 
indicated  frequencies,  (water  depth  80m,  source  depth  76m).  Black  curve  is  the  predicted  TL  using  seabed  properties 
taken  from  MFP  inversion  (MFP  inversion  data  was  collected  at  2  km  source  receiver  separation).  Green  curves  were 
generated  using  archived  seabed  parameters. 
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Figure  01 -B. 5  Matched-field 
beamformed  images  of  acoustic 
source  at 2  km  and  10  km  range. 
To  generate  the  image,  model 
based  processing  was  used  with 
the  inverted  seabed  parameters 
from  data  collected  at  2  km.  On 
the  right  are  zoomed  pictures  of 
the  source  location  for  each  multi- 
tone  transmission  (200-700  Hz  in 
100  Hz  increments)  taken  at  4 
min  intervals  over  6  h  (2  km  data) 
and  18  h  (10  km  data).  Note  the 
overall  stability  of  the  source 
position,  with  many  data  points  in 
the  same  location.  There  is  a 
larger  spread  of  predicted  source 
positions  at  10  km  due  to  time- 
variable  ocean  processes. 
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Geoacoustic  inversion  using  drifting  sonobuoys 

Drifting  sonobuoys  which  feature  large  dynamic  range  signal  conditioning  and  precise  GPS 
positioning,  were  utilized  in  the  rapid  assessment  of  littoral  environments  to  resolve  the 
range-dependence  of  seabed  properties.  The  buoys  drift  away  from  a  stationary  source, 
receiving  pings  at  known  positions.  Single  hydrophone,  matched  waveform,  inversion 
techniques  are  applied  to  the  received  signals  to  invert  for  the  geoacoustic  parameters. 

A  field  of  eight  prototype  buoys  was  successfully  tested  on  the  western  Sicilian  shelf  in 
EnVerse97.  A  follow-up  geophysical  survey  in  1998  provided  ground  truth  seabed  data  to 
support  and  validate  the  buoy  field  inversion.  A  unique  feature  of  the  survey  is  that  the 
seabed  was  characterized  at  particularly  high  spatial  resolution  using  Swath  multibeam 
system,  seismic  boomer,  in  situ  sound  speed  bottom  probes  (Forchetta)  and  a  large  number 
of  water-sediment  cores.  Figure  01-B.6  shows  several  layers  of  a  GIS  map.  Of  particular 
interest  is  the  thickness  distribution  of  the  sand  cover,  automatically  reconstructed  from  a 
fine  grid  of  digitized  seismic  profiles.  The  locations  of  water-sediment  cores  and  bottom- 
penetrated  acoustic  probes  are  shown.  CD-ROM  versions  with  associated  data  such  as 
seismic  reflection  or  sediment  sound  speed  profiles  are  available.  The  geometry  of  an 
EnVerse  97  acoustic  run  is  overlaid,  including  the  positions  of  a  bottom-moored  sound 
source  and  the  tracks  of  the  current-driven  acoustic  buoys.  The  zoomed  regions  illustrate 
the  degree  of  detail  achieved  in  producing  this  unique  two  dimensional,  ground  truth, 
geoacoustic  model. 


Satellite  remote  sensing  and  ambient  noise  estimation 

Algorithms  have  been  developed  for  extracting  estimates  of  shipping  density  and  wind 
speed  from  satellite  remote  sensors  (Synthetic  Aperture  Radar).  Shipping  noise  and  wind 
speed  are  two  environmental  parameters  which  are  required  in  directional  ambient  noise 
prediction  models  such  as  RANDI. 


—  16  — 


DePMm>  Depth  (m) 


The  physical  processes  contributing  to  ambient  noise  variability  include  the  effects  of  surface 
wind  and  wave  activity,  rainfall,  wave-wave  interactions,  monomolecular  surface  films,  seismo- 
acoustic  noise,  marine  organisms  and  anthropogenic  noise.  Here  we  demonstrate  the 
methodologies  for  an  ocean  acoustic  ambient  noise  prediction  system  from  satellite  derived 
information. 

The  surface  processes  which  dominate  ambient 
noise  between  500  to  50,000  Hz  can  be 
measured  from  a  number  of  satellites,  including 
the  special  sensor  microwave  imager  (SSM/I) 
on  the  defense  meteorological  satellite  program 
(DMSP),  an  ideal  sensor  for  wide-area,  daily 
coverage  of  critical  environmental  parameters. 

SSM/I  is  passive  microwave  with  a  swath  width 
of  1400  km  and  spatial  resolution  of  50  km.  In 
Fig.  01  -B.7  the  surface  wind  speed  is  estimated 
by  measuring  sea  surface  microwave  emissions 
which  are  related  to  wave  structure  and  foam, 
which  in  turn  are  related  to  the  surface  winds. 

The  wind  retrieval  method  uses  a  linear 
combination  of  measured  SSM/I  brightness 
temperatures  that  vary  with  frequency  and 
polarization.  The  SSM/I  accuracy  specification 
for  wind  speed  retrieval  is  ±2  m/s  over  the  range 
of  3-25  m/s.  Rainfall  rates  can  also  be  measured 
from  SSM/I  using  attenuation-based  and 
scattering-based  inversion  techniques.  More 
recent  inversion  algorithms  use  regression 
analysis  based  on  ground-based  measurements 
of  rainfall  and  brightness  temperatures  from  the 
SSM/I  channels. 

For  resolving  the  spatial  variability  of  the  wind  field  in  the  coastal  regions  at  spatial  resolutions 
of  5  km  and  below,  we  utilize  the  ERS-2  synthetic  aperture  radar  (SAR),  an  active  microwave, 
all-weather  sensor.  The  oceanic  wind  vector  is  retrieved  by  inverting  a  scatterometer  model 
function  which  relies  on  the  “Bragg  resonant"  reflections  from  the  slopes  of  wind-generated 
gravity-capillary  short  waves.  In  Fig.  01-B.8  the  CMOD4  scatterometer  model  function  is 
used  to  convert  between  the  SAR  image  intensity  or  the  normalized  radar  cross-section 
and  the  wind  speed.  Because  SAR  responds  only  to  the  line-of-sight  component  of  the 
wind  vector,  the  true  wind  direction  must  be  measured  independently  or  derived  from  image 
processing  techniques.  Here  we  use  the  two-dimensional  image  spectra  of  the  wind-wave 
field  as  a  means  of  deriving  the  wind  direction.  An  empirical  model  is  used  to  convert  wind 
speeds  to  wind-induced  source  levels  at  the  acoustic  frequency  of  300  Hz. 

In  the  frequency  band  between  20  to  500  Hz,  the  noise  spectrum  is  dominated  by  shipping. 
To  model  ship-generated  noise,  shipping  density  must  be  determined  as  a  function  of 
range  and  azimuth  from  the  receiver.  Figure  01-B.9  shows  the  processing  steps  involved  in 


Figure  01 -B. 7  Wind  speed  (m/s)  and  rainfall 
(mm/hr)  and  estimated  noise  levels. 


Figure  01  -B.6  (pages  1 8  and  1 9)  GIS  Map  with  overlays  showing  the  area  and  thickness  of  the  sand  cover,  acoustic 
tracks  geometry  and  positions  of  cores  and  bottom-penetrating  acoustic  probes 
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deriving  shipping  information.  Artificial  neural  network  architectures  for  hard  target  detection 
in  spaceborne  SAR  imagery,  produce  a  return  which  is  a  function  of  radar  frequency,  incident 
angle,  sea  state,  wind  speed  and  look  direction.  In  automatic  detection  a  2-D  SAR  image 
is  classified  using  four  1  -D  profiles.  For  each  profile,  the  network  outputs  a  decision  that  the 
detected  profile  belongs  to  a  predefined  shape.  The  individual  opinions  are  fused  to  derive 
the  final  classification  results.  The  second  and  third  steps  are  measurements  of  the  ship 
heading  and  speed  which  are  critical  to  directional  ambient  noise  calculations.  For  ship 
heading  estimation,  we  rely  on  the  visibility  of  the  wake  patterns,  which  depends  on  the 
choice  of  radar  parameters  and  environmental  conditions.  Detecting  ship  wakes  in  a  SAR 
image  automatically  involves  several  steps  which  include  de-speckling,  edge  detection, 


Figure  01-B.9  Ambient  noise  level  estimation  from  SAR-derived  ship  distributions. 
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transformation  of  the  edge-detected  image  to  the  Hough  parameter  space,  local  maxima 
selection,  inverse  Hough  transform  and  line  detection.  The  Hough  transform  algorithm  is  a 
robust  method  for  automatic  detection  of  lineaments  in  imagery.  For  estimating  ship  speeds 
we  rely  on  the  Doppler-offset  or  the  displacement  of  the  wake  from  the  hard  target  in  the 
image.  In  the  final  step,  source  levels  associated  with  different  ship  types  and  sizes  are 
estimated  using  empirical  data. 


Shallow  water  expendable  environmental  profiler  (SWEEP)  and 
Shallow-water  environmental  profiler  in  trawl-safe ,  real-time 
configuration  (SEPTR)  REA  buoy  development 

The  SEPTR  provides  additional  support  to  environmental  modelling  and  rapid  environmental 
assessment,  permitting  real  time  reporting  of  ADCP,  tide/wave  gauge  and  directional  acoustic 
data,  in  addition  to  SWEEP  water 
column  profiler  data,  over  a  three 
month  period.  Dual  development  of 
SEPTR  and  SWEEP  is  made  practical 
by  designing  common  electronic  and 
sensor  systems.  Sensor  modules  for 
the  first  SWEEP  prototype  were  built 
and  tested  with  hardware  and  software 
for  acoustic  ambient  noise 
measurements.  Submerged  field  tests 
of  the  SWEEP  prototype  and  an  early 
SEPTR  prototype  were  conducted. 

Hardware  and  software  development 
were  commercially  subcontracted. 

Development  of  a  multi-element 
acoustic  subsystem  for  SEPTR 
acoustic  noise  directionality  and  sound 
propagation  experiments  was  initiated 
for  use  on  SWEEP  in  single  element  form.  Figure  01-B.10  shows  the  SWEEP  prototype 
buoy  on  the  bottom,  prior  to  automatic  ballast  release.  Divers  were  used  to  observe  and 
film  buoy  dynamics. 

Water  column  environmental  sensors  were  incorporated  in  a  shallow  water  environmental 
profiler  in  a  trawl-safe  and  real  time  configuration.  The  sensor  system  intended  for  long 
duration  deployments,  with  real  time  data  return,  can  be  controlled  via  cellular  phone  or 
satellite  transmission.  An  expendable,  non-trawl  safe  version  is  under  development. 

Development  tests  included: 


Orbcomm  Low  Earth  Orbit  Satellite  communicator 

GSM  cellular  phone  communicator  and  antenna 

DGPS  navigation  performance  testing  using  GSM  correctors 

CTD  calibration 

Battery  system 

Accelerometer  sensor 

Acoustic  subsystem 

Pressure  case  overpressure  and  motor  seals 
Motor  and  winch  performance 
Submerged  field 


Figure  01-B.10  SWEEP  prototype  prior  to  automatic 
anchor  release. 
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Standard  SWEEP  and  trawler  safe  SEPTR  profilers  are  being  developed.  By  the  fourth 
quarter  2000  the  pre-production  SWEEP  prototype  should  be  complete  and  commercial 
production  undertaken  of  a  series  of  SWEEP  buoys  for  rapid  environmental  assessment 
demonstrations  in  2001  and  2002.  By  September  2000,  two  SEPTR  units  will  be  complete 
and  field  trials  undertaken  as  part  of  environmental  modelling  programmes. 

Autonomous  underwater  vehicle  (AUV)  for  rapid  environmental 
assessment 

The  Centre  specified  and  procured,  through  international  competitive  bidding,  an  AUV  to 
conduct  experiments  in  support  of  REA  and  MCM  studies.  The  AUV  selected  is  a  modified 
Ocean  Explorer  (OEX),  manufactured  by  Florida  Atlantic  University  (Fig.  01-B.11).  OEX  is 
a  modular  vehicle  which  features  a  PC-104  Pentium,  executing  the  supervisory  control 
structure,  interfaced  to  a  distributed  communication  and  control  network,  based  on  the 
open  LONTalk  network  protocol. The  modular  design  results  in  a  field  reconfigurable  vehicle 
which  is  well  suited  to  experimental  work 

The  aft  section  contains  the  propulsion  unit,  control  surface  motors,  navigation  sensors, 
computer,  RF  and  acoustic  communication  systems  and  the  nickel  metal  hydride  (Ni-MH) 
battery  packs.  The  navigation  relies  on  a  standard  suite  of  sensors  (tri-axial  flux  gate 
compass,  precision  rate  gyros  and  accelerometers  and  a  Doppler  Velocity  Logger),  which 
provide  relative  and  geodetic  positional  information.  The  modular  design  facilitates  the 
inclusion  of  a  variety  of  geodetic  navigational  instruments  including  a  DGPS,  long  baseline 
(LBL)  and  ultra-short  base  line  (USBL)  acoustic  tracking  systems. 

The  payload  section  which  contains  the  sensors  is  attached  to  the  aft  portion  by  a  bayonet 
mechanism.  Elements  in  the  payload  are  simply  additional  nodes  on  the  LONTalk  network 
and  are  integrated  with  the  tail  section  with  a  power  conductor  and  a  control  line. 

The  nose  section  contains  a  video  camera,  strobe,  emergency  drop  weight  and  the  vehicle 
recovery  line  device.  OEX  can  be  tracked  acoustically  via  USBL  and  communicate  via  the 
acoustic  modem.  Launching  and  recovery  can  be  achieved  in  seas  up  to  sea  state  4  without 
diver  intervention. 

The  REA  payload  section  includes  a  Sea  Scan  PC  (SSPC)  side  scan  sonar  system, 
developed  by  Marine  Sonic  Technology  Ltd.  The  system  transmits  broad  band  pulses  (5-7 
cycles  in  length)  at  150  and  600  kHz  and  receives  the  scattered  signals  with  very  low  noise 
pre-amplification  circuitry.  The  SSPC  software  allows  automatic  sonar  data  acquisition.  Any 
item  in  the  sonar  imagery  can  be  correlated  with  geographical  position. The  software  provides 
complete  playback  of  sonar  images  acquired.  The  geo-referenced  sonar  data  files,  rectified 
for  orientation  and  track  curvature,  may  be  imported  into  a  GIS  platform  such  as  ArcView. 
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AUV  Parameter 


Shape 


Diameter 


Vehicle  Specification 


Gertler  Series  58  Model  5154E 


0.533  m 


2.66  m  -  3.91  m 


Fig  ure  01  -B .  1 1  Detail  Of  the  Ocean  Explorer  AUV. 
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Thrust  03  Mine  Countermeasures  (MCM) 

Project  03-A:  Detection  and  classification  of  buried 
mines  (Manning  days  13) 


Operational  relevance 

A  mine  may  become  buried  in  marine  sediment  as  a  result  of  impact  in  soft  muddy 
bottoms,  scouring  and  ridge  migration  in  hard  sandy  bottoms  or  the  mine  may  be  designed 
to  be  self  burying.  /As  the  major  mine  influences,  (magnetic,  low  frequency  acoustic, 
pressure)  and  their  explosive  potential  are  unaffected  by  burial,  buried  mine 
countermeasures  area  widely  recognized  NATO  capability  shortfall,  identified  in  the  NATO 
study  MO  2015,  to  which  this  project  has  contributed. 


Mine  sweeping  is  the  only  active 
measure  against  all  types  of  buried 
mines,  but  arming  delays,  ship 
counts  and  pressure  sensors 
reduce  its  effectiveness. 

Conventional  mine  hunting  sonar  is 
ineffective  against  mines  which  are 
completely  buried  in  sediment,  in 
which  sound  absorption  is  very  high 
because  at  the  usual  minimum 
operating  frequency  of  80  kHz, 
sediment  sound  absorption  reduces 
the  level  of  the  target  echo  to  well 
below  that  of  the  seabed. 

Buried  mine  sonar  must  therefore  be  designed  to  minimize  sediment  sound  absorption  and 
the  loss  of  resolution  which  affects  sonar  performance  in  reverberation  limited  conditions 
(Fig.  03-A.1 ).  In  hard  bottoms,  in  which  sediment  sound  velocity  is  greater  than  in  the  water 
column,  penetration  of  sound  into  the  sediment  is  limited  to  high  grazing  angles  by  the 
critical  angle  of  total  reflection.  The  operating  range  is  therefore  limited  to  approximately 
twice  the  height  of  the  sonar  above  the  seafloor,  a  significant  limitation,  especially  in  the 
very  shallow  waters  where  mine  burial  is  most  probable.  The  concept  does  become  feasible 
however,  if  the  sonar  is  deployed  by  means  of  autonomous  underwater  vehicles. 

GOATS 

To  experimentally  quantify  the  detection  performance  of  such  a  sonar  system,  a  parametric 
source  was  deployed,  generating  sound  in  the  2-15  kHz  band  and  a  receiver  consisting  of 
a  segment  of  towed  array  12  m  in  length,  deployed  statically  7.5  m  above  the  seafloor 
(Fig.  03-A.2).The  source  was  mounted  on  a  24  m  underwater  rail  and  placed  in  the  middle 
of  the  water  column,  in  20  m  of  water.  Thin-walled  spherical  shells  1  m  in  diameter  SI ,  S2 
and  S3,  were  buried  at  depths  of  0  m,  0.5  m  and  1  m  respectively  from  the  centre  of  the 
spheres. 
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Figure  03-A.2  Experimental  configuration. 


The  data  were  processed  in  various  frequency  bands  and  the  signal-to-reverberation  ratio 
for  S2  computed  as  a  function  of  the  incident  grazing  angle  (Fig.  03-A.3).There  is  significant 
signal  excess  at  the  higher  incident  grazing  angles  and  detection  becomes  more  difficult  as 
the  angle  is  reduced. 

It  is  apparent  that  detection  is  better  in  the  2-8  kHz  band  than  in  the  2-15  kHz  band,  which 
highlights  another  characteristic  of  sonar  detection  of  buried  objects,  namely  that  increasing 
bandwidth  does  not  necessarily  improve  detection  performance  in  reverberation.  This  is 
because  the  reduction  of  the  seafloor  reverberation  level  is  offset  by  the  increased  attenuation 
of  the  subbottom  echo  at  the  higher  end  of  the  frequency  band.  The  opposite  situation  is 
obtained  in  the  2-5  kHz  band,  where  the  increase  in  the  level  of  the  subbottom  echo  fails  to 
compensate  for  that  of  seafloor  reverberation,  due  to  the  loss  in  bandwidth  and  angular 
resolution  of  the  receiver  array. 

As  a  preliminary  assessment  of  the  feasibility  of  buried  mine  classification,  high  resolution 
images  of  a  buried  1  m  air-filled  spherical  shell  and  a  buried  water-filled  cylinder,  2  m  in 
length  and  0.5  m  in  diameter,  were  formed  with  the  12  m  towed  array  (Fig.  03-A.4).  The 
array  was  focused  at  short  range  (~7.5  m)  to  produce  images  which  reveal  interesting 
target  features.  Sound  is  scattered  specularly  by  the  front  face  of  the  sphere  and  by  the 
front  and  back  faces  of  the  cylinder.  In  addition,  in  the  case  of  the  spherical  shell,  a  second 
highlight  is  visible  which  corresponds  to  a  resonant  effect  created  by  waves  travelling  within 
the  metallic  shell. 

Future  work  will  focus  on  reproducing  similar  detection  performance,  with  more  compact 
sonar  designs  suited  to  AUV  operations  and  on  classification. 


Figure  03-A.3  Signal  to 
reverberation  ratio  (dB)  as  a  function 
of  grazing  angle  and  frequency  fora 
buried  sphere  (TS=- 16  dB). 
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Figure  03-A.4  High  resolution  images  of  buried  sphere  and  cylinder. 
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Results  from  (terminated  Project  031-2)  on  temporal 
classification  contributing  to  Projects  03-A  and  03-B 

In  order  to  validate  mine  classification  techniques  exploiting  the  aspect  dependence  of  the 
acoustic  signature,  theTASCOE  (TArget  Scattering  in  controlled  Environment)  experiment 
was  conducted  at  Marciana  Marina  (Elba).  The  acoustic  devices  are  installed  on  a  tower 
mounted  on  a  24  m  linear  rail  (Fig.  03-A. 5).  At  the  top  of  the  tower,  a  frame  which  may  be 
horizontally  and  vertically  positioned  supports  a  high  frequency  (325  kHz)  interferometric 
sonar  and  a  low  frequency  sonar  with  high  fractional  bandwidth  (2-16  kHz).  The  high 
frequency  sonar  images  the  acoustic  shadow  of  the  objects  at  several  aspects  (coarse 
sampling).  The  low  frequency  sonar  acquires  the  aspect  dependence  of  echoes  with  an 
extremely  fine  sampling  in  azimuth.  The  interferometric  capability  of  the  high  frequency 
sonar  is  used  to  compute  the  fine  scale  bathymetry  of  the  experimental  site. 


Figure  03-A.6  shows  the  improvement  of  shadow  based  mine  classification  when  combining 
three  views  at  intervals  of  45  °.  An  exercise  mine  and  a  rock  were  acquired  at  low  frequency 
with  an  aspect  resolution  of  1  °.The  difference  between  the  two  targets,  clearly  seen  in  the 
frequency  domain  using  spectrum  variation  when  circumnavigating  the  object,  is  also  visible 
in  the  spatial  domain  when  tomographic  reconstruction  is  applied  (Fig.  03-A.7).  Another 
way  to  represent  the  aspect  variations  of  the  target  signature  is  to  reconstruct  a  temporal 
signal  from  the  aspect  varying  spectrum.  When  the  reconstruction  is  carried  out  in  the 
audio  frequency  band,  a  rock  and  a  geometrical  object  produce  different  sounds  (see 
Fig.  03-A. 8). 

When  the  object  is  classified  as  man-made  and  considered  as  a  potential  mine,  additional 
geometrical  and  elastic  properties  can  be  estimated  from  the  multi-aspect  resonance 
scattering  analysis,  as  demonstrated  in  previous  work  on  water-filled  steel  cylinders.  The 
same  methodology  is  applied  to  an  exercise  mine  with  a  fibre-glass  shell  which  is  shown  to 
generate  both  diffraction  (Fig.  03-A.9a)  and  elastic  features  (Fig.  03-A.9b).The  feasibility 
of  classification  schemes  based  on  resonance  scattering  analysis,  developed  and  tested 
on  free-field  and  proud  targets  was  demonstrated  on  buried  targets,  from  analysis  of 
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Figure  03-A.8  Reconstruction 
of  time  series  in  the  audio  band 
from  the  aspect  varying 
spectrum  of  a  rock  (red)  and  a 
cylinder  (blue).  The  input  data 
are  generated  by  2D  time 
domain  finite  difference  (TDFD) 
modelling  tools. 


Figure  03-A.9  Rigid  and  elastic  feature 
identification  of  a  cylindrical  exercise  mine  proud 
on  a  sandy  bottom,  (a)  Diffraction  feature  analysis: 
axial  (dashed)  and  radial  (solid)  mode  loci,  (b) 
Elastic  feature  preliminary  analysis:  the  curves 
corresponding  to  shear  (dashed)  and  S0  Lamb- 
type  (solid)  wave  modes  are  obtained  by 
connecting  the  detected  energy  peaks  (shown  by 
circles). 


b 


(a) 


(b) 


Figure.  03-A.10  GOATS  98  data-model  comparison  of  time  and  frequency  response  of  1  m0  air-filled  sphere  in 

free  space  (a)  and  flush  buried  (b).  The  model  assumes  the  sphere  in  the  free  field  even  when  buried.  Main  results  of 
elastic  wave  analysis  are  included. 
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Project  03-B  Remote  large  area  search  in  coastal  waters 
(Alliance  days  -13,  Manning  days  -  3) 


Operational  relevance 

Autonomous  underwater  vehicles  (AUV)  equipped  with  high  resolution  sonar  and  precision 
navigation  sensors  and  systems  will  contribute  significantly  to  the  reduction  of  NATO 
MCM  capability  shortfalls  identified  in  MO20 1 5 


In  order  to  improve  navigational  accuracy,  critical  to  successful  AUV  operations,  data  from 
traditional  navigation  sensors  is  combined  with  that  of  high  resolution  mine  hunting  sonar 
on  the  AUV.  The  traditional  sensor  package  consisting  of  an  inertial  navigation  system 
(INS),  is  converted  into  Aided  INS  (AINS)  by  the  addition  of  Doppler  velocity  log  and 
intermittent  DGPS  fixes,  the  data  from  which  are  fused  by  an  extended  Kalman  filter.  Although 
the  recently  acquired  strapdown  INS  is  one  of  the  best  available,  the  accuracy  of  positioning 
remains  insufficient,  in  particular  for  synthetic  aperture  sonar,  (SAS). 

AINS  performance  analysis,  using  a  navigation  simulator/estimator  provided  by  the 
Norwegian  Defence  Research  Establishment  (NDRE),  demonstrated  that  the  attitude  error, 
shown  in  Fig.03-B.1  (a)  is  negligible,  due  to  very  low  sensor  noise  of  the  selected  INS.  The 
position  error  shown  in  Fig.  03-B.1  (b)  is  dominated  by  quadratic  drift  due  to  accelerometer 
bias.  The  SAS  motion  compensation  technique  is  expected  to  significantly  reduce  this  drift. 

The  proposed  approach  of  data-driven  motion  compensation  is  the  Displaced  Phase  Centre 
Antenna  (DPCA)  technique,  which  has  been  demonstrated  to  be  precise  to  a  fraction  of  the 
acoustic  wavelength  (1  cm  for  a  150  kHz  sonar).  The  operationally  significant  feature  of 
DPCA  is  that  it  is  unaffected  by  seafloor  topography. 


The  100  kHz  programmable  sonar  was  integrated  into  a  passive  towfish  which  was 
successfully  deployed  from  R/V  Alliance  in  June  during  the  MASAI’99  (Multi-Aspect  Synthetic 
Aperture  Imaging)  experiment.  This  system  (Fig.  03-B.2)  is  being  used  to  experimentally 


Standard  deviation  of  position  displacement  error  (x  component) 
over  60  runs. 


Figure  03-B.1  Analysis  of  simulated  AINS  short-term  residual  error  of  attitude  (a)  and  position  displacement  Figure 
03-B.  1  (b)  shows  only  the  x  components,  but  analogous  results  have  been  obtained  for  the  other  two  components. 
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Figure  03-B.2  The  100  kHz 

programmable  sonar  deployed  from 
R/V  Alliance  in  June  99  during  the 
MASAI’99. 


validate  a  multi-aspect  SAS  concept  (Fig.  03-B.3).  In  addition  to  the  extremely  high  cross¬ 
range  resolution  provided  by  SAS,  which  allows  shadow  classification  to  be  performed  at 
increased  ranges,  better  target  recognition  is  expected  to  result  from  the  acquisition  of 
multiple  high  resolution  looks.  The  benefits  of  multi-aspect  processing  are  illustrated  by 
images  of  a  vessel  sunk  during  WWII,  obtained  with  the  PS  during  MASAI’99  (Fig.  03-B.4). 
SAS  processing  is  continuing.  The  challenge  is  that  the  motion  compensation  has  to  be 
performed  by  DPCA  alone,  due  to  the  lack  of  INS  in  the  MASAI’99  experiment.  Preliminary 
DPCA  estimates  are  shown  in  Fig.  03-B.5.  The  standard  deviations  are  close  to  theoretical 
expectations  but  the  attitude  rates  are  offset  for  reasons  that  are  as  yet  unexplained. 


The  AUV  specified  and  procured  in  support  of  REA  and  MCM  studies,  is  a  modified  Ocean 
Explorer  (OEX),  of  modular  structure,  manufactured  by  Florida  Atlantic  University.  The  rear 
propulsion  module  houses  the  components  common  to  all  operations,  including  power, 
propulsion,  navigation  and  control.  One  or  more  forward  payload  sections  are  dedicated  to 
mission  specific  sensors  A  common  nose  cone  section  includes  a  video  camera  and  a 
recovery  line  system. The  MCM  payload  (Fig  03-B-6)  will  contain  the  recently  procured  INS 
and  a  wideband  interferometric  array  for  synthetic  aperture  sonar  studies.  The  array  will  be 
specified  and  procured  in  2000  and  installed  in  the  OEX  in  2001 . 


Squinted 
backwards  SAS 


Broadside 

SAS 


Squinted 
forward  SAS 


Figure  03-B.3  Geometry  of  the  multi-aspect 
synthetic  aperture  sonar. 


—  36 


ping  to  ping  surge(m) 


Figure  03-  B.  4  Images  of  a  vessel  sunk  during  WWII  obtained  at  multiple  aspects  with  the  100  kHz  programmable 
sonar  during  MASAI’99. 


Figure  03-B.5  DPCA  micronavigation  of  MASAI’99  (carrier  frequency  100  kHz,  20 
kHz  bandwidth,  128  channels  spaced  at  1/2)  for  three  different  ranges  giving  three 
independent  estimates  of  the  motion  in  the  slant-range  plane.  The  yaw  estimate  is 
biased,  possibly  because  of  calibration  errors  of  the  physical  array. 


The  1 00  kHz  programmable  sonar  deployed 
from  R/V  Alliance  in  June  99  during  the 
MASAT99. 
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Figure  03- B. 6  The  MCM  autonomous  underwater  vehicle  consists  of  a  propulsion  section  which 
houses  the  common  navigation  components,  a  payload  sections  with  the  INS  and  the  interferometric 
wide  band  sonar  and  a  nose  cone  section  with  a  video  camera  and  a  recovery  line  system. 
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Ocean  Engineering  Group  since  1982,  working  mainly  on  buoy  technology 
and  towed  body  design. 


After  graduating  from  the  University  of  Bath  with  a  BSc  in  physics,  Luigi 
Troiano  was  a  design  engineer  for  a  multibeam  imaging  sonar  system,  at 
Ulvertech  Ltd.  Since  joining  SACLANTCEN  in  1987  he  has  provided 
engineering  support  to  scientific  sea-trials  in  the  fields  of  low  frequency 
active  sonar,  mine  countermeasures  and  environmental  acoustics. 
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Project  03-C:  Mine-ship  interaction 


Operational  relevance 

During  a  short  term  operation,  reduce  the  risk  to  a  target  transiting  a  minefield,  contain¬ 
ing  unknown  mine  mechanisms,  to  less  than  5%. 


Mine  ship  interaction 

Research  focused  on  mine  jamming  with  a  subsidiary  task  to  develop  planning  algorithms 
for  use  in  Target  Sweeping  Mode  operation. 

Mine  sweeping 

Many  of  the  mines  likely  to  pose  a  threat  to  NATO  forces  are  influence  mines,  equipped 
with  sensors  to  detect  the  presence  of  a  suitable  target  vessel  from  the  various  influence 
fields  (such  as  magnetic  and  acoustic)  radiated  by  ships. 

NATO  nations  have  the  ability  to  sweep  these  mines  using  equipment  to  generate  the 
required  influence  fields.  This  equipment  is  towed  through  the  water  by  surface  or  airborne 
assets  to  actuate  mines  by  fooling  their  sensors  and  logic.  Traditionally,  influence  sweeping 
has  exploited  particular  weaknesses  in  mine  design  and  has  assumed  a  known  mine  threat. 
This  is  commonly  termed  minesweeping  mode  (MSM). 


Mine  jamming 

Minefields  are  rendered  safe  by  time  consuming  minesweeping  or  mine  hunting  operations. 
There  are  occasions  when  vessels  are  required  to  traverse  a  mined  area  before  mine 
clearance  operations  can  be  completed.  Project  03-C  is  studying  techniques  which  will  jam 
the  mine’s  decision  logic  or  algorithm  allowing  the  target  to  proceed  safely. 
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When  time  is  unconstrained  and  the  level  of  conflict  is  low,  time  can  be  spent  on  surveys 
and  surveillance.  As  time  becomes  more  constrained  and  as  the  conflict  intensity  increases, 
effort  will  be  spent  on  rapid  area  assessment  and  covert  surveillance.  Subsequently  the 
MCM  operation  includes  remote  mine  hunting  and  minesweeping,  followed  by  mine 
avoidance  in  support  of  an  in-stride  operation.  Finally,  when  conflict  intensity  is  high  and 
time  very  limited,  mine  jamming  becomes  important. 

Several  mine  jamming  concepts  were  proposed  during  the  pre  study  phase  of  this  work1 
and  computer  models  written  to  simulate  these  concepts. 


Figure  03-C.1  illustrates  the  performance  of  one  mine  against  a  frigate  with  an  outline  of 
the  ship  and  a  damage  circle.  The  mine  actuations  are  represented  by  the  grey  symbols.  It 
is  clear  that  the  mine  presents  a  significant  risk  to  the  target.  When  a  mine  jammer  protected 
this  same  target,  the  mine  actuated  in  less  than  0.05%  of  the  engagements  and  no  mine 
actuations  presented  a  threat  to  the  target.  The  trial  data  (Fig.  03-C.2)  confirmed  these 
findings. 


Fig  u  re  03-C.  2  Trial  data  showing  the  target  ship  protected  by 
a  jamming  signal. 


The  Joint  Research  Project  partners  met 
in  September  to  discuss  progress  and 
the  way  ahead.  A  second  trial  was 
suggested  for  late  2001  at  the  Coastal 
Systems  Station  in  Panama  City. 

Most  of  the  mine  jamming  techniques 
considered  are  novel  applications  of 
existing  equipment.  Initial  indications 
from  the  trial  show  that  these  techniques 
are  promising  and  if  successful,  will 
augment  the  MCM  armoury  without 
significant  additional  cost  to  the  NATO 
Nations. 


Target  sweeping  mode 
planning 

NATO  nations  have  developed  an 
influence  minesweeping  capability 
against  an  unknown  mine  threat  by 
configuring  the  sweeping  equipment  to 
emulate  the  characteristics  of  a  ship’s 
influence  signatures.  This  concept  is 
termed  target  sweeping  mode  (TSM). 


The  objective  of  mine  countermeasures  (MCM)  operations  planning  is  to  minimize  the  risk 
to  subsequent  target  shipping,  the  time  taken  and  the  risk  to  the  MCM  forces.  The  current 
NATO  doctrine  for  the  planning  and  evaluation  of  influence  minesweeping  operations  is 
based  on  a  priori  knowledge  of  performance,  described  in  terms  of  sweep  width  and  actuation 
probability  of  the  minesweeping  system  against  the  threat  mine  type(s).  As  it  is  not  possible 
to  quantify  sweep  performance  against  an  unknown  mine  threat,  this  approach  may  no 
longer  represent  an  optimum  tactic  for  the  employment  of  a  TSM  minesweeping  system. 


1  SACLANTCEN  SR-271 
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The  aim  of  this  study  is  to  investigate  whether  the  assumptions  underlying  the  current 
NATO  algorithms  remain  valid  for  TSM  and  to  develop  methodologies  for  the  optimum 
employment  of  TSM  minesweeping  systems  and  the  evaluation  of  risk  reduction  to  the 
target  vessel. 


During  1999,  a  technical  solution 
was  developed,  which  is  an 
extension  of  the  current  NATO 
MCM  planning  and  evaluation 
procedures  implemented  within  the 
MCM  EXPERT  tool.  The  planning 
algorithm  uses  either  exhaustive 
search  or  simulated  annealing  to 
optimize  the  minesweeping  effort. 
The  selection  of  technique  depends 
on  the  number  of  possible 
combinations  in  the  scenario  and 
the  speed  of  the  PC. The  evaluation 
algorithm  determines  the  remaining 
risk  as  a  function  of  the  aggregate 
actuation  width  of  the  mine  against 
the  target  vessel.  An  example  of  the 
evaluation  algorithm  output  is 
shown  in  Figure  03-C.3. 


Figure  03-C.3  TSM  evaluation:  the  risk  to  the  target  plotted  against 
the  aggregate  actuation  width  of  the  mine  against  the  target. 


The  algorithms  have  been  implemented  in  EXCEL  97  VBA  for  development  and  testing.  At 
a  meeting  with  the  project  sponsor,  Belgium,  in  November,  the  concepts  were  discussed 
and  the  algorithms  and  prototype  user  interface  demonstrated  for  a  standalone  prototype 
TSM  planning  and  evaluation  software  tool  (Figure  03-C.4).  It  was  agreed  that  the 
development  of  the  software  (Planning  and  Evaluation  for  TSM  -  PET)  should  be  continued. 


This  task  is  scheduled  to  be  completed  at  the  end  of  May  2000  with  delivery  to  the  sponsor 
of  the  PET  prototype  software,  source  code,  user  guide,  help  files  and  a  SACLANTCEN 
report  describing  the  algorithms  implemented. 


Figure  03-C.4  Example  of  prototype  TSM  planning  and 
evaluation  software  user  interface. 


3i  coastal  bystems  Station,  Panama  City,  FL,  October  1999. 


Glynn  Field  joined  the  Admiralty  Underwater  Weapons  Establishment ,  AUWE 
in  1971.  During  his  early  years  he  wrote  finite  element  software  to  predict  the 
dynamic  behaviour  of  towed  arrays.  He  later  developed  computer  models  to 
predict  the  hydrodynamic  properties  of  flexible  sonar  domes.  In  1981,  he 
joined  the  Mine  Warfare  Department  where  he  was  responsible  for  research 
into  Mine  sensors.  In  1985 he  became  responsible  for  work  on  Mine  Algorithms 
and  Simulation .  One  of  his  many  responsibilities  was  the  Total  Mine  Simulation 
System,  TMSS.  TMSS  has  become  the  focus  of  international  collaboration. 
Versions  of  TMSS  are  installed  in  Australia,  Canada,  France,  Germany, 
Netherlands,  New  Zealand,  Norway,  Sweden,  US  and  SACLANTCEN,  where 
Glynn  Field  has  been  Project  Leader  of  Mine/Ship  Interaction  since  1994 
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Project  03-D:  Environmental  impact  of  MCM  on  sonar 
design  and  performance 
(Manning  days  -  36) 


Operational  relevance 

An  increased  understanding  of  high  frequency  acoustic  interaction  with  the  seabed  in 
shallow  water  to  support  MCM  modelling  requirements. 

In  shallow  water  at  high  frequencies,  the  rapid  spatial  changes  in  water  column  and 
seabed  properties  are  of  significance  in  the  understanding  and  prediction  of  high- 
resolution  sonar  response ,  particularly  synthetic  aperture  sonar  (SAS).  Of  particular 
concern  are  the  levels  and  spatial  coherence  of  high  frequency  acoustic  signals  scat- 
tered  from  the  seabed  in  high-energy  environments.  Understanding  of  near  critical 
angle  penetration  of  the  seabed  and  scatter  from  buried  targets  is  essential  to  optimal 
design  of  sonars  for  countering  a  buried  mine  threat.  The  need  to  map  seabed  areas 
where  mines  are  likely  to  become  buried  is  necessary  for  planning  purposes  and  is 
being  addressed  through  the  inversion  of  acoustic  data  particularly  that  from  wide 
swath  sonars . 


Acoustic  fluctuations  due  to  water  column  effect 

The  effect  of  sea-surface  waves  on  the  travel  time  of  acoustic  signals  in  shallow  water  has 
been  investigated  theoretically  and  experimentally  in  the  context  of  synthetic  aperture  sonar. 
Vertical  and  horizontal  displacement  of  water,  a  function  of  depth  and  azimuth  results  from 
sea-surface  waves.  The  horizontal  displacement  velocity  produces  a  vector  addition  to  the 
mean  sound  speed.  Vertical  displacement,  in  the  presence  of  vertical  gradients  in 
temperature  and  salinity  produces  a  scalar  addition  to  the  mean  sound  velocity.  For  one 
way  transmission,  scalar  and  vector  effects  will  influence  the  signal.  A  backscattered  signal 
will  be  affected  only  by  the  scalar  effect  if  it  can  be  assumed  that  the  time  scale  associated 
with  the  variability  is  much  longer  than  the  two  way  travel  time.  Such  effects  have  been 
modelled  and  verified  against  experiment.  When  this  knowledge  is  applied  to  SAS  in  which 
displaced  phase  centre  antenna  (DPCA)  techniques  are  used  to  micronavigate  the  synthetic 
aperture,  the  study  has  confirmed  that  fluctuations  of  the  medium  can  be  treated  as  an 
additional  motion  of  the  AUV  and  thus  corrected  for,  provided  the  ping  repetition  period 
remains  small  compared  to  the  temporal  coherence  time. 


Physics  of  seabed  acoustic  scatter  and  penetration 

Repeated  experiments  have  demonstrated  “anomalous”  penetration  of  acoustic  energy 
into  the  seabed  below  the  critical  angle.  A  major  unanswered  question  is  whether  this 
energy  can  be  exploited  by  MCM  systems  for  detecting  buried  mines.  An  experiment 
(APEx99)  was  designed  to  address  questions  that  recent  SACLANTCEN  investigations 
have  failed  to  resolve: 

•  the  arrival  angle  of  acoustic  energy 

•  the  mechanism  for  subcritical  penetration  at  higher  frequencies 

•  the  modelling  technique  which  is  most  applicable  to  the  interface  scattering  component. 
The  main  experimental  goal  of  APEx99  was  to  measure  acoustic  penetration  into  a  sandy 
seafloor,  above  and  below  the  critical  grazing  angle,  over  a  wide  frequency  range  and  the 
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Accelerometers  Hydrophone  Cardioid 
Pairs  Arrays 


Figure  03-D.1  Sensor  hydrophone  array. 


Figure  03-D.2  Particle  motion  measured  by  the 
buried  accelerometers  for  a  2  m,  11  kHz ,  CW pulse. 
The  beam,  which  is  incident  on  the  seafloor  at  37° , 
is  largely  just  above  critical  angle  and  is  refracted, 
after  penetrating  the  seafloor,  to  22  °.  The  output  of 
each  accelerometer ;  displayed  in  the  top  two  graphs 
are  combined  to  give  particle  motion  and  thus  the 
signal  arrival  direction.  The  colour  coding  allows 
spurious  effects  at  the  start  and  end  of  each  pulse 
to  be  identified. 


Arrival  Arx^e 
30° 


angle  of  arrival.  The  latter  is  an  important  clue  to  the  penetration  mechanism  (e.g.,  Biot 
slow  wave  or  interface  roughness  scattering).  Acoustic  measurements  in  conjunction  with 
ground  truth  measurements  will  allow  discrimination  between  candidate  mechanisms  for 
penetration  and  a  test  of  the  robustness  of  the  penetration  prediction  models,  which  have 
been  developed. 

The  Apex99  experiment  used  a  diver-buried  vertical  array  of  accelerometers  in  a  sand 
seabed  in  conjunction  with  two  small  hydrophone  arrays  of  3  sensors  each,  closely  positioned 
at  the  same  depth  (Fig.  03-D.1).  Short  acoustic  pulses  were  transmitted  from  conventional 
sources  mounted  on  a  tower  equipped  with  a  pan  and  tilt  mechanism.  Different  acoustic 
energy  incidence  angles  on  the  seafloor  were  achieved  by  moving  the  tower,  tilting  the 
transmitter  head  and  by  raising  the  height  of  the  telescopic  tower.  This  configuration  allowed 
the  directionality  of  the  incoming  acoustic  waves  and  transmitted  pressure  levels  to  be 
determined  as  a  function  of  frequency  (3-30  kHz)  and  angle  (5°-35°).The  buried  directional 
sensors  allowed  the  propagation  direction  of  the  penetrating  acoustic  waves  to  be  accurately 
identified  (Fig.  03-D.2).  This  technique  is  a  considerable  improvement  on  measurements 
made  using  pressure  sensors,  as  the  acoustic  arrival  angles  are  difficult  to  obtain  with  a 
sparse  array.  The  use  of  linear  as  opposed  to  parametric  sources  will  simplify  the  analysis 
of  the  penetration  data.  Accurate  measurements  of  seafloor  roughness  were  made  daily 
with  the  new  digital  photogrammetry  system  and  cores  were  analyzed  for  their  physical 
properties. 

High  frequency  volume  scattering  experiment  (31  May  -18  June 

1999) 

This  trial  was  the  second  milestone  of  the  Defence  Evaluation  and  Research  Agency  (DERA) 
Joint  Research  Project  Wideband  acoustic  interaction  with  the  seabed  (April  1998  -  April 

2000) .  The  emphases  were  on  measuring  high-frequency  volume  scattering  from  different 
seafloors  (from  silt  to  coarse  sand)  and  extraction  of  statistical  properties  of  backscatter  as 
a  function  of  ensonified  area,  grazing  angle  and  frequency  in  relation  to  seabed  properties. 
The  DERA  wideband  research  sonar  (1 1 0-190  kHz)  and  the  SACLANTCEN  Simrad  EM3000 
(300  kHz)  bathymetry/backscatter  measurement  system  were  deployed.  Progress  in  this 
area  will  be  of  direct  benefit  to  the  interpretation  and  use  of  high-resolution  sonar  imagery. 


Figure  03-D.3  Scattering 
strength  versus  range  and 
azimuth  at  140  kHz  for  a 
hard,  rough  bottom  (left)  and 
a  soft,  inhomogeneous 
bottom  (right). 
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An  example  of  acoustic  scattering  data  obtained  with  the  wideband  system  is  shown  in 
Fig.  03-D.3.The  stronger  returns  from  the  soft  bottom  provide  evidence  of  the  importance 
of  volume  scattering  even  at  frequencies  above  1 00  kHz.  Acoustic  data  was  complemented 
by  the  EM3000  raw  data  at  300kHz  over  many  beams  with  swath  bathymetry  and  the 
sidescan  sonar  operating  at  the  same  frequency  as  the  EM3000.  Supporting  environmental 
data  including  seabed  roughness  spectra  derived  from  stereo  photography,  grab  samples, 
sediment  cores,  conductivity,  temperature  depth  (CTD)  and  video.  The  seafloor  ground 
truth  cores  were  analyzed  using  an  x-ray  computed  tomography  (CT)  scan  technique  which 
provided  important  information  on  the  internal  three-dimensional  structure  of  the  sediment 
volume,  such  as  density  variations  and  the  distribution  of  shell  detritus  and  worm  tubes 
(Fig.  03-D.4). These  features  have  a  significant  impact  on  high-frequency  volume  scattering 
(Fig.  03-D.3)  and  their  quantification  is  an  important  step  in  being  able  to  model  seafloor 
acoustic  response  which  causes  strong,  high  frequency  volume  scattering. 
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Fusion  of  acoustic  data  for  sea-bed  mapping 

A  procedure  based  on  simultaneous  operation  of  four  acoustic  sensors  is  under  development 
to  provide  in  a  rapid,  quasi-automatic  and  unsupervised  manner,  complementary  information 
about  an  area  of  interest.  A  300  kHz  multibeam  echosounder  provides  high  resolution 
bathymetry  and  backscattering  strength  mapping.  A  sidescan  sonar  operating  at  the  same 
frequency  provides  information  about  low  grazing  backscattering  and  textural  imaging. 
Signals  from  a  vertical  echosounder  are  processed  to  discriminate  between  seabed  types 
and  complement  information  on  bathymetry.  A  subbottom  profiler  provides  vertical  information 
on  seabed  substructure.  The  combination  of  these  four  sensors  with  the  use  of  unsupervised 
segmentation  algorithms  produces  high  resolution  mapping.  Combining  bathymetric,  textural 
and  discriminating  information  improves  analysis  and  reduces  ambiguity  in  the  classification 
process.  An  interesting  aspect  of  this  procedure  is  its  use  of  commercially  available 
instruments.  The  procedure  has  been  successfully  tested  at  different,  well  known  sites 
using  video  imaging,  sediment  grab  samples  and  photography.  Figure  03-D.5  shows  an 
example  of  complementary  data,  unsupervised  segmentation  of  100kHz  sidescan  sonar 
with  supervised  segmentation  based  on  angular  dependence  of  backscatter  coefficients 
from  swath  bathymetric  sonar  at  300kHz.  The  further  development  of  segmentation  and 
classification  using  fusion  from  several  sensors  will  be  conducted  with  the  objective  of 
producing  the  MCM  ATP24  classification  using  an  autonomous  underwater  vehicle  as  the 
sensor  platform. 


Figure  03-D.5  Segmentation  of  the  seabed 
off  the  northern  coast  of  Corfu.  Upper  is 
unsupervised  segmentation  of  100kHz 
sidescan.  The  lower  is  supervised 
segmentation  using  the  angular 
dependence  of  the  backscatter  coefficient 
derived  from  the  combined  bathymetry- 
backscatter  mode  of  the  300kHz  Simrad 
EM3000 
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classification  and  geographic  information  systems. 
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Project  03-E:  High  fidelity,  broadband  structural 
acoustic  scattering  models  for  buried  and  partially 
buried  mines 


Operational  relevance 

Accurate  computer  models  of  scattering  of  sonar  signals  from  seabed  mines  are  essential 
for  significantly  improving  the  effectiveness  of  acoustic  minehunting  operations.  The  models 
improve  the  predictive  capabilities  for  detection  and  classification  systems,  thus  providing 
an  essential  component  of  predictive  sonar  performance  models  for  MCM.  Computer  models 
also  increase  understanding  of  how  various  physical  features  in  the  mines  and  seabed 
affect  the  scattering  of  sonar  signals,  thereby  enabling  the  design  of  acoustically  stealthy 
mines  to  counter  an  adversary's  MCM. 


Phase  one  of  this  project  (1995-1998)  focused  on  two-dimensional  scattering  from  simple, 
highly  idealized  mine  shapes  (spheres  and  cylinders). 

The  objective  of  the  current  phase  is  to  employ  state-of-the-art  finite  element  (FE)  techniques 
to  develop  a  more  realistic  and  useful  computer  modelling  capability  for  use  in  all  research 
activities  related  to  mine  hunting  sonar  concept  developments  and  performance  predictions. 

The  following  system  requirements  for  the  computer  model  have  been  established: 

•  High  fidelity:  Highest  priority  is  the  ability  to  model  real  (3-D)  mines  to  sufficient 
structural  detail  (e.g.,  internal  and  external  structural  features,  electronics,  explosives, 
ballast,  etc.)  to  provide  high  accuracy  relative  to  real  world  conditions,  thereby  providing 
good  agreement  with  experimental  data  measured  at  sea. 

•  Broadband:  As  different  types  of  information  are  revealed  across  the  frequency 
spectrum,  it  is  necessary  to  be  able  to  model  from  low  to  very  high  frequencies. 

•  Computationally  efficient:  The  preceding  requirements  impose  severe  demands 
on  the  computational  efficiency  of  the  mathematics,  numerical  algorithms  and 
computer  hardware.  Therefore  state-of-the-art  techniques  are  essential  to  providing 
a  modelling  tool  which  is  affordable  on  in-house  computers. 

•  User-friendly:  The  code  is  intended  to  be  used  by  scientists  who  are  not  specialists 
in  FE  analysis,  e.g.,  acousticians,  signal  processors,  oceanographers.  Therefore  the 
code  should  be  highly  automated,  in  addition  to  having  a  good  graphical  user  interface. 

To  meet  all  four  requirements,  two  types  of  technology  are  being  employed: 

•  a  patented1  mathematical  acoustic  technique  using  computationally  efficient  “infinite 
elements,” 

•  a  unique,  commercial  FE  technology  which  is  considerably  more  advanced  than  that 
in  conventional  FE  codes. 


1  US  Patents  5604891 , 5604893, 5963459 
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Mathematical  acoustic  technique 

The  acoustic  environment  surrounding  a  mine  is  the  ocean  and  seabed,  both  of  which  are 
extremely  large  spaces  compared  to  the  dimensions  of  a  mine,  too  large  to  be  modelled  by 
conventional  FE  methods.  Many  mathematical  techniques  for  modelling  such  large  spaces 
have  evolved,  but,  prior  to  the  early  ‘90s,  none  had  been  numerically  efficient.The  patented 
“infinite  element”  technique  is  more  than  400  times  faster  than  competing  techniques  in 
achieving  the  same  solution,  to  the  same  degree  of  accuracy  in  benchmark  tests  and  has 
been  described  as  “one  of  those  rare  breakthroughs  that  caused  the  entire  field  of 
computational  acoustics  to  re-evaluate  itself  and  proceed  in  new  directions”2. 

Figure  03-E.1  illustrates  how  finite  and  infinite  elements  will  be  applied,  using  a  partially 
buried  Manta  mine.  Finite  elements  will  model  the  mine  (e.g.,  casing,  internal  structure, 
mechanisms,  explosive  and  ballast)  and  a  small  volume  of  adjacent  water  and  seabed  out 
to  an  ellipsoidal  surface,  which  closely  circumscribes  the  mine.  For  the  remaining  space 
“out  to  infinity”,  a  single  layer  of  the  new  infinite  elements  is  sufficient  to  describe  the  full  3- 
D  physics  of  any  scattered  field  to  any  desired  accuracy. 


Figure  03-E.1  Computer 
modelling  approach  for 
solving  mine  scattering 
problems  by  finite  and  infinite 
elements. 


Finite  element  technology 

Although  the  infinite  element  technology  is  licensed  to  two  commercial  acoustic  code 
developers  and  available  in  their  current  code  releases,  those  codes  lack  other  capabilities 
for  meeting  the  above  system  requirements.Therefore,  the  Centre  has  acquired  commercial 
FE  software  packages,  which  employ  a  unique  state-of-the-art  general  FE  technology, 
which  is  considerably  more  sophisticated  than  in  other,  more  conventional  commercial 
codes.  Following  installation  and  initial  validation  with  vendor-supplied  test  problems, 
validation  is  continuing  with  Centre  test  problems. 

The  software  employs  full  3-D  physics  for  all  modelling,  eschewing  the  traditional  2-D  and 
1-D  engineering  approximations  such  as  plates,  shells  and  beams.  The  additional 
computational  costs,  which  normally  accompany  a  richer  physics  are,  nevertheless, 
dramatically  reduced  through  the  use  of  sophisticated  error  estimation  and  adaptive 
“meshing”  algorithms.  Thus,  the  software  quantitatively  estimates  the  discretization  error 


2  J.T.  Oden,  Director,  Texas  Institute  for  Computational  and  Applied  Mathematics. 
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throughout  a  model,  i.e.,  the  departure  from  the  exact  solution,  which  is  used  to  automatically 
modify  the  “mesh”  (computational  grid)  until  a  near-optimal  mesh  is  created,  which  achieves 
a  user-specified  error  at  the  lowest  cost  in  computer  resources.  The  result  is  a  model  with 
near-maximum  computational  efficiency.This,  in  turn,  permits  modelling  to  higher  frequencies 
as  computational  costs  increase  exponentially  with  frequency. 

Operationally,  the  analyst’s  task  is  limited  to  the  creation  of  a  first,  coarse  mesh  of  the  mine, 


casing  of  a  Manta  mine. 


seabed  and  water,  using  a  state-of-the-art  mesh  generator  program.  Then,  upon  specifying 
a  desired  error  (e.g.,  ±  5%  relative  to  the  exact  solution),  the  program  takes  over  and 
follows  the  automated  procedure  until  a  solution  with  that  error  is  found. 

This  process  is  illustrated  by  the  following  FE  analysis  of  a  Manta  mine.  Figure  03-E.2 
shows  a  mesh  for  the  hard-plastic  casing  of  the  mine;  only  180°  of  the  360°  rotationally 


symmetric  structure  is  depicted,  in  order  to  reveal  inner  parts  of  the  casing.  Figure  03-E.3 
shows  the  casing  as  well  as  the  other  internal  structural  components:  the  priming  device 
(green),  TNT  charge  (light  yellow-brown),  booster  (orange)  and  water  (blue),  with  appropriate 
material  properties  assigned  to  each  different  component.  This  model  was  then  subjected 
to  a  purely  structural  load:  a  pressure  exerted  over  a  small  area  on  the  top  and  side  surfaces, 
while  the  mine  sits  on  a  hard,  flat  ground.  The  analyst  requested  that  the  software  produce 
a  solution  in  which  the  maximum  error  (in  strain  energy)  in  any  of  the  finite  elements  be 
less  than  7%.  Figures  03-E.4  and  03-E.5  show  the  resulting  solution,  the  colour  contours 
representing  the  magnitude  of  the  computed  displacement  at  every  point  on  the  surface, 
due  to  the  deformation  of  the  mine  under  the  action  of  the  applied  pressure  load.  The  view 
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Figure  03-E.4  Colour  contours  show 
displacement  magnitude  on  the 
outside  surface  of  the  mine  casing 
(red  is  maximum;  blue  is  minimum). 


in  Fig.  03-E.5  is  of  an  interior  vertical  plane  bisecting  the  mine;  the  plane  contains  the 
central  axis  and  bisects  the  area  of  application  of  the  pressure  load.  A  comparison  of  the 
final  mesh  in  Figs.  03-E.4  and  03-E.5  with  the  initial  mesh  in  Figs.  03-E.2  and  03-E.3 
reveals  some  of  the  mesh  refinement  that  was  performed  automatically  by  the  software  to 
achieve  the  requested  error. 

Figures  03-E.2  to  03-E.5  illustrate  a  so-called  static  analysis,  i.e.,  the  equilibrium  response 
of  a  structure  to  a  load  held  constant  in  time. The  software  also  performs  dynamic  analyses, 
i.e.,  the  transient  response  to  a  time-varying  load,  with  solutions  displayed  as  animated 
sequences  of  pictures  like  Figs.  03-E.4  or  03-E.5,  depicting  the  propagation  of  stress  waves 
through  the  mine. 

During  the  year  2000  we  shall  develop  acoustic  finite  elements  (for  fluids),  which  employ 
the  above  error-adaptive  technology.  The  code  will  perform  in  the  same  manner  as  the 
above  structural  code,  except  in  fluids  rather  than  elastic  solids,  i.e.,  generating  animated 
solutions  of  acoustic  wave  propagation  through  fluids. 

This  acoustic  technology  will  then  be  integrated  with  the  structural  technology  to  achieve 
the  above-stated  objective  of  this  project:  a  structural  acoustics  scattering  code,  which  will 
display  animated  sequences  of  pictures  depicting  the  propagation  of  acoustic  waves  as 
they  interact  with  the  mine  and  seabed.  The  code  will  be  directly  applicable  (without 
modification)  to  scattering  from  general  targets,  e.g.,  submarines.  Indeed,  the  acoustic 
infinite  elements  were  originally  developed  for  ASW  applications. 
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Project  03-E  publications  and  presentations 

Burnett,  D.S.  Modelling  scattering  from  objects  near  boundaries  using  infinite  elements.  Acoustical  Society  of  America 
conference,  Berlin,  March,  1999,  invited  paper.  Journal  of  the  Acoustical  Society  of  America,  105, 1999:1 166:3aUWa4. 

Burnett,  D.S.  Modelling  scattering  from  objects  near  boundaries  using  infinite  elements.  4th  International  Conference 
on  Theoretical  and  Computational  Acoustics,  Trieste,  May,  1999,  invited  paper. 

Burnett,  D.S.  Numerical  solution  of  the  exterior  radiation  problem  with  infinite  elements.  Machine  Acoustics  ’99,  Wiesloch, 
Germany,  Sept.  1 999,  invited  plenary  lecture. 

Burnett,  D.S.  Development  of  high  fidelity,  broadband  3-D  structural  acoustic  codes  for  scattering  and  radiation  using 
acoustic  infinite  elements.  Defence  Evaluation  and  Research  Agency,  (DERA),  Rosyth,  Edinburgh,  Oct.  1 999,  invited 
presentation. 

Burnett,  D.S.  Acoustic  infinite  elements,  University  of  Durham,  Mathematics  Department,  October  1 999,  invited  seminar. 


David  Burnett  received  B.S.  and  M.  Eng.  degrees  in  Engineering  Physics 
from  Cornell  University  in  1962,  an  M.S.  degree  in  Engineering  Science 
from  the  California  Institute  of  Technology  in  1963  and  a  Ph.D.  degree  in 
Theoretical  Mechanics  from  the  University  of  California,  Berkeley  in  1969. 
He  worked  more  than  28  years  at  Bell  Laboratories  ( Lucent  Technologies, 
formerly  part  of  AT&T),  primarily  in  undersea  R&D  for  the  U.S.  Navy, 
specializing  in  theoretical  and  computational  mechanics  (elasticity  and 
acoustics).  In  the  ‘80s  and  ‘90s  he  was  a  group  technical  leader  for  the 
development  of  3-D  structural  acoustics  FE  codes  for  AS  W  applications. 
Dr.  Burnett  taught  many  courses  in  finite  element  analysis  in  the  Bell  Labs 
In-Hours  educational  program.  In  1983  he  received  the  title  of  Distinguished 
Member  of  Technical  Staff  and  in  1996  the  title  of  Fellow,  Bell  Laboratories’ 
highest  honour.  Dr.  Burnett  holds  several  patents  in  the  field  of 
computational  acoustics  and  is  the  author  of  Finite  element  analysis:  from 
concepts  to  applications.  He  joined  SACLANTCEN  in  1998. 
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Project  03-F:  Minehunting  sonar  performance  model 


Operational  relevance 

To  provide  NATO  with  an  improved  performance  prediction  tool  for  minehunting,  particu¬ 
larly  in  shallow  water  and  to  provide  the  sonar  performance  parameters  required  as 
inputs  by  the  NATO  MCM  planning  and  evaluation  tool  MCM  EXPERT. 


The  aim  of  Project  03-F  is  to  provide  NATO  with  a  minehunting  sonar  performance  prediction 
tool  which  is  capable  of  interfacing  to  the  NATO  MCM  planning  tool,  MCM  EXPERT  and 
environmental  data  supplied  by  Rapid  Environmental  Assessment  (REA)  and  other  sources. 
The  tool  will  produce  p(y)  curves  (probability  of  detection  as  a  function  of  across  track 
distance)  and  calculate  A  (characteristic  detection  width)  and  Bd  (characteristic  detection 
probability)  values  as  required  by  MCM  EXPERT.  Figure  03-F.1  illustrates  where  the  tool  to 
be  developed  by  Project  03-F  fits  with  MCM  EXPERT  and  REA  data  (or  environmental  data 
from  other  sources). 


Sonar  Performance  Model 


Figure  03-F.1  Interfaces  between  the  minehunting  sonar  performance  model,  environmental  and  target  data  and  the 
NATO  MCM  planning  and  evaluation  tool  MCM  EXPERT. 

This  project  aims  to  incorporate  the  best  acoustic  sub-models  available  from  SACLANT CEN 
and  other  research  centres.  The  research  at  the  Centre  on  target  modelling  (Project  03-E) 
and  seabed  scattering  (Project  03-D)  is  particularly  relevant.  The  Centre’s  acoustic 
propagation  expertise  is  being  used  to  select  a  suitable  propagation  model,  which  includes 
multipaths,  particularly  important  in  shallow  water  environments. 

A  study  is  in  progress  to  define  the  sub-models  which  will  be  used  and  to  develop  a 
preliminary  design.  Visits  to  government  laboratories  in  the  US,  Germany,  UK  and 
Netherlands  have  been  made,  to  establish  which  minehunting  performance  tools  and  sub¬ 
models  exist  and  to  ensure  that  the  tool  is  state-of-the-art.  The  project  will  make  use  of 
existing  tools  and  sub-models  to  avoid  “re-inventing  the  wheel”. 
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The  preliminary  design  work  has  focused  on  adapting  different  acoustic  sub-models  and 
environmental  descriptions,  without  having  to  re-write  code.  This  is  achieved  by  defining  an 
interface”  for  each  major  element  in  the  model.  “Adapter  Classes”  may  be  used  to  convert 
between  incompatible  interfaces.  These  are  analogous  to  electrical  adapters  which  provide 
a  means  of  connecting  incompatible  electrical  plugs  and  sockets.  For  example,  if  the 
environmental  parameters  describing  the  seabed  do  not  match  the  parameters  required  by 
the  bottom  scattering  sub-model,  it  may  be  possible  to  use  an  “Adapter  Class”  to  provide  a 
conversion. 

Software  prototyping  is  being  used  during  the  study  to  experiment  with  methods  of  allowing 
acoustic  sub-models  to  be  added  without  changing  code  in  a  “plug  &  play”  style.  Prototyping 
has  also  been  used  to  develop  aspects  of  the  user  interface.  Figure  03-F.2  is  an  example 
from  the  prototype  software,  which  is  used  to  display  the  output  of  various  acoustic  sub¬ 
models  showing  seabed  scattering  strength  in  dB,  from  the  McKinney  and  Anderson 
algorithm,  as  a  function  of  grazing  angle  (abscissa  axis)  and  bottom  type  (ordinate  axis).. 


Thrust  04  Tactical  Active  Sonar  (TAS) 

Project  04- A:  Advanced  active  sonar 
(Alliance  days  -31,  Manning  days  -  2) 


Operational  relevance 

To  improve  sonar  performance  in  shallow  water,  reverberation-limited  environments  by 
exploiting  the  right/left  ambiguity  resolving  capability  of  advanced  towed  arrays,  broadband 
source/receiver  technology  and  advanced  processing  techniques. 


Background 

Towed-array  based  active  sonar  (ATAS)  systems  have  been  the  subject  of  evaluation  and 
study  at  the  Centre  since  1981 .  The  major  challenge  is  to  better  discriminate  between  the 
threshold  crossings  caused  by  false  targets  and  true  target  returns.  The  data  contain  more 
target  information  than  can  be  extracted  with  available  algorithms  for  signal/data  processing 
and  sonar  information  extraction.  Future  systems  will  have  the  potential  to  ensure  exploitation 
of  improved  understanding  of  underwater  acoustics,  computer  technology  and  the 
development  of  new  algorithms. 

The  Active  Towed  Array  Sonar  (ATAS)  was  upgraded  by  implementing  a  broadband  R/L 
directional  cardioid  towed  array  comprising  two  octaves,  high  capacity  data  acquisition  and 
programmable  real  time  processing  system  (CDAS).  When  the  broadband,  high  power 
sound  source  is  delivered  in  March  2000,  the  Centre  will  be  well  positioned  to  study  advanced 
system  concepts  for  signal/data  processing  and  information  extraction,  which  in  the  short 
to  medium  term  could  be  exploited  during  the  development  and  subsequent  phases  of 
national  mid-life  conversion  and  upgrade  programmes. 

Research  has  concentrated  on 

•  investigating  the  potential  of  the  cardioid  array 

•  investigating  the  potential  of  broadband  sonar  signals 

•  defining  a  Measure-of-Performance  (MOP)  applicable  to  algorithms  aiming  at  clutter 
reduction 

•  planning  and  executing  the  Joint  Research  Project  (JRP)  sea  trial  MERCURY  99  in 
conjunction  with  the  Defence  Evaluation  and  Research  Agency  (DERA). 

Cardioid  arrays 

Conventional  active  towed-array  sonar  suffers  from  the  problem  of  ambiguity.  It  is  impossible 
to  distinguish  between  returns  from  port  and  starboard.  Although  operational  methods  for 
overcoming  this  problem  exist,  they  are  time  consuming,  cannot  be  used  for  single  ping 
applications  and  do  not  provide  additional  processing  gain  against  reverberation. 
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France,  Germany,  the  Netherlands  and  the  United  Kingdom  are  investigating  equipment 
based  solutions  to  the  problem.  One  proposal  is  to  use  an  array  of  triplet  hydrophones. 
Such  an  array,  also  called  a  cardioid  array,  using  state-of-the-art  sensor/telemetry  technology 
covering  two  octaves  has  been  procured.  The  array  comprises  126  hydrophone  triplets 
implemented  in  a  nested  configuration. 
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Figure  04-A.1(a)  Range-Bearing  Detection 
display  (B-scan)  derived  from  a  conventional 
linear  array  AH  features  are  symmetrically 
distributed  around  the  180°  degree  bearing 
that  corresponds  to  endfire-aft.  The  R/L 
ambiguity  problem  is  evident  04-A.  1(b)  The 
same  data  analyzed  using  the  signals  from 
the  hydrophone  triplets.  Features  have  been 
resolved  to  port  or  starboard.  Important 
features  have  been  labelled. 


Reverberation 

Feature 
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The  decision  to  procure  a  cardioid  array  does  not  imply  that  it  is,  in  principle,  superior  to 
other  concepts  such  as  the  twin  array  approach.  It  was  primarily  specified  as  a  research 
tool  to  be  used  at  low  speeds,  for  investigating  broadband  underwater  acoustics  (clutter/ 
reverberation  up  to  3600  Hz)  and  related  signal/data  processing  issues,  more  efficiently 
than  would  be  possible  with  a  traditional  R/L  ambiguous  array. 

Investigation  of  the  new  technology  and  algorithms,  (one  of  which  can  be  used  with  the 
real-time  system  over  a  broad  range  of  frequencies)  has  been  facilitated  by  the  use  of  data 
acquired  during  BACCHUS’98  (a  JRP  with  DERA)  using  the  DERA  40  triplet  array.  Run 
geometries  between  Malta  and  Sicily  were  combined  with  known  bathymetric  and 
reverberation  features,  which  would  allow  investigation  of  left-right  discrimination. 

The  following  results  from  the  BACCHUS’99  trial  demonstrate  the  versatility  of  the  triplet 
array  which  can  be  used  to  improve  sonar  performance  and  investigate  environmental 
acoustic  problems  such  as  reverberation  and  clutter. 


Figure  04-A.1  shows  typical  B-scans  from  BACCHUS’98.  Figure  04-A.1  (a)  shows  the  range¬ 
bearing  sonar  display  as  it  would  be  derived  from  a  traditional  ambiguous  array.  Features 
are  replicated  to  starboard  (180°-  360°)  and  to  port  (0°-  180°).  The  bow  shaped  returns 
from  a  reverberation  ridge,  a  feature  which  has  been  measured  and  observed  frequently, 
obscures  back-beam  targets.  Figure  04-A.1(b)  shows  the  same  data  after  the  2-step  cardioid 
BF  algorithm  has  been  applied.  The  reverberation  ridge  has  been  resolved  to  starboard  (in 
agreement  with  known  ground  truth)  and  the  bright  features  between  8  and  10  km  have 
been  separated  to  port  and  starboard.  Similarly,  targets  of  opportunity  have  been  resolved 
(denoted  as  target  1  and  target  2). 

A  more  quantitative  measure  of  R/L  ambiguity  resolution  is  the  suppression  ratio  of  the 
echo  level  of  the  true  beam  versus  the  echo  level  in  the  ambiguous  beam.  It  can  be  derived 
by  analyzing  the  received  signal  levels  as  a  function  of  bearing,  i.e.,  by  looking  at  the 
directional  A-scan.  In  Fig.  04-A.2,  three  A-scans  are  shown  containing  returns  from  three 
different  targets  (1  ).Two  of  the  three  targets  can  also  be  detected  in  the  ambiguous  directions 
(2).  In  the  top  scan  the  return  (1)  is  due  to  the  Campovega  oil  platform.  In  the  opposite, 
ambiguous  beam,  the  respective  echo  level  (2)  is  significantly  lower,  resulting  in  a  suppression 
ratio  of  about  12  dB.The  middle  scan  shows  the  reverberation  ridge  clearly  displayed  to 
starboard,  with  no  corresponding  signal  in  the  opposite,  ambiguous  beam.  In  the  third 
scan,  a  bright,  bottom  feature  can  be  seen  to  starboard. The  corresponding  signal  (2)  in  the 
ambiguous  bearing  has  been  reduced  by  about  1 0  dB.  One  extra  feature  has  been  indicated 
in  this  figure:  passive  ship  noise  has  been  resolved  to  starboard. This  highlights  the  possible 
use  of  the  triplet  array  in  passive  mode  as  well  as  for  active  systems  and  demonstrates  the 
advantage  of  R/L  directivity  on  the  receiver  rather  than  on  the  source. 

Higher  suppression  ratios  than  those  mentioned  have  been  observed  for  returns  from  “point¬ 
like”  scatterers  in  specific  scenarios.  It  has  been  possible  to  show  that  this  performance 
can,  by  judicious  choice  of  algorithms,  be  extended  to  scan  angles  up  to  60°  from  broadside. 
The  suppression  ratio  is  not  however  the  only  performance  parameter  which  influences  the 
beam  former  design.  In  the  presence  of  poorly  correlated  noise,  the  resolution  of  R/L 
ambiguity  is  at  the  expense  of  sensitivity.  This  loss  of  sensitivity  increases  with  decreasing 
frequency  or  decreasing  intra-triplet  hydrophone  spacing.  Consequently,  for  some  scenarios, 
such  as  flow  noise  limited  operations,  the  performance  of  R/L  suppression  algorithms  may 
be  severely  degraded.  What  may  be  required  is  a  BF  concept  which  allows  adaptation  to 
the  characteristics  of  the  noise  at  the  triplet  hydrophones  and  to  the  directional  distribution 
of  reverberation/clutter. 

To  further  illustrate  the  processing  gain  against  reverberation,  a  target  track  based  on  the 
same  triplet  hydrophone  data  set  is  shown  in  Fig.  04-A.3  as  a  40-ping  history,  applying  the 
conventional,  linear-ambiguous  beamformer  (a)  and  the  cardioid  beamforming  algorithm 
(b).  The  abscissa  and  ordinate  show  range  and  beam  number  (discrete  bearing  values), 
respectively.  Each  of  the  40  pings  has  been  interleaved  on  a  beam  to  beam  basis.  In 
Fig.  04-A.3(a)  returns  from  the  target  and  the  reverberant  ridge,  although  from  opposite 
directions,  cannot  be  separated.  The  reverberation  returns  almost  mask  the  target  track.  In 
Fig.  04-A.3(b),  the  target  track  and  the  reverberation  ridge  have  been  resolved  to  port  and 
starboard,  respectively.  Consequently,  a  much  brighter  target  track  can  be  observed. 


Broadband  processing 

Sonar  signals  covering  a  broad  frequency  spectrum,  e.g.  1-2  octaves,  have  potential  to 
improve  sonar  performance  under  reverberation  limited  conditions  for  two  reasons.  First, 
large  bandwidths  result  in  very  small  range  resolution  cells.  Application  of  proper  “post¬ 
detection”  integration  reduces  fluctuations  in  target  echo  and  reverberation  level,  resulting 
in  a  better  SNR  at  the  detector  input.  Second,  large  bandwidth  signals  are  suitable  for 
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Figure  04-A.2  Three  examples  of 
directional  A-scans  i.e.  received  signal 
levels  as  a  function  of  bearing  using  data 
from  the  cardioid  array.  Returns  from 
three  targets  ( 1).  Two  of  the  three  targets 
can  also  be  detected  at  lower  levels  in 
the  ambiguous  directions  (2). 
Suppression  ratios  for  target  amplitudes 
in  the  true  versus  the  ambiguous 
bearings  range  from  10  to  12  dB. 


a 


b 
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Figure  04- A. 3  Forty  pings  of  data 
interleaved  to  give  a  history  of  a  target 
track.  Multiple  tracks  in  adjacent  beams 
are  due  to  sidelobe  leakage.  The  data 
were  obtained  for  a  target  ‘ apparently ' 
passing  through  reverberation ,  i.e., 
physically  the  target  and  reverberation 
patches  were  in  a  pair  of  ambiguous 
beams.  Traditional  processing  (Fig  04- 
A.3a)  makes  them  appear  to  be 
physically  co-located.  Only  a  weak  target 
track  can  be  seen  due  to  the  ambiguous 
reverberation  background  overlaying  the 
target  returns,  (b)  Using  the  triplet  array 
the  reverberation  returns  and  the  target 
echoes  are  separated  into  their 
respective  beams  and  the  target  can  be 
seen  clearly  over  the  entire  track. 
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applying  sub  band  processing  techniques,  which  exploit  the  frequency  dependence  of 
detection  performance. There  are  indications  that  detection  performance  (e.g.  SNR)  across 
frequency  bands  varies  significantly  with  time. The  reason  for  this  is  that  many  of  the  terms 
in  the  sonar  equation  such  as  transmission  loss,  boundary  backscattering  strength,  true/ 
false  target  strengths,  scattering  statistics  and  transmission  loss  statistics,  depend  on 
frequency,  time  and  space. The  challenge  is  to  design  a  fusion  concept  which  combines  the 
SNR  from  all  frequency  bands  in  such  a  way  that  the  resulting  global  detection  performance 
is  maximized. 


Fig04-A-4  AF  =  120Hz 
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Figure  04-A.4  Matched  filter  output  signals  of  five 
sub-bands  arranged  with  increasing  centre  frequency 
from  bottom  FI  to  top  F5 .  Considerable  fluctuations 
in  echo  shape  and  amplitude  as  a  function  of  space, 
time  and  frequency  are  visible  and  suggest  that  the 
underlying  statistics  may  be  used  for  choosing  an 


Broadband  (1200  Hz)  linear  FM  signals  of 
12  s  duration  centred  around  3  kHz  were 
transmitted.  A  sub-band  matched  filter 
processing  scheme  was  applied  to  divide 
the  band  into  ten  sub-bands.  Individual  echo 
signals  from  5  different  sub  bands  (Fig.  04- 
A.4)  and  the  signal-to-noise  ratios  for  each 
sub-band  signal  (Fig.  04-A.5)  are  shown  for 
a  sequence  of  50  and  140  pings, 
respectively.  There  is  no  single  sub-band  in 
which  the  average  signal  shape  and  SNR  is 
distinctly  different  from  other  sub-bands. 
Ping-to-ping  comparison,  however,  reveals 
rapid  changes  of  echo  shape  and  SNR. The 
detection  scheme  will  depend  on  the 
underlying  fluctuation  statistics.  The  final 
scheme  will  be  selected  after  the  features 
corresponding  to  false  targets  have  been 
examined. 

Another  advantage  of  FM  sub-band 
processing  is  that  the  relative  velocity 
between  source  and  receiver  (range  rate) 
can  be  estimated  based  on  a  single  ping 
observation,  i.e.  without  ping  history  which 


Figure  04-A.5  SNR  as  a  function  of  the  centre  frequency  for  10  consecutive  120  Hz  sub  band  signals  in  a  reverberation 
limited  environment.  For  each  sub-band,  the  average  SNR  for  the  entire  run  is  indicated.  Although  the  average  SNR  is 
not  differing  much,  it  is  not  obvious  that  simple  cross  band  averaging  is  the  best  detector  strategy  especially  around 
the  low  SNR  situation  around  ping  100. 
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Figure  04-A.6  Normalized  intra-ping  sub  band 
echoes  signals  with  SNR  estimates  for  each  sub¬ 
band  indicated  on  the  right  of  each  signal.  The 
differential  time  delay  between  the  sub-bands  is 
shown.  The  top  plot  represents  the  normalized 
average  of  the  individual  outputs.  Note  that  the 
echoes  from  different  sub  bands  are  uncorrelated. 
SNR  differences  of  up  to  10  dB,  the  Doppler-induced 
differential  delays  and  the  relatively  low  SNR  of  the 
sum  signal  (top  of  this  figure)  suggest  detector 
schemes  different  from  simple  averaging. 


would  allow  range  rate  estimates  from 
inter-ping  delay  measurements.  Figure 
04-A.6  shows  the  decreasing  intra-ping 
delay/travel  times  between  the  echoes 
from  the  lowest  and  the  highest  sub-band 
received  during  the  closing  part  of  the 
run.  This  is  explained  by  combining  the 
facts  that  the  high  frequencies  are 
transmitted  last  (up-sweep  pulse)  and 
that  during  transmission  time  the  target 


TIME  (sec) 

moves  towards  the  receiver.  Quantitative  analysis  of  the  measured  delays  has  to  take  into 
account  the  well  known  Range-Doppler  ambiguity  of  the  FM  pulse.  It  has  been  shown  that 
there  exists  a  mathematical  relationship  between  the  true  range  rate  and  the  sub  pulse 
time  delays.  Figure  04-A.7  presents  the  comparison  between  the  actual  (inter-ping)  and 
the  intra-ping  rate  both  without  (a)  and  with  (b)  Doppler  compensation.  Without  Doppler 
compensation  the  range  rate  which  is  derived  from  the  differential  intra-ping  time  delays 
differs  from  the  nominal  value  by  a  factor  of  2. 


Figure  04-A.6  reveals  three  other  interesting  facts:  First,  the  echo  signatures  in  each  sub 
band  are  uncorrelated.  Second,  the  S/N  ratio  is  varying  by  as  much  as  10  dB.  Third,  simple 
averaging  of  the  sub  band  signals  (top  signal  in  Fig.  04-A.3)  provides  a  SNR  of  25  dB, 
which  is  much  less  than  the  best  SNR  of  the  sub  band  signals.  This  suggests  two  detector 
strategies:  k  out  of  n  rule  or  sub  band  averaging  after  motion  compensation  of  the  differential 
delaytimes. 


Figure  04-A.7  Comparison  between  inter-ping  (red)  and  intra-ping  range  rate  estimation  without  (Fig.  4-A.7a)  and 
with  (Fig.  04- A.  7b)  Doppler  compensation.  Range  rate  estimation  from  sub  band  delay  time  measurements  within  one 
ping  is  feasible,  but  requires  a  Doppler-independent  correction  factor. 
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Figure  04-A.8  SNR  for  a  sequence  of  40  target 
echoes  as  simultaneously  observed  in  two 
different  frequency  bands.  The  data  were  taken 
during  the  BACCHUS  98  trials  and  demonstrate 
that  the  lower  band  delivers  SNRs  which  are 
on  average  about  3  dB  higher  than  in  the  higher 
band.  Both  frequency  bands  are  100  Hz  wide 
and  separated  by  500  Hz. 


An  example,  which  may  prove  that  detection  performance  is  frequency  dependant  was 
taken  from  BACCHUS  98  data.  Figure  04-A.8  depicts  the  SNR  for  a  sequence  of  40  echo 
signals  processed  at  two  frequency  bands  100  Hz  wide  and  separated  by  500  Hz.  The 
lower  frequency  band  provides  on  average  3  dB  more  SNR,  although  the  reverberation 
patch  size  is  larger  than  for  the  higher  frequency  band.  A  detection  scheme  for  optimally 
combining  the  detection  results  from  both  bands  is  still  a  challenge.  The  solution  again  will 
depend  on  the  in-band  SNR  fluctuation  statistics.  More  data  have  to  be  analyzed  to  better 
determine  these  characteristics. 
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Figure  04-A.9  Unnormalized  Range-Bearing  sonar  image  generated  from  the  receive  signal  amplitudes  of  each 
beam  (left),  the  same  data  set  after  normalization  (middle),  the  same  data  set  after  thresholding  (right).  Rectangles 
arelxl  km.  This  series  of  displays  allows  an  operator  to  classify  alarms  by  reviewing  the  normalized/unnormalized 
sonar  displays  to  determine  if  a  contact  could  be  due  to  e.g.,  a  reverberation  edge,  similar  to  the  bottom  four  contacts 
in  the  right  hand  image. 


Measure  of  performance  (MOP) 

Algorithms  which  are  implemented  for  the  detection  of  active  sonar  echoes  in  antisubmarine 
warfare  are  evaluated  against  measure  of  performance. The  classical  ROC  curve  measures 
probability  of  detection  against  probability  of  false  alarm  per  range  cell.  For  highly  non¬ 
stationary  data  or  after  application  of  a  cluster  algorithm,  the  probability  of  false  alarm  per 
bin  becomes  less  meaningful. 

We  have  defined  a  new  measure  of  performance  which,  in  a  ping  of  data,  makes  objects  of 
connected  pixels.  The  aiming  point  of  the  object,  which  is  usually  the  data  mass-center, 
defines  the  estimated  location  of  the  object.  An  object  within  a  certain  distance  from  the 
submarine  is  considered  a  detection,  otherwise  it  is  counted  as  a  false  alarm.  The  number 
of  false  alarms  is  now  per  ping. 
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During  the  evaluation  phase  of  the  MOP  a  very  useful  graphical  tool  for  classifying  alarms 
was  implemented.  The  underlying  concept,  which  has  been  used  successfully  for  some 
time,  is  explained  in  Fig.  04-A.9.  The  raw  sonar  image  generated  from  the  receive  signal 
amplitudes  of  each  beam  (left),  the  same  data  set  after  normalization  (middle),  the  same 
data  set  after  thresholding  (right).  The  images  from  left  to  right  contain  progressively  less 
information.  The  operator  is  alerted  by  the  alarms  in  the  easy-to-scan  detection  display  at 
the  right  and  can  revisit  the  preceding  images  concentrating  on  an  area  (e.g.  1  x  1  km) 
around  each  alarm.  Thus  he  can  regain  information  needed  for  classification  which  has 
gradually  been  removed  during  the  process  of  generating  the  detection  display.  For  the 
data  shown  in  Fig.  04-A.9  the  operator  would  conclude  that  the  bottom  four  contacts  could 
be  residuals  from  the  returns  from  a  reverberation  ridge.  The  system  also  allows  the  operator 
to  examine  the  relevant  sections  of  the  original  A-Scan  signals. 


Page  test 

The  Page  test  algorithm  was  designed 
to  detect  the  beginning  and  ending  of  a 
signal  in  a  stationary  background.  The 
Page  test  has  not  yet  been  compared 
with  the  standard  threshold  detector  on 
a  large  real  data  set,  but  it  has  been 
compared  to  the  threshold  detector  for 
real  data,  consisting  of  164  pings  of 
active  FM  transmissions,  using  the 
measure  of  performance  described 
above.  The  result  can  be  seen  in  Fig.  04- 
A-IO.The  red  curve  shows  the  Page  test 
detector,  the  blue  curve  the  threshold 
detector.  The  Page  test  shows  superior 
performance  over  a  large  region  of 
sensible  parameter  values.  For  example, 
at  200  false  alarm  objects,  the  Page  test 
detector  has  20%  higher  probability  of 
detection  than  the  threshold  detector. 


Figure  04- A.  10  Comparison  of  Page  test  detector,  Markov 
random  field  segmentation  and  threshold  detector.  The  red 
curve  relates  to  the  Page  test  detector,  the  blue  curve  to  the 
threshold  detector  and  the  symbols  to  the  Markov  random 
field  segmentation.  The  Page  test  shows  superior 
performance  over  a  large  region  of  sensible  parameter 
values.  The  Markov  random  field  segmentation  results  are 
between  the  Page  test  detector  and  the  threshold  detector. 


The  interplay  between  the  measure  of  performance  (object  counting)  and  the  Page  test  are 
shown  in  Fig.  04-A.1 1  .The  fragmented  reverberation  returns  belonging  to  the  same  bottom 
feature,  cause  many  detections  in  the  threshold  detector,  but  are  seen  to  be  combined  by 
the  Page  test  into  one  detection. 


Interbeam  clutter  reduction 

Operating  a  low  frequency  active  sonar  in  shallow  water  results  in  a  large  number  of  target 
like  clutter  returns.  An  automatic  method  of  image  segmentation  using  Markov  Random 
Field  modelling  has  been  used  to  reduce  clutter.  The  method  looks  at  detections  over 
range  and  bearing.  It  removes  small  objects  which  do  not  exhibit  the  right  signature  over 
beams.  Separate  detections  corresponding  to  one  large  object  are  combined  to  form  one 
single  object.  Objects  too  large  to  be  a  submarine  can  then  be  removed. 

The  Markov  Random  Field  model  used  is  based  on  the  physical  and  probabilistic  knowledge 
of  the  sonar  picture.  It  assumes  that,  statistically,  a  target  has,  on  average,  a  larger  SNR 
and  that,  on  a  local  scale,  the  sonar  display  exhibits  homogeneity.  The  model  is  tuned  to 
remove  as  much  clutter  as  possible  while  retaining  the  target. 


—  68  — 


Beam  number  Beam  number 


Figure  04- A.  1 1  Effect  of  Page  test  detection  (at  the  right)  compared  to  threshold  detector  (at  the  left)  on  a  bottom 
feature.  The  threshold  detector  has  fragmented  the  bottom  feature  into  many  separate  detection  objects.  The  Page 
test  detector  interprets  them  as  one  single  object. 

In  Fig.  04-A.10the  result  of  the  Markov  random  field  segmentation,  denoted  by  the  symbols, 
is  compared  to  the  Page  test  and  threshold  detector.  For  the  single  ping  case,  performing 
Markov  random  field  segmentation  has  no  advantage  over  just  performing  the  Page  test 
detector. 

The  Markov  random  field  segmentation  algorithm  was  designed  to  merge  separate  returns 
from  one  geophysical  object.  As  the  detection  of  large  geophysical  objects  is  irreconcilable 
with  the  detection  of  a  small  submarine,  use  of  a  single  algorithm  is  not  indicated.  The 
Markov  random  field  segmentation  can  still  be  used  to  detect  large  geophysical  objects, 
which  can  be  removed,  providing  that  size  or  geological  position  do  not  change. 


MERCURY  99 

The  JRP  with  DERA  Broadband  Active  Multistatic  Sonar,  shares  equipment,  resources, 
technology  and  acoustic  data  for  the  purpose  of  investigating  issues  related  to  broadband 
active  sonar,  right/left  directional  receive  arrays  and  active  multistatics.  A  series  of 
collaborative  at-sea  experiments  are  being  conducted  with  one  major  annual  sea  trial.  The 
objective  of  MERCURY  ’99  was  to  evaluate  broadband  transmissions,  technology  and 
processing  techniques. 


SACLANTCEN  technicians  installing 
transducers  and  hydrophones  on  the 
Italian  submarine  Sauro  in  Agusta  for 
the  Joint  Research  Project, 
Broadband  Active  Multistatic  Sonar 
experiment 
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The  experiment  was  conducted  primarily  in 
the  Malta  Plateau,  south  of  Sicily. The  water 
depth  in  this  area  is  100-200  m,  with  high 
levels  of  reverberation  and  clutter. 
Participating  assets  included  the  R/V 
Alliance  (broadband  acoustic  source,  two 
broadband  towed  receive  arrays  and 
command  and  control  functions),  the  ITS 
Tavolara  (echo  repeater)  and  the  ITS  Sauro 
(submarine  target).  DERA  provided  the 
broadband  acoustic  source  and  a 
broadband  receive  array  (including  real-time 
processor).  SACLANTCEN  provided  the 
broadband  R/L  directional  receive  array 
(including  real-time  processor),  echo 
repeater. 

Data  was  recorded  with  broadband  (1000- 
3400  Hz)  active  transmissions  on  two 
separate  broadband  towed  arrays.  Figure 
04-A.12  shows  data  from  a  single 
transmission.  The  broadband  phone  data 


from  the  SACLANTCEN  array  were  processed  over  contiguous  frequency  bands  of  500  Hz 
bandwidth.  One  hundred  and  sixty  beams  are  formed  from  forward  to  aft  and  the  smoothed 
(mean)  power  levels  are  shown  out  to  35  km  with  the  same  relative  dB  power  scale.  Looking 
at  the  overall  colour  distribution  one  recognizes  features  which  are  different  for  each 
frequency  band,  such  as  the  background  noise  level  and  the  reverberation  returns  from  the 
Ragusa  Ridge  at  around  15  km  and  beams  65-80.  Diffuse  reverberation,  immediately 
following  the  direct  blast  with  similar  intensity  and  decay  rate  show  little  dependence  over 
the  four  frequency  bands.  Other  passive  and  clutter  features  are  also  evident  in  these  data. 
Data  such  as  this  will  allow  evaluation  of  the  frequency  dependence  of  some  of  the  most 


prominent  parameters  of  the  sonar  equation. 


The  photographs  on  this  page 
show  the  hydrophones 
installed  on  the  Italian 
submarine  Sauro 
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chaos  theory,  he  moved  to  DERA  (Malvern)  in  1993,  where  he  worked  on 
advanced  algorithms  for  radar  and  sonar  applications.  Since  1997  he  has 
worked  at  SACLANTCEN  in  the  areas  of  signal  and  information  processing. 


Rene  Laterveer  received  the  Ph.  D.  in  theoretical  physics  from  the 
University  of  Amsterdam  in  1993  on  a  subject  in  elementary  particle 
physics.  From  1992  to  1995  he  was  at  TNO  Physics  and  Electronics 
Laboratory  in  the  Hague,  the  Netherlands,  working  on  active  low  frequency 
sonar.  He  has  been  a  SACLANTCEN  since  1996. 


Marcel  van  Velzen  received  a  masters  degree  in  theoretical  physics 
from  the  University  of  Amsterdam  in  1987  and  until  1989  worked  at  the 
Netherlands  National  Institute  for  High  Energy  Physics  (NIKHEFH).  From 
1 989  to  1 998  he  was  at  the  TNO  Physics  and  Electronics  Laboratory  in 
the  Hague,  the  Netherlands,  where  he  worked  on  Synthetic  Aperture 
Satellite  Radar  processing  and  from  1991  worked  on  real-time 
processing  and  data  analysis  for  the  Netherlands  Low  Frequency  Active 
Sonar  program.  In  1999  he  started  at  SACLANTCEN  as  a  senior  scientist 
working  in  the  areas  of  signal  processing  and  data  analysis  related  to 
Active  Sonar. 
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Project  04-B.  Deployable  Underwater  Surveillance 
Systems  (DUSS) 

(Alliance  days  -16,  Manning  days  -  4) 


Operational  relevance 

An  active  sonar  concept  comprising  a  network  of  small,  autonomous  sources  and  receivers 
designed  to  optimize  multistatic  operation  against  small  targets  in  shallow  and  coastal 
waters,  with  heavy  shipping  traffic  and  strong  reverberation. 


Concept  demonstration  and  performance  study 

The  following  goals  were  achieved: 


•  Demonstration  of  the  DUSS  concept  at  sea. 

•  System  performance  estimation  of  monostatic  and  multistatic  receivers  in  a  wide 
range  of  conditions. 

•  Assessment  of  tradeoffs  between  performance  and  sonar  resources. 

•  Assessment  of  coverage  extension  via  additional  multistatic  receivers  and  target 
aspect  diversity. 

Echo  fluctuations  and  diversity 

Significant  random  fluctuations  of  echoes  from  real  targets  and  the  Echo  Repeater  were 
observed,  using  CW  and  FM  pulses.  A  statistical  analysis  of  echo  signal-to-noise  ratio 
(SNR)  was  effected,  using  calibrated,  aspect-independent  returns  from  the  Echo  Repeater. 
Fluctuations  were  decomposed  into  a  slowly-evolving,  low-pass  component  (Fig.  04-B.1) 
and  a  white  random  process.  (Fig.  04-B.2).The  latter  component  strongly  affects  detection 
performance  of  individual  DUSS  receivers.  It  was  demonstrated  to  be  uncorrelated  in  space 
(buoy-to-buoy),  time  (ping-to-ping)  and  frequency  (CW  pings  in  adjacent  bands).  The 


PING  number 

Figure  04-B.  1  Long-term  components  of  signal 
fluctuations. 


Figure  04-B. 2  Short-term ,  uncorrelated  components  of 
signal  fluctuations 
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Figure  04-B.3  A  zoom  view  of  echo  signals  shows  that  multipath  effects  come  and  go  from  ping  to  ping. 


Figure  04- B. 4  A  comparison  between 
measured  propagation  loss  (red)  and  C- 
SNAP  model  (black)  shows  the  impact  of 
propagation  loss  patterns  and  target 
movements  on  signal  fluctuations 
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3  small  receivers 
3  x25  hydrophones 
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Single  large  receiver 
75  hydrophones 

Single  small  receiver 
1  x  25  hydrophones 


Figure  04-B.5  System  detection  performance  enhancement  when  the  receivers  are  integrated  into 
a  multistatic  system. 
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Fraction  of  collected  data  samples 


relationships  between  multipath  (Fig.  04-B.3)  and  measured  versus  predicted  propagation 
loss  (Fig.  04-B.4)  were  analyzed  and  detection  performance  improvements  were  quantified 
in  terms  of  equivalent  SNR  gain  and  ROC  curves  (Fig.  04-B.5).  Such  “diversity”  gain  is 
additional  to  the  typical  “aspect  diversity”  of  multistatic  sonar  with  real  targets. 


Contact  consistency,  localization,  fusion,  tracking 


Run  01 03  Receiver  n.2  Acoustic  vs  DGPS  target  position 


Run  0103,  Rec.  2,  Empirical  distrib.  of  localiz.  errors 


Figure  04-B.6  Contact  localization  consistency  of 
DUSS  networks  was  analyzed  as  a  function  of  system 
characteristics. 


The  DUSS  network  of  deployed  transmitters 
and  receivers  integrates  and  reduces  acoustic 
data  into  a  comprehensive,  coordinated 
tactical  picture  of  the  monitored  waters.  This 
is  possible  because  the  contacts  from  the 
buoys  are  geometrically  consistent  (Fig.  04- 
B.6).  Contact  fusion  is  therefore  possible,  with 
enhanced  detection  performance  and  target 
localization  performance.  The  advantages  of 
multistatic  systems  were  demonstrated  by 
comparing  passive  listening  (bearings  only), 
active  sonar  (time  delay)  and  autonomous 
underwater  vehicle  support  (AUV)  (pinger 
pulse  delay)  (Fig.  04-B.7).  The  option  of  in¬ 
buoy  tracking  and  data  reduction  prior  to  data 
fusion  is  being  studied. The  evident  advantage 
consists  of  reduced  and  therefore  more  robust 
data  flow  between  buoys.  The  study  will 
quantify  overall  detection  performance, 
diversity  gain  and  localization  precision. 


Development  of  a  new  DUSS 
prototype  receiver. 

A  prototype  receiver,  designed  and 
constructed  at  the  Centre  completed 
engineering  tests  in  December.  The  system 
architecture  definition  has  been  finalized.  A 
prototype  was  deployed  during  SIRENA  ’99 
and  GEOSCAT  ’99. 


Bearings-only  Time-only 


Figure  04-B.7  Contact  localization  performance  of  DUSS  networks  in  passive  (bearings-only)  mode  and  in  active 
sonar  (time-only)  mode.  Colour  maps  show  expected  localization  errors  versus  target  position. 
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Elevated  communications 
rotoy  station  (o.g.  tethered 
biBoon.  antenna  ashore, 
support  ship) 


Local  Area  Networks 


Figure  04-B.8  Radio  network  techniques 
implement  the  “sonar  Internet’  of  DUSS 
units  for  experiments. 


The  system  parameter  tradeoffs  and  environmental  factors  affecting  long-range,  digital 
radio  links  were  analyzed  with  the  support  of  E/M  propagation  modelling. 

A  theoretical  study  was  carried  out  on  the  applicability  and  potential  of  spread-spectrum 
network  techniques  to  DUSS  buoys.  Sub-systems  were  procured  in  order  to  implement 
tests  at  sea  and  validate  the  model-based  predictions. 

Specification ,  bid  and  procurement  of  DUSS  transmitter  and 
receivers. 

An  advance  performance  specification  has  been  forwarded  to  potential  bidders.  The 
specification  will  be  finalized  by  the  end  of  the  year.  Bids  will  be  requested  by  January  2000 
and  offers  will  be  finalized  by  July  2000.  Contract  award  will  be  completed  by  summer  2000. 
A  detailed  plan  of  delivery  and  acceptance  tests  of  the  sonar  units  is  in  preparation. 

The  acquisition  and  real-time  processing  system  were  defined.  Computer  hardware 
architecture  and  specifications  were  finalized  and  the  procurement  process  was  initiated. 
The  software  architecture  and  specifications  were  defined.  Key  modules  of  the  existing 
system  were  tested  on  an  Alpha  workstation. 

.4  Average  Bit  Error  Rate  2.28  GHz,  7.8m  33m 
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Figure  04- B. 9  Average  bit  error  rate 
obtained  with  digital  radio  link  during  10* 
DUSS’98  tests. 
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Project  04-B  publications  and  presentations 

Berni,  A.,  Mozzone,  L.The  application  of  spread  spectrum  communications  to  REA  tactical  networks  and  deployable 
surveillance  systems,  SACLANTCEN  SM-367. 

Mozzone,  L.,  Bongi,  S.  Diversity  in  deployable  underwater  surveillance  systems,  SACLANTCEN  SR-31 8. 

Mozzone,  L.,  Bongi,  S.  Diversity  in  multistatic  active  sonar.  In:  IEEE  Oceans‘99  Conference,  1999,  Seattle,  US. 

Mozzone,  L.,  Guerrini,  P.,  Berni,  A.  Long  range,  large  throughput  radio  data  link  for  DUSS  (Deployable  Underwater 
Surveillance  Systems),  SACLANTCEN  SM-360. 

Mozzone,  L.,  Lorenzelli,  P.,  Bongi,  S. Target  localization  with  multiple  sonar  receivers,  SACLANTCEN  SR-317. 


Lorenzo  Mozzone  received  the  Masters  degree  in  electronic  engineering 
from  the  University  of  Genova  in  1989.  He  worked  for  four  years  with  Alenia 
S.p.a.,  Genova  and  has  been  at  SACLANTCEN  since  1 994.  His  activity  is 
focused  on  sonar  systems  research  with  specific  reference  to  deployable 
underwater  surveillance  systems  (DUSS).  He  has  published  papers  on 
signal  processing  techniques ,  real  time  implementation  and  operational 
research  studies  on  activated  towed  arrays. 
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Project  04-C:  Low  frequency  shallow  water 
reverberation  and  propagation:  adaptation  to  large 
bandwidths 
(Alliance  days  -  25) 


Operational  relevance 

Low  frequency  (1-4  kHz)  active  sonar  systems  in  shallow  water  must  contend  with 
seafloor  reverberation,  which  limits  and  confounds  the  ability  to  detect  and  classify 
small  submarines.  NATO  systems  under  development  are  beginning  to  rely  on  frequency 
diversity,  which  recognizes  that  the  frequency  of  optimum  system  performance  varies 
according  to  environment.  Which  frequency  operates  best  in  which  environment  is  a 
function  of  a  number  of  parameters  including  the  bottom  scattering  function,  which 
remains  poorly  understood.  The  current  NATO  ‘database’ for  bottom  scatter  is  a  single 
coefficient,  independent  of  frequency  and  area.  One  of  the  goals  of  this  project  is  to 
develop  experimental  techniques  for  measuring  bottom  scattering  in  shallow  water, 
which  in  conjunction  with  advances  in  modelling  will  improve  performance  prediction 
of  frequency  dependent  sonar  systems. 


One  of  the  important  issues  in  bottom  scattering  is  to  know  where  the  scattering  is  coming 
from,  that  is  whether  the  scattering  occurs  at  the  water-bottom  interface  or  from  within  the 
sediment.  The  ability  to  distinguish  the  scattering  location  has  important  implications  for 
performance  prediction,  model  and  database  development  and  the  validation  of  rapid 
environmental  assessment  techniques. 


Results 

A  key  result  in  1999  was  the  development  of  the  capability  to  distinguish  the  origins  of 
scattering.  Figure  04-C.1  shows  results  from  a  shallow  water  area  north  of  the  Island  of 
Elba  (designated  Site  S).  Measurements  at  1800  Hz  and  3600  Hz  are  shown  in  Figs.  04- 
C.la  and  1b  respectively.  The  vertical  axis  represents  vertical  angle  and  the  grazing  angle 
of  beams  13  to  22  is  from  5°  to  65°. 

Corresponding  model  results  presented  in  Fig.04-C.1c  and  04-C. Id  predict  the  scattered 
field  of  water-bottom  interface  scattering.  At  3600  Hz,  the  arrival  structure  (Fig.  04-C.1d)  is 
similar  to  that  observed  in  the  data  (Fig.  04-C.1b).This  result  indicates  that  the  scattering 
at  3600  Hz  occurs  in  approximately  the  upper  10  m  of  the  seabed. 

At  1800  Hz  (Fig.  04-C.1a),  it  is  apparent  that  there  is  a  later  scattering  arrival  (starting  at 
0.08  sec  in  the  highesl/steepest  beam),  not  seen  in  the  modelled  result  (Fig.  04-C.1c). 
From  timing  considerations,  this  late  arrival  must  correspond  to  a  scattering  horizon  at 
about  25  m  below  the  water-sediment  interface. 

Independently  acquired  geoacoustic  data  at  this  site  show  a  sub-bottom  feature  at  24  m. 
Figure  04-C.2  shows  the  sound  speed  structure  versus  depth  in  the  bottom  obtained  from 
a  high-resolution  geoacoustic  analysis.  It  shows  a  randomly  layered  seafloor;  the  low  speed 
layers  are  silt-clay  and  the  high-speed  layers  are  sand-shell  matrix.  Figure  04-C.3  shows 
shells  recovered  from  a  10  cm  layer  at  almost  5  m  depth  sub-bottom.  The  shell  layer  at 
25  m  sub-bottom  gives  rise  to  the  scattering  observed  in  Figure  04-C.1a. 
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Figure  04-C.1e  shows  mode!  predictions  at  1 800  Hz  including  this  deeper  layer;  a  comparison 
with  the  measurements  in  Fig.  04-C.1  a  appears  to  confirm  its  existence  and  general  features. 
Note  that  the  model  shows  a  time  separation  between  the  scattering  from  each  interface, 
whereas  the  measured  data  do  not.This  is  interpreted  as  further  evidence  that  the  scattering 
along  the  monostatic  branch  is  caused  by  the  shell-sand  layers  at  about  0.5  m,  5  m  and 
15  m.  Scattering  from  the  sub-bottom  layer  at  24  m  is  not  observed  at  3600  Hz  (Fig.  17b) 
presumably  due  to  increased  attenuation  and/or  reduced  transmissivity  through  the  overlying 
sedimentary  layers. 


Time  (s) 

Figure  04-C.1  Comparison  of  measured  data  at:  a)  1800  Hz  and  b)  3600  Hz  with  modelled  results  for  water-sediment 
interface  scattering  c)  1800  Hz,  d)  3600  Hz  and  e)  1800  Hz  modelled  results  with  water-sediment  interface  scattering 
and  interface  scattering  at  25  m  sub-bottom.  The  lettering  sequence  follows  a  left  to  right,  top  to  bottom  order.  The 
quantity  plotted  is  received  level  minus  source  level. 


These  data  can  be  processed  to  yield  a  scattering  strength  versus  angle  as  shown  in 
Fig.  04-C.4  (red  curve).  Also  shown  in  the  figure  are  bottom  scattering  data  at  another  site 
about  10  km  away  (blue  curve)  which  show  a  markedly  different  angular  and  frequency 
dependency.  Site  S  bottom  scattering  strength  data  (red  curve)  are  strongly  dependent 
upon  frequency.  The  monotonically  increasing  frequency  dependence  is  perhaps  due  to 
the  fact  that  the  mean  scatterer  (sub-bottom  shell  and  coral  fragments)  dimensions  are 
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Figure  04-C.2  Site  S  sediment  parameters.  Parameters  from  Holland  and  Osier  (1998) 
(solid  line)  and  that  used  in  the  model  predictions  of  Figure  04-C.  1  (dotted  line). 


—  82  — 


Figure  04-C.3  Shell  material  corresponding  to  high-speed  layer  in  Figure  04- 
C.2  near  5  m  depth  sub-bottom. 
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Fig  ure  04-C.4  Scattering  strength  results  at  two  sites:  silty  clay  with  interstitial  sandy  layers 
(red)  and  sandy  with  thin  silt  layer  at  interface  (blue).  Also  shown  is  the  NATO  standard 
(dotted  line). 


1800  Hz 


Figure  04-C.5 

results  at  site  2. 


Scattering  strength 


Grazing  Angle  (deg) 

less  than  a  wavelength. The  observed  angular  dependence  is  flatter  than  the  NATO  standard, 
particularly  below  3600  Hz;  this  may  also  be  due  to  sub-bottom  scattering.  Overall  scattering 
strengths  are  less  than  or  equal  to  the  NATO  standard. 


Figure  04-C.5  shows  the  measured  scattering  strength  data  at  1800  Hz  at  Site  S  with  two 
modelling  results.  The  red  dots  indicate  model  predictions  using  interface  scattering.  The 
blue  dots  include  both  interface  and  sub-bottom  scattering.  Thus,  low  angle  scattering  at 
1800  Hz  is  strongly  influenced  by  sub-bottom  scattering:  in  this  area  down  to  25  m  sub¬ 
bottom. 
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Project  04-C  publications  and  presentations 

Holland,  C.,  Hollett  R.,  Troiano  L.  Bottom  scattering  measurements  in  shallow  water,  (i)  SACLANTCEN  SM-320,  (ii) 
Journal  of  the  Acoustical  Society  of  America. 

Holland,  C.W.,  McDonald,  E.  Shallow  water  reverberation  from  a  time  reversed  mirror,  SACLANTCEN  SR-326. 
Osier,  J.C.,  Algan,  O.  A  high  resolution  seismic  sequence  analysis  of  the  Malta  Plateau,  SACLANTCEN  SR-31 1 . 

Osier,  J.C.,  Algan,  O.,  Holland,  C.W.  A  high  resolution  seismic  sequence  analysis  of  the  Malta  Plateau.  EOS 
Transactions,  AGU,  1999. 

Osier,  J.C.,  Gualdesi,  L.,  Michelozzi,  E.,  Tonarelli,  B.  Piston  coring  capabilities  at  SACLANTCEN:  minimizing  and 
assessing  core  disturbance,  SACLANTCEN  SR-321 . 
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Charles  W.  Holland  received  the  MS  and  Ph.  D.  degrees  in  acoustics 
from  the  Pennsylvania  State  University  in  1985  and  1991  respectively,  in 
1985,  he  began  working  for  Planning  Systems  Inc.,  Virginia,  on  various 
projects  including  geoacoustic  modelling,  seafloor  classification 
techniques,  high  frequency  seafloor  acoustic  penetration  and  low  to  mid 
frequency  bottom  loss  and  bottom  scattering  measurement/modelling 
techniques.  One  of  the  models  which  he  developed,  treating  reflection 
from  a  stochastic  layered  seafloor,  is  employed  in  the  AESS  NATO 
performance  prediction  systems.  From  1995-1 996 he  served  as  Chairman 
of  the  Critical  Sea  Test  Bottom  Interaction  Working  Group ,  a  consortium 
of  scientists  from  universities,  government  laboratories  and  private  industry. 
Since  1996  he  has  been  a  project  leader  at  SACLANTCEN,  leading 
research  activities  in  shallow  water  low  frequency  propagation  and 
reverberation. 


B.  Edward  McDonald  received  his  Ph.  D.  in  physics  in  1970  at  Princeton, 
after  which  he  joined  the  Naval  Research  Laboratory  (NRL)  Plasma  Physics 
Division.  From  1970  to  1980  he  carried  out  numerical  investigations  of 
ionospheric  plasma  processes  related  to  high  altitude  nuclear  weapons 
effects  and  to  naturally  occurring  plasma  turbulence  affecting  satellite 
communication.  From  1980  to  1990  he  worked  for  the  Naval  Ocean  R&D 
Activity/ Naval  Oceanic  -  Atmospheric  Research  Laboratory  developing 
theory  and  numerical  solution  techniques  for  fluid  dynamics  and  nonlinear 
acoustics.  He  joined  the  NRL  Acoustics  Division  in  1990,  where  he 
developed  theory  and  computer  models  for  prediction  and  interpretation 
of  ocean  experiments.  In  1997 he  joined  SACLANTCEN  where  he  carries 
out  theoretical  and  numerical  investigations  related  to  acoustic 
oceanography,  acoustic  modelling  and  fluid  dynamic  processes  affecting 
ocean  acoustics.  McDonald  has  published  papers  in  the  fields  of  solar 
physics,  fluid  dynamics,  plasma  physics,  numerical  analysis,  oceanography 
and  ocean  acoustics.  He  has  been  a  member  of  the  Acoustic  Oceanography 
Technical  Committee  and  the  Underwater  Acoustics  Technical  Committee 
of  the  American  Acoustical  Society.  He  holds  NRL  publication  awards, 
NRL  Division  and  Directorate  best  product  awards  and  is  a  Fellow  of  the 
Acoustical  Society  of  America. 


John  C.  Osier  was  bom  in  Montreal,  Quebec,  Canada  on  July 9, 1964.  He 
received  an  honours  B.Sc.  degree  in  solid  earth  geophysics  from  McGill 
University,  Montreal,  in  1986  and  the  Ph.D.  degree  in  geological 
oceanography  from  Dalhousie  University,  Halifax,  in  1993.  He  was  a  visiting 
fellow  at  the  Defence  Research  Establishment  Atlantic,  Dartmouth,  from 
1993  to  1996,  where  he  conducted  research  in  ocean  bottom  seismology 
and  environmental  acoustics.  In  1996,  he  joined  the  SACLANTCEN  as  a 
scientist  in  the  Large  Scale  Acoustics  and  Oceanography  Group,  where 
he  has  been  involved  in  a  variety  of  projects  in  the  area  of  seabed  acoustics. 
His  research  interests  include  techniques  for  determining  geoacoustic 
properties  of  marine  sediments,  sub-critical  penetration  of  acoustic  energy 
into  the  seabed  and  Quaternary  marine  geology.  Dr.  Osier  is  a  member  of 
the  American  Geophysical  Union,  the  Acoustical  Society  of  America,  the 
Canadian  Geophysical  Union  and  the  Canadian  Acoustical  Association. 
The  Canadian  Acoustical  Association  awarded  him  the  Edgar  and  MiHicent 
Shaw  Postdoctoral  Prize  in  1994  and  the  Directors’  Award  in  the 
“Professional  over  30” category  in  1997. 
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Project  04-D:  Broadband  modelling  of  AS  W  monostatic 
and  multistatic  propagation  and  reverberation 


Operational  relevance 

Broadband  propagation  and  reverberation  models  are  required  to  determine  the 
performance  prediction  and  system  concept  assessments  related  to  low-frequency 
active  sonars.  Accurate  models  are  essentia!  for  significantly  improving  the 
effectiveness  of  acoustic  ASW  operations.  The  models  improve  the  predictive 
capabilities  for  both  detection  and  classification  systems,  thus  providing  an  essential 
component  of  predictive  sonar  performance  models  for  ASW.  Computer  models  also 
increase  understanding  of  how  various  physical  features  of  the  ocean  and  the  seabed 
affect  the  propagation  and  scattering  of  sonar  signals. 


Broadband  models  and  shallow  water  acoustic  variability 

The  PROSIM  (PROpagation  channel  SIMulator)  project  finished  in  April  1999.  The  project 
was  carried  out  with  partial  funding  from  the  European  Union  as  a  collaborative  effort 
within  the  Marine  Science  and  Technology  initiative,  MAST-lll.  The  other  partners  were 
Thomson-Marconi  Sonar  in  Nice,TNO-FEL  in  The  Hague,  Heriot-Watt  University  in  Edinburgh 
and  University  of  Wales  in  Bangor. 

The  scope  of  the  project  was  to  develop  a  software  package  for  signal  modelling  in  shallow- 
water  areas  with  both  spatial  and  temporal  variability.  The  completed  package  consists  of 
12  numerical  models:  3  oceanographic  models,  5  acoustic  models  and  4  auxiliary  tools 
sharing  a  user-friendly  graphical  interface  for  input  and  output  from  the  various  models. 
The  graphical  user  interface  (GUI)  and  six  of  the  oceanographic  and  acoustic  models  were 
newly  developed  during  the  project.  A  user’s  guide  describing  the  background  of  the  models 
and  how  to  create  data  inputs  and  run  the  models  was  also  provided. 

Two  major  experiments  were  conducted  during  the  project  resulting  in  extensive 
environmental  and  acoustic  data  from  shallow-water  regions  in  the  Clyde  Sea  and  the 
Mediterranean. The  objective  of  the  experiments  was  to  assess  the  degree  of  time-variability 
of  acoustic  signals  received  over  fixed  propagation  paths  and  to  correlate  the  acoustic 
variability  with  changes  in  the  environment.  The  experimental  data  have  been  used  to 
validate  the  numerical  models  in  the  PROSIM  software  package. 

The  primary  contributions  from  Project  04-D  to  the  PROSIM  development  were: 

•  a  computationally-efficient,  broadband,  range-dependent,  deterministic  propagation 
model  based  on  normal  modes 

•  a  model  for  converting  geoacoustic  bottom  properties  to  a  plane-wave  reflection 
coefficient 

•  oceanographic  and  acoustic  data  from  the  Mediterranean  Sea  experiment. 

The  work  in  1 999  concentrated  on  finalizing  the  data  analysis  and  model  validation,  preparing 
written  documentation  for  the  final  report  to  the  European  Union  and  constructing  a  web 


—  87  — 


site  hosted  by  SACLANTCEN  advertising  the  final  software  product  (http:// 
www.saclantc.nato.int/mast/prosim/mainpg/).  This  web  site  contains  a  detailed  description 
of  the  PROSIM  project,  including  all  project-related  publications  in  electronic  form.  Moreover, 
visitors  may  download  executables  of  the  developed  channel  simulator  for  a  variety  of 
computer  platforms. 


In  April/May  1999  SACLANTCEN  (Project  01 -B)  conducted  another  shallow-water 
experiment  (ADVENT’99),  which  provided  additional  data  on  time-variability  effects  on 
acoustic  signal  transmission  over  a  fixed  path. The  experiment  took  place  on  the  Adventure 


Figure  04-D.1  Measured 
sound-speed  profiles 
(towed  CTD-chain)  along 
a  10  km  track.  Some 
spatial  variability  is  evident 
and  most  features  were 
seen  to  change  both 
position  and  shape  with 
time . 


Sound  speed  profile  (m/s) 


Range  (km) 


Bank  off  the  west  coast  of  Sicily,  where  an  acoustic  track  with  a  weakly  range-dependent 
bathymetry  was  chosen.  Three  ranges  of  2,  5  and  10  km  were  considered  for  fixed-path 
propagation,  where  broadband  signals  (200-1600  Hz)  were  transmitted  every  minute  for 
up  to  18  h.  Extensive  oceanographic  data  (current,  wave  height,  sound-speed  profiles) 
were  acquired  during  the  acoustic  transmissions  in  order  to  correlate  the  oceanographic 
and  acoustic  variability.  A  conductivity-temperature-depth  (CTD)  chain  was  towed  along 
the  propagation  paths  resulting  in  a  set  of  fully  range-dependent  sound-speed  profiles 
approximately  every  hour.  An  example  of  the  sound-speed  variation  along  the  10  km 
propagation  path  is  shown  in  Fig.  04-D.1. 

Data  Model 


Frequency  (Hz)  Frequency  (Hz) 


Figure  04-D.2  Correlation  between  the  first  and  the  subsequent  received  signals  across  the  vertical 
array  as  a  function  of  acoustic  frequency  and  transmission  time  for  a  10-km  fixed  path  propagation. 
Note  in  the  left  panel  that  signals  de-correlate  faster  at  higher  frequencies,  with  the  de-correlation 
time  being  only  around  1  h  for  frequencies  above  650  Hz.  Similar  effects  are  seen  in  the  modelled 
result  in  the  right  panel,  even  though  the  overall  correlation  at  lower  frequencies  is  much  higher  here 
than  in  the  (noisy)  data. 
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A  measure  of  the  time  stability  of  the  received  acoustic  signals  is  established  by  correlating 
the  first  received  signal  with  the  subsequent  signals  as  transmission  time  progresses.  The 
correlator  is  the  normalized  linear  Bartlett  processor  which  adds  the  signals  coherently 
across  the  vertical  array  for  individual  frequencies. The  resulting  measure  of  signal  similarity 
has  a  value  between  0  and  1 ,  with  0  meaning  no  similarity  of  signals  received  across  the 
array  between  two  transmissions  and  1  meaning  two  identical  signals. 

The  correlation  of  the  experimental  data  (Fig.  04-D.2,  left  panel)  varies  considerably  with 
time,  indicating  changes  in  the  propagation  conditions  as  well  as  in  the  background  noise 
level.  This  decorrelation  of  received  signals  is  stronger  at  higher  frequencies.  The  low 
correlation  (blue  horizontal  stripes)  occurring  every  2  h  across  the  entire  frequency  band  is 
caused  by  noise  from  the  ship  towing  the  CTD-chain  when  it  approaches  the  vertical  array. 

The  effect  of  a  time-varying  water  column  on  signal  transmission  was  simulated  by  the 
SACLANTCEN  propagation  model  C-SNAP.The  measured  range-dependent  sound-speed 
profiles  over  an  18  h  time  period  were  included  in  the  calculations.  The  correlation  of  the 
simulated  signals  (Fig.  04-D.2,  right  panel)  is  high  at  the  lower  frequencies  for  the  entire 
time  period.  A  fast  decorrelation  of  the  signals  is  seen  for  the  higher  frequencies  as  observed 
in  the  experimental  data,  indicating  substantially  correct  modelling  of  the  influence  of  the 
varying  sound-speed  structure  on  sound  propagation  at  these  frequencies.  The  difference 
between  model  and  data  at  the  lower  frequencies  may  be  due  to  the  lack  of  noise  in  the 
modelled  signals. 


Propagation  model  validation 

The  broadband  range-dependent  normal-mode  model  PROSIM  has  been  applied  to  a  set 
of  test  cases  developed  for  the  Shallow-Water  Acoustic  Modelling  (SWAM)  workshop  held 
in  Monterey,  CA,  in  September  1 999.  One  of  the  cases  considered  was  simulation  of  sound 
propagation  across  a  shelf  break  from  deep  to  shallow  water  with  range-dependent 
bathymetry  and  sound-speed  profile  (Fig.  04-D.3). 

Three  different  propagation  models  were  applied  to  this  problem:  A  parabolic  equation 
model  (RAM),  a  coupled-mode  model  (C-SNAP)  and  the  recently-developed  adiabatic 
normal-mode  model  (PROSIM).  The  results  are  in  very  good  agreement  as  the  signal  is 
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Figure  04-D.3  Transmission  loss  as  a  function  of  range  (source  depth  =  30  m,  receiver  depth  =  35  m, 
frequency  =  1000  Hz)  calculated  by  using  C-SNAR  PROSIM  and  RAM.  Very  good  agreement  between 
the  models  is  seen  out  to  a  range  of  around  14.5  km  for  this  shelf-break  environment.  Then  levels 
drop  off  rapidly  and  only  RAM  provides  accurate  results  on  the  shelf  itself. : 
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propagating  up-slope  to  a  range  of  around  14.5  km  (Fig.  04-D.3,  right  panel).  Hence  the 
adiabatic  approximation  is  sufficiently  accurate  to  this  range  for  a  frequency  of  1000  Hz. 
Close  to  the  shallow-water  area,  the  predictions  from  C-SNAP  and  PROSIM  deviate 
significantly  from  RAM  (correct  result),  because  of  the  approximations  introduced  in  the 
mode  models. 

Bistatic  reverberation  model  development 

The  reverberation  modelling  in  Project  04-D  produces  models  which  can  be  used  to  support 
the  interpretation  and  understanding  of  active  sonar  datasets  acquired  under  Project  04-A. 
It  also  supports  the  scattering  mechanism  extraction  efforts  in  Project  04-C.  In  both  cases 
the  emphasis  is  on  modelling  the  reverberation  process  to  the  maximum  fidelity  possible, 
taking  advantage  of  the  state  of  the  art  in  forward  propagation  modelling  to  and  from  the 
scattering  region  and  utilizing  physically  reasonable  models  of  the  scattering  processes. 
The  approach  is  motivated  by  the  desire  to  accurately  capture  in  the  model  predictions  the 
complexity  and  richness  of  the  angular  and  temporal  evolution  of  shallow  water  reverberation. 
The  fact  that  reverberation  is  not  “White  Gaussian  Noise”  is  caused  to  a  large  degree  by 
complexity  and  spatial  variability  of  the  shallow  water  sound  channel,  as  well  as  by  the 
underlying  distribution  parameters  of  the  scatterers  themselves.  For  this  reason,  high  fidelity 
modelling  contributes  to  a  better  understanding  of  the“non-Rayleighness”  of  the  reverberation 
power,  which  causes  the  high  false  alarm  rate  of  shallow  water  active  systems. 

In  1999  two  new  reverberation  models  were  created,  1)  R-SNAP  for  modelling  long-range 
monostatic  reverberation  in  range-dependent  waveguides  and  2)  a  volume  inhomogeneity 
scattering  enhancement  to  OASES  which  is  useful  for  modelling  early-time  subbottom 


Figure  04-D.  4  Stacked 
reverberation  power  measured 
during  the  SWAC-1  experiment 
conducted  on  the  Malta  plateau. 
Unusually  strong  reverberation  is 
measured  from  the  vicinity  of  the 
Ragusa  ridge  in  the  Malta 
channel  and  from  southern  Sicily 
and  Malta 


reverberation  in  range-independent  waveguides.  The  R-SNAP  model  has  been  used  to 
support  the  interpretation  and  understanding  of  the  active  reverberation  datasets  collected 
on  the  Malta  plateau. 

In  Fig.  04-D.4  a  stacking  of  reverberation  power  collected  during  the  SWAC  1  experiment 
over  the  Malta  plateau  is  superimposed  on  the  bathymetry.  Strong  scattering  is  seen  from 
the  vicinity  of  the  Ragusa  ridge  in  the  eastern  part  of  the  Malta  channel  as  well  as  from  the 
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Figure  04-D.5  Model  predictions 
of  monostatic  reverberation 
versus,  azimuth  for  rough 
sediment  scattering  on  the  Malta 
Plateau.  The  predictions  do  not 
show  the  strong  variability  with 
azimuth  visible  in  the  dataset. 


Figure  04-D.6  Reverberation 
predictions  for  the  Malta  plateau 
with  both  rough  sediment  and 
rough  basement  scattering.  The 
predictions  show  stronger 
reverberation  localized  on  the 
Ragusa  ridge  and  off  the 
southern  shore  of  Sicily,  in  better 
agreement  with  the  experimental 
measurements.  Boomer  surveys 
indicate  that  the  sediment 
thickness  goes  to  zero  in  these 
regions,  enabling  the  rougher 
rock  basement  to  be  directly 
illuminated  by  the  incident  field. 
The  boomer  survey  did  not 
extend  to  the  northern 
escarpment  of  Malta,  thus  the 
prediction  is  less  accurate  in  this 
region. 


southern  coast  of  Sicily  and  the  northern  slope  of  Malta.  In  Figs.  04-D.5  and  04-D.6 
monostatic  reverberation  predictions  from  R-SNAP  are  shown  for  rough  sediment  (Fig.  04- 
D.5)  and  a  combination  of  rough  sediment  and  rough  basement  (Fig.  04-D.6)  scattering 
hypotheses.  The  rough  basement  predictions  were  obtained  using  sediment  thickness 
estimates  obtained  from  a  boomer  survey  of  the  Malta  plateau  conducted  by  Project  04-C. 
The  survey  showed  that  the  sediment  thickness  was  essentially  zero  over  the  Ragusa 
ridge  and  on  parts  of  the  southern  escarpment  of  Sicily,  exposing  the  underlying  rough 
basement  to  the  acoustic  energy  of  active  systems.  R-SNAP  scattering  predictions  based 
on  rough  sediment  scattering  alone  are  too  isotropic  in  azimuth,  but  when  the  basement 
scattering  mechanism  is  added,  more  features  of  the  data  are  captured.  It  is  expected  that 
the  R-SNAP  model  will  prove  useful  for  understanding  and  predicting  the  spatial-temporal 
characteristics  of  these  types  of  reverberation  data  sets  in  shallow  water,  where  significant 
spatial  diversity  in  the  reverberation  may  be  caused  by  the  inhomogeneity  of  the  underlying 
geophysical  description  of  the  bottom. 
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The  volume  scattering  OASES  module  was  developed  to  provide  an  independent  tool  for 
modelling  the  vertical  transmit  receive  array  (VETRA)  scattering  datasets  collected  in  Project 
04-C.  These  data  are  usually  represented  as  beam-time  evolutions  of  backscatter  (Fig.  04- 
D.7).  From  these  data  it  is  desired  to  resolve  the  scattering  mechanism  (rough  surface 
versus  sediment  volume  scatter)  for  the  various  layers  in  the  subbottom,  as  well  as  the  total 
scattering  strength  versus  grazing  angle  and  frequency.  The  OASES  modelling  capability 
will  provide  data  against  which  the  estimation  process  can  be  benchmarked.  Beam-time 
evolutions  of  backscatter  from  rough  surface  and  sediment  volume  inhomogeneity  scatterers 
are  illustrated  in  Figs.  04-D.8  and  04-D.9.The  differences  between  these  results  indicate 
that  these  two  mechanisms  should  be  resolvable  from  one  another  using  the  VETRA 
measurement  apparatus. 


Figure  04-D.7  Beam-time  history  of  bottom  and 
subbottom  reverberation  collected  by  the  VETRA 
apparatus  in  the  Capraia  basin  north  of  Elba.  At 
early  time  the  reverberation  arrives  on  the  Vi A  from 
steep  angles,  approaching  the  horizontal  at  late 
time.  To  the  direct  measurement  of  the  monostatic 
bottom  scattering  are  superimposed  bistatic 
multipath  which  give  additional  information  about 
the  bottom  scattering  kernel,  but  also  cause 
ambiguity  at  times  after  approximately  0.2  s. 
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Figure  04- D.  8  Model  predications  of  the  time-angle  history  of  backscatter  for  rough  surface  scattering  fora  VETRA 
experimental  geometry  The  various  trajectories  of  the  scattered  field  are  identified  in  the  accompanying  schematic. 
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Figure  04-D.9  Model  predictions  of  the  time-angle  history  of  backscatter  from  volume  inhomogeneities  in  the  sediment 
layer  below  the  VETRA  apparatus.  The  volume  scatterers  contribute  for  a  range  of  angles  at  all  times,  thus  yielding 
data  which  are  distinguishable  from  rough  surface  scattering  shown  in  Figure  04-D.8. 
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Thrust  05  Command  Support  (CDS) 

Project  05-A:  Command  and  NATO  support 


Operational  relevance 

Force  Employment  Studies  are  a  major  part  of  the  Centre’s  Operational  Research 
and  Command  Support  programme.  The  studies  make  specific  recommendations  to 
operational  commanders,  NATO  and  the  nations  on  the  optimal  use  of  undersea 
warfare  resources.  The  studies  vary  in  scope  from  advising  an  at-sea  commander  on 
the  best  use  of  sensors  to  achieve  his  mission,  to  advising  senior  NATO  and  national 
decision  makers  on  the  most  cost  effective  way  to  proceed  with  future  procurement 
and  research  and  development.  Another  key  part  of  this  work  is  operational  planning 
which  supports  NATO  crisis  response.  This  work  aids  operational  commanders  in 
determining  which  forces  they  must  commit  to  accomplish  a  particular  mission.  In  the 
past,  studies  have  typically  been  quite  specific  in  terms  of  the  scenario,  the 
geographical  area,  the  NATO  objectives  and  force  composition  and  the  threat.  These 
specific  studies  are  no  longer  appropriate  given  the  transformation  of  NATO’s  defence 
posture.  The  present  approach  is  to  produce  computer  based  decision  aids  which 
have  the  flexibility  to  cope  with  any  situation  that  could  arise.  Because  of  the  flexibility 
of  these  decision  aids,  they  can  also  be  used  to  support  higher  level  force  requirements. 


Implications  of  new  technologies  for  maritime  operations  in  2015 

The  aim  of  this  study,  initiated  by  SACLANT  and  sponsored  as  a  Long  Term  Scientific 
Study  by  the  NATO  Research  and  Technology  Board,  is  to  help  NATO  and  the  Nations  with 
long-term  defence  planning  and  new  systems  requirements,  in  the  context  of  a  transformed 
security  environment  and  diminishing  defence  budgets. 

Phase  I  of  the  study,  which  was  completed  in  1 995,  provided  a  detailed  evaluation  of  potential 
security  challenges,  ranking  critical  shortfalls  in  maritime  force  capabilities  and  assessing 
the  most  promising  and  affordable  systems  for  development  up  to  2005.  Phase  II,  which 
started  in  1997  comprises  two  Advanced  Concept  Studies  (ACS),  which  extend  the  time 
frame  to  2015,  Antisubmarine  Warfare  (led  by  France)  and  Mine  Warfare  (led  by  Germany). 

ASIV  advanced  concepts  study 

The  critical  shortfalls  in  ASW  identified  in  Phase  I,  provide  the  focus  for  the  ASW  Advanced 
Concepts  Study.  The  aim  of  this  study  is  to  recommend  the  most  promising  directions  for 
Allied  research  and  development  leading  to  affordable  system  packages  in  the  2005-2015 
time  frame,  which  reduce  the  most  significant  ASW  shortfalls.  SACLANTCEN  has  provided 
support  to  MO2015  activities  in  ASW  from  the  outset  and  is  a  member  of  the  ASW  ACS 
Core  Team  for  Phase  II. The  ASW  ACS  is  supported  by  Canada,  Denmark,  France,  Germany, 
Italy,  Netherlands,  Norway,  Spain,  Turkey,  United  Kingdom  and  the  United  States. 

The  first  step  of  Phase  II  was  to  review  and  quantify  the  most  significant  ASW  shortfalls 
identified  in  Phase  I  projected  to  persist  in  2005-2015,  in  order  to  assess  their  magnitude, 
ascertain  their  principal  causes  and  identify  areas  for  R&D.  Step  1  was  completed  in  April 
1 998  with  a  report. 
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Step  2  of  Phase  II  which  started  in  April  1998  was  concerned  with  formulating  solutions  to 
the  most  significant  ASW  shortfalls  from  Step  1 .  The  solutions  generated  by  means  of 
workshops  and  brainstorming  sessions  were  either  evolutionary  in  nature  or,  in  some  cases, 
revolutionary.  Step  2  was  completed  with  a  Mid  Course  Review  (MCR)  in  April  1999,  at 
which  the  concepts  were  reviewed  and  a  sub-set  selected  for  more  detailed  evaluation  in 
Step  3.  Step  3  culminated  in  a  Concept  Development  Workshop  (CDW)  in  November  1999, 
where  a  final  selection  of  concepts  took  place  for  refinement  and  evaluation.  Step  4  is 
concerned  with  high-level  operational  assessment  of  the  concepts  in  preparation  for  the 
Multi-National  Exercise  (MNE)  in  January  2000.The  MNE  will  comprise  two  panels,  a  Military 
Worth  Panel  and  an  R&D  Panel,  which  will  provide  the  final  evaluation  of  the  concepts  and 
decide  on  an  R&D  road  map. 

SACLANTCEN  contributions  to  the  ASW  ACS  in  1999  include  the  following: 

•  Leadership  of  the  Acoustic/Distributed  Deployable  SystemsTechnical  Provider  Group 

•  Organization  of  the  Distributed  Deployable  Systems  (DDS)  Workshop,  held  at 
SACLANTCEN  in  February  1999 

•  Writing  of  concept  data  sheets  for  the  MCR  and  CDW 

•  Performance  evaluation  of  selected  concepts  for  the  CDW  and  MNE. 

MW  advanced  concepts  study 

The  leaders  of  SACLANTCEN  projects: 

•  03-A  Buried  mine  detection  and  classification 

•  03-B  Large  area  remote  search  in  coastal  waters 

•  03-C  Mine  ship  interaction 

•  provided  baseline  data  sheets  on  UUV-based  minehunting  and  mine  jamming 
concepts  which  were  reviewed  at  the  evaluation  meeting  held  in  September  1999. 

NATO  defence  requirements  review 

At  the  beginning  of  the  year  the  Centre  supported  both  NATO’s  Strategic  Commanders 
(SCs)  in  their  assessment  of  force  requirements.  Anti-submarine  warfare  studies,  started 
in  1998,  were  completed  in  1999  and  transmitted  to  the  staff  of  the  SCs.  The  studies 
examined  the  capability  of  different  force  levels  to  conduct  area  searches,  screening  and 
barrier  operations. 

Exercise  support 


Operational  relevance 

The  overall  aims  of  the  Centre’s  participation  in  NATO  exercises  are: 

•  To  facilitate  the  flow  of  information  and  expertise  between  the  Centre  and  the 
operational  community  and  Commands  including  SACLANT,  SACEUR  and 
their  subordinate  commands. 

i 

•  To  provide  expert  assistance  to  the  NATO  maritime  analysis  agencies  (the 
Permanent  Analysis  Team  at  Northwood,  UK  and  the  Independent  Maritime 
Analysis  Team  at  Naples,  IT). 

\ 

•  To  present  Centre  scientists  with  the  opportunity  to  test  new  concepts  in  an 
operational  situation  and  for  naval  forces  to  receive  experience  and  training 
in  the  use  of  these  new  techniques . 
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Exercise  Support  comprises  participation  in  the  planning,  analysis  (and  sometimes  conduct) 
of  exercises  and  in  the  provision  of  specialized  exercise  support  tools. 

Exercise  planning  and  analysis 

DOGFISH 

ORD  made  a  significant  contribution  to  the  analysis1  by  carrying  out  all  the  necessary 
reconstruction.  An  existing  EXCEL  based  system,  originally  developed  to  facilitate  the  in- 
house  analysis  of  low  frequency  active  sonar  (LFAS)  data,  was  modified  to  provide  more 
extensive  reconstruction  and  analysis  capabilities.  The  success  of  the  system,  which  was 
used  during  the  analysis  conference  resulted  in  a  formal  COMNAVSOUTH  proposal  for  the 
SACLANTCEN  Scientific  Programme  of  Work  (SPOW)  2001  to  further  develop  the  system. 

BLUE  GAME 

ORD  provided  support  for  the  planning  of  the  exercise  and  fitted  participating  fast  patrol 
boats  with  the  Electronic  Minefield  Referee  (EMIR)  system  (see  below).  Members  of  ORD 
attended  the  analysis  conference  and  carried  out  a  detailed  quantitative  analysis  of  the 
performance  of  a  number  of  operational  NATO  minehunting  systems  within  a  “controlled” 
environment2  .The  results,  which  have  been  presented  at  a  number  of  NATO  MW  meetings, 
generated  considerable  interest. 

BELL  BOTTOMS 

ORD  attended  the  synthetic  computer  war  game  and  contributed  to  the  exercise  analysis3 . 

Exercise  support  tools 

Maritime  Electronic  Log  (MEL) 

MEL  was  developed  for  Exercise  DOGFISH  1999.  Its  purpose  is  to  replace  the  paper 
versions  of  FORMEX’s  110  and  114  with  an  accurate  electronic  log.  It  was  successfully 
employed  during  DOGFISH  on  board  8  submarines  and  1  surface  ship.  MEL  features  a 
GIS  interface  and  an  interface  with  commercial  NMEA  0183  compliant  GPS  systems. 

Electronic  Minefield  Referee  (EMIR) 

EMIR  simulates  the  interaction  between  a  synthetic  minefield  and  surface  and  sub-surface 
exercise  participants.  During  the  exercise  planning  phase  ‘paper’  minefields  are  coded  into 
the  EMIR  software  using  the  dedicated  Planner  software.  During  the  exercise  a  laptop 
computer  is  installed  on  participating  units  which  contains  the  EMIR  software  with  the 
embedded  coded  minefields. Throughout  the  exercise  the  EMIR  software  checks  the  position 
of  the  participating  unit,  compares  it  to  the  position  of  the  minefields  and  provides  a  warning 
if  a  mine  is  deemed  to  have  detonated  and  damaged  the  unit.  Originally  units  had  to  enter 
their  positions  throughout  the  exercise.  Now,  for  surface  units,  it  is  possible  to  connect  a 
handheld  GPS  unit  which  automatically  feeds  positional  data  to  EMIR.  In  February,  the 


1  NATO  Southern  Region  Independent  Maritime  AnalysisTeam  (IMAT).  DOGFISH  99  exercise  analysis  report. 
NATO  COMNAVSOUTH  1 999  (NATO  CONFIDENTIAL).  [COS/AA/99], 

2  NATO  Maritime  Central  AnalysisTeam  (MCAT).  Exercise  BLUE  GAME  99  analysis  report.  NATO 
CINCGERFLEET,  1999.  (NATO  CONFIDENTIAL). 

3  NATO  COMNAVNORTHWEST.  Exercise  BELL  BOTTOMS  99  analysis  report. 
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MEL  main  screen  with  GIS  display  window 


JMC  961  exercise  was  supported  and  EMIR  installed  on  board  2  UK  Navy  units.  This 
exercise  provided  useful  feedback,  which  led  to  several  improvements  in  the  software.  In 
April,  during  exercise  BLUE  GAME  1999,  EMIR  was  installed  on  board  20  Fast  Patrol 
Boats  (FPBs)  and  proved  to  be  effective  in  enforcing  monitoring  of  the  mine  threat.  BLUE 
GAME  2000  will  be  supported.  A  new  improved  version  of  the  planner  program  has  been 
developed,  which  will  include  a  track  analysis  module. 

The  UK  has  established  an  EMIR  User  Group  the  first  meeting  of  which  was  held  in 
Northwood,  UK  during  November  1999. 


EMIR  planner  main  screen 
and  GPS  configuration 
window 
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GPS  PAD 
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Tactical  Decision  Aid  for  Planning  ASW  Barrier  Operations  in  the  Strait  of 
Gibraltar  (STROG  TDA) 

The  objective  of  this  task  is  to  develop  a  prototype  Tactical  Decision  Aid  which  permits  the 
optimal  allocation  of  available  units  to  antisubmarine  barriers  in  the  Strait  of  Gibraltar.  The 
project  was  initiated  in  January  1999.  A  document  containing  the  user  requirements  was 
submitted  to  the  Spanish  Navy  and  approved.  A  prototype  of  the  man-machine  interface 
was  designed  and  delivered  to  the  Spanish  Navy. 

A  genetic  algorithm  has  been  developed  to  optimally  allocate  available  units  to  the  barrier. 
The  algorithm  needs  to  be  provided  with  the  performance  of  every  unit  on  every  cell  in  the 
area. The  solution  can  be  a  connected  or  disconnected  barrier. The  time  to  obtain  a  solution 
for  an  area  containing  500  cells  and  5  units  is  considered  to  meet  user  requirements.  The 
algorithm  designs  the  barrier  taking  into  account  the  coastline  and  a  preferred  transit  axis 
for  the  submarines. 

Work  continues  to  make  the  algorithm  more  general  and  applicable  to  any  geographic 
area.  The  way  ahead  is  to  integrate  the  developed  solver  with  the  previously  delivered  user 
interface,  to  document  the  work  done  and  deliver  the  end-product  by  June  99.  The  format 
for  the  data,  which  defines  the  performance  of  each  unit  in  the  area,  will  be  clearly  defined. 

Mine  Countermeasures  Exclusive  Planning  and  Evaluation  Tool  (MCM 
EXPERT) 

MCM  EXPERT  is  an  experimental,  PC  based  software  tool  for  the  planning  and  evaluation 
of  NATO  MCM  operations.  MCM  EXPERT  was  developed  by  the  NC3A  and  incorporates 
algorithms  developed  by  a  NATO  Ad  Hoc  Working  Group  consisting  of  representatives 
from  NATO  agencies,  NATO  Commands  and  national  research  laboratories.  As  part  of  this 
working  group,  SACLANTCEN  is  responsible  for  the  standard  NATO  method  to  calculate 
the  threat  posed  by  mines  to  a  defined  target  vessel.  The  software  is  now  maintained  by  the 
NATO  Integrated  Software  Support  Centre  (ISSC).  During  1999,  SACLANTCEN  continued 
to  support  the  MCM  EXPERT  User  Group  (MEUG)  meetings. 
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MCM  Expert  meeting  in  the  SACLANTCEN  main  conference  room 


Decision  Aid  for  Risk  Evaluation  (DARE) 

During  1998  the  Centre  sponsored  its  first  NATO  Experimental  Tactic  (EXTAC)  -  DARE 
(Decision  Aid  for  Risk  Evaluation)  which  assists  Maritime  Mine  Countermeasures  (MCM) 
Commanders  to  assess  the  risk  to  follow-on  traffic/naval  operations  arising  from  uncountered 
mines  after  the  completion  of  a  MCM  operation. The  software  has  continued  to  be  supported 
by  resolving  queries  arising  from  its  use  during  operations  and  exercises.  The  tactic  has 
now  been  issued  to  1 1  Nations  and  Commands. 

Search  and  Screening  Working  Group 

This  is  a  working  group  under  the  auspices  of  the  Maritime  Tactical  Working  Group.  The 
Centre  has  been  a  regular  participant  in  the  proceedings  of  the  group  in  preparation  for 
2000,  when  the  Centre  will  assume  the  responsibility  for  developing  an  area  search  tactical 
planning  aid. 


Ken  Pye  holds  a  B.Sc  from  Sheffield  University  England  and  an  M.Sc  in 
Operational  research  from  Strathclyde  University,  Scotland.  His  career  in 
operational  research  spans  some  30  years.  In  the  1970's  he  worked  on 
weapon  systems  development  and  conducted  studies  at  the  then  SHAPE 
Technical  Centre,  the  Netherlands.  In  1981  he  joined  the  NATO  E-3A  Com¬ 
ponent  at  Geilenkirchen,  as  Chief  of  the  Plans  and  Analysis  Office  in  the 
Software  Support  Centre.  In  the  1990’s  he  was  the  Operational  Research 
Scientist  Scientific  Advisor  to  the  Former  CINCHAN  and  then  to 
CINCEASTLANT/COMNAVNORTHWEST.  He  joined  the  Centre  in  1998  as 
Head  of  the  Operations  Research  Department. 


Project  05-B  -  Planning  support  for  underwater  warfare 
operations 

The  objective  of  this  project  is  to  provide  decision  support  aids  to  NATO  Commands  and 
Nations  for  planning  underwater  warfare  operations  for  future  crisis  operations,  the  location 
and  timing  of  which  are  unpredictable.  Activities  have  concentrated  on  providing  opera¬ 
tional  planning  tools  for  ASW  operations.  Mediterranean  Operational  Planning  Software 
(MEDOPS)  evolved  from  WESTMED,  which  was  limited  to  the  western  Mediterranean  Sea. 
Both  projects  were  sponsored  by  Commander  Naval  Forces  Southern  Europe 
(COMNAVSOUTH). 

PLANET  (Planning  Expert  Tool),  sponsored  by  SACLANT,  will  increase  the  geographical 
coverage  of  the  ASW  planning  tool.  Figure  1  shows  the  relationship  and  progression  of 
these  planning  tools. 
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Evolution  of  operational  and  defence  planning  decision  aids  at  SACLANTCEN. 


MEDOPS  -  Mediterranean  Operational  Planning  Tool 

MEDOPS  is  a  decision  aid  which  combines  high  quality  sonar  modelling  with  sophisticated 
search  algorithms.  Mission  planners  can  quickly  determine  the  most  efficient  force  deploy¬ 
ment  with  software  running  on  a  standard  PC.  Intuitive  interfaces  and  extensive  on-line 
help  allow  the  configuration  of  complex  operations  without  specialist  training  or  computer 
knowledge. 

MEDOPS  provides  easy  access  to  generic  sets  of  ASW  platforms,  reflecting  NATO  ASW 
capability.  These  assets  can  be  assigned  to  any  mission  and  stored  in  a  database  for 
subsequent  use. 

High-resolution  vector  maps  are  provided  to  allow  the  user  to  select  the  proper  mission 
area.  Purpose  -  built,  high  -  speed  zooming  tools  ensure  that  accurate  scenario  definition  is 
possible. 

Results  are  provided  in  report  and  graphical  format.  Acoustic  sensor  performance  predic¬ 
tions  at  high  and  low  frequencies  cover  12  months. 

Fast,  high  order  interpolation  and  graphic  algorithms  have  been  developed  to  ensure  that 
the  maps  have  a  high  degree  of  smoothness. 

A  Steering  Group,  consisting  of  Command  and  National  representatives,  has  been  estab¬ 
lished  to  guide  this  project  and  the  follow-on  work  on  PLANET 
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The  major  accomplishments  of  1999  are: 

•  MEDOPS  evaluation  version  released  in  January  1999. 

•  Steering  group  meeting  in  March  1 999.  (Principle  activity:  review  of  MEDOPS  evaluation 
version) 

•  MEDOPS  Version  1 .0  released  in  June  1999  —  based  on  Steering  Group  review. 

•  Steering  group  meeting  in  September  1999.  (Review  of  MEDOPS  Version  1.0) 

•  MEDOPS  Version  1.1  released  in  October  1999  (Included  further  software 
enhancements  based  on  September  Steering  group  review  which  had  not  been 
identified  in  the  March  meeting). 

•  MEDOPS  Version  1.1  at-sea  validation  exercise  conducted  in  conjunction  with  UKN. 
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Thrust  06  Exploratory  Research  (EXR) 

Project  06-B:  Focused  acoustic  field  (FAF)  studies 
(Alliance  days  -  23,  Manning  days  -  24) 


Operational  relevance 

Although  this  work  is  at  the  early  concept  demonstration  stage,  it  is  already  apparent  that 
the  operational  applications  of  FAF  include  multistatic  active  sonar,  low  frequency  active 
sonar  (LFAS),  acoustic  communication  systems  with  submarines,  air  deployable/fixed 
surveillance  systems  and  the  detection  of  low  target  strength  objects  and  buried  mines. 


A  high  frequency  (3.5  kHz)  phase  conjugation  experiment  was  conducted  jointly  by 
SACLANTCEN  and  the  Marine  Physical  Laboratory  (MPL)  of  the  Scripps  Institution  of 
Oceanography,  University  of  San  Diego,  in  July.  A  new,  high  frequency  vertical  array  (SRA) 
of  29  source/receive  transducers  operating  nominally  in  the  3-4  kHz  band,  an  underwater 
pressure  case  containing  the  source/receive  electronics,  a  surface  buoy  providing  battery 
power,  system  control  and  wireless  local  area  network  (LAN)  connectivity,  were  designed 
and  fabricated  by  MPL. 


1999  Sea  trial 

This  third  experiment  (FAF  99)  of  the  Joint  Research  Project  was  designed  to  exploit  the 
concept  of  phase  conjugation  (PC)  to  acoustically  focus  in  the  water  column.  The  signal 
from  the  probe  source  is  received  by  the  SRA  where  it  is  time  reversed  (the  time  domain 
equivalent  of  phase  conjugation)  and  retransmitted.  The  focused  field  is  measured  by  the 
vertical  receive  array  (VRA)  which  is  contiguous  (there  may  be  an  out  of  plane  displacement) 
with  the  probe  source.  The  first  two  experiments  used  50  m,  450  Hz  pulses.  Sharp  focus 
was  observed  at  ranges  of  up  to  30  km.  Methods  were  developed  and  verified  to  move  the 
focus  in  range  away  from  the  PS  position.  The  PC  process  was  shown  to  be  robust  and 
stable.  The  experiment  was  carried  out  adjacent  to  Formiche  di  Grosseto  Island  and  north 
of  the  Island  of  Elba  (Fig.  06-B.1).The  former  provided  a  link  to  experiments  conducted  in 
1996  and  1997  while  the  latter  provided  a  brief  opportunity  to  explore  a  new  environment 
with  different  conditions  resulting  in  four  major  accomplishments: 

•  High  frequency  (3.5  kHz)  phase  conjugation  focusing  at  ranges  up  to  14  km  in  flat 
and  sloping  coastal  regions. 

•  Acquisition  of  data  to  investigate  the  short-term  temporal  and  spatial  properties  of 
the  focal  region  in  two  different  bottom  regimes. 

•  Demonstration  of  phase-conjugation  processing  in  acoustic  communications. 

•  Demonstration  of  source/receive  array  technology  operating  as  a  node  on  a  wireless 
local  area  network  (LAN). 
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44=00'  |  A  schematic  of  FAF-99  in  the 

vicinity  of  Formiche  showing 
data  (10  msec  pulse  centred  at 
3.5  kHz)  is  given  in  Fig.  06-B.2. 
A  probe  source  (PS)  co-located 
with  the  vertical  receive  array 
(VRA)  ensonifies  the 
waveguide.The  resulting  signal 
received  on  the  source/receive 
array  (SRA)  is  time-reversed 
and  retransmitted  by  the  same 
transducers.  The  degree  to 
which  this  retransmitted  energy 
refocuses  at  the  source  (as 
observed  by  the  VRA)  is  used 
as  a  measure  of  phase 
conjugation  processing. 

Range  and  time 
dependence: 

Additional  examples  of  time 
reversal  focusing  as  a  function 
of  range  are  given  in  the 
composite  illustration  in  Fig.  06- 
B.3.  The  variation  in  vertical 
10=30'  „w  no20.  focusing  is  due  to  slight  offsets 

between  the  PS  and  VRA  at  the 
four  ranges.  Similar  results  are 
available  as  a  function  of  PS 
depth  in  the  waveguide.  In  addition  to  exploring  the  range  and  depth  characteristics  of  the 
focal  region,  we  also  investigated  the  temporal  stability.  Although  clearly  sensitive  to 
environmental  fluctuations,  focusing  was  observed  to  be  stable  at  least  in  the  short-term. 
Figure  06-B.4  shows  a  sequence  of  10  msec  pulses  observed  by  the  VRA  from  SRA 
retransmissions  of  the  same  pulse  over  a  10  min  period  at  a  range  of  13  km. 
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Figure  06-B.1  Experimental  area. 


Figure  06- B. 2  The  FAF-99  experiment 
in  the  vicinity  of  Formiche  illustrating  the 
phase  conjugation  process  with  actual 
data.  A  probe  source  (PS)  co-located 
with  the  vertical  receive  array  (VRA) 
ensonifies  the  waveguide.  The  resulting 
signal  received  on  the  source/receive 
array  (SRA)  is  time-reversed  and 
retransmitted  by  the  same  transducers. 
The  retransmitted  energy  refocuses  at 
the  source  as  observed  by  the  VRA . 


Range  (kin) 


Figure  06- B. 3  Time  reversal 
focusing  as  a  function  of  range. 
The  variation  in  vertical  focusing 
is  due  to  slight  offsets  between 
the  PS  and  VRA  at  the  four 
ranges. 


Figure  06- B. 4  Sequence  of  10  m  pulses  observed  by  the  VRA  from  SR  A  retransmissions  of  the  same  pulse  over 
1 0  min  period  at  a  range  of  13  km. 


Working  in  the  vicinity  of  North  Elba  provided  an  opportunity  to  investigate  focusing  in  a 
different  bottom  environment  and  to  acquire  range-dependent  data.  A  particularly  interesting 
result  was  our  ability  to  demonstrate  focusing  over  a  10  km  propagation  path  from  112  m 
deep  water  (SRA)  into  32  m  deep  water  (PS  and  VRA)  as  shown  in  Fig.  06-B.5.  The  focus 
has  a  depth  extent  of  ~2m. 


Frequency-shifting  for  refocusing 

In  earlier  450  Hz  FAF  experiments,  frequency-shifting  algorithms  were  developed  to  move 
the  focus.  It  has  been  demonstrated  that  these  algorithms  can  be  used  at  3.5  kHz.  Figure 
06-B.6  below  shows  a  pulse  being  “refocused”  400  m  from  its  natural  focus  and  brought 
back  into  focus.  The  experiment  is  accomplished  by  frequency  shifting  the  same  received 
SRA  pulse  as  the  NR/V  Alliance  (with  the  SRA)  drifts. 


Pulse  sequences 

Data  were  acquired  to  study  applications  of  phase  conjugation  to  underwater 
communications.  This  first  coded  sequence  was  originally  transmitted  as  a  total  sequence 
from  the  probe  source  and  then  time  reversed  and  retransmitted  from  the  SRA.  In  this 
sequence,  a  2  msec  probe  signal  is  transmitted  from  the  probe  source.  The  time  reversed 
signal  at  the  SRA  is  a  result  of  convolving  a  50  bit  random  sequence  with  the  time  reversed 
probe  source  signal,  resulting  in  a  100  msec  sequence.  In  Fig.  06-B.7,  the  two-way  time 
reversed  sequence  (on  the  right)  is  compared  with  a  one-way  non-time-reversed  control 
sequence  (on  the  left). 
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Figure  06-B.5  Upslope  focus  near  Elba.  The  SRA  was  at  a 
depth  of  112  m  and  the  PS  and  VRA  (MFA)  were  co-located 
at  32  m  depth.  The  focus  has  a  depth  extent  of -2  m. 


Figure  06- B. 6  Frequency  shifting  for  refocusing. 


Figure  06-B.7  A  50  bit  random 
sequence  of  2  m  pulses.  Left  column: 
One  way.  Right  column:  two-way  time 
reversed. 
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Figure  06-B.8  Multipath  spread  ofa2mPS  pulse  (source  depth  40  m)  observed  at  the  SRA  over  a  relatively  range- 
independent,  7.2  km  propagation  path  at  a  depth  oft  110  m  near  Elba  (left).  Pulse  received  at  the  VRA  after  time 
reversal  and  retransmission  ( right). 

The  multipath  spread  for  a  2  m  pulse  transmitted  over  a  relatively  range-independent,  7.2  km 
propagation  path  in  1 10  m  deep  water  north  of  Elba  was  of  the  order  of  30  m  (Fig.  06-.B.8). 
After  time  reversal  and  retransmission,  the  pulse  was  compressed  back  to  ~2  m  duration. 
Packets  of  50  bits  were  then  generated  by  convolving  a  random  sequence  of  {+1  ,-1  }’s  with 
the  time-reversed  reception  at  the  SRA  of  a  2  m  probe  source  pulse.  The  time  series  of 
these  packets  received  at  the  VRA  is  being  compared  with  the  received  time  series  of 
packets  generated  at  the  SRA  by  convolving  the  random  sequence  with  a  2  m  pulse. 


2D  Tranmtssion  Loss  vs.  Range  &  Depth  SD=80m 


Modelling 

During  1999,  acoustic  propagation 
models  were  modified  to  include 
back  propagation  for  1999 
experimental  data  interpretation. 
Figure  06-B.9  shows  a  comparison 
between  data  and  model  for 
transmission  loss.  Figure  06-B.10 
shows  the  simulated  focal  area  for 
an  80  m  deep  source  at  5  km 
distance.  The  results  confirm  that 
the  focal  structure  observed  is 
consistent  with  theory. 


Figure  06 

depth. 


B.9  TL  data-  model  comparison  for  80  m  source 


Figure  06- B.  10  Single  frequency  time  reversed 
spatial  focus  versus  range  and  depth.  The  source 
depth  was  80  m  and  range  5  km. 
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Project  06-C:  Sound,  oceanography  and  living  marine 
resources  (SOLMR) 


Operational  relevance 

To  formulate  a  Marine  Mammal  (Acoustic)  Risk  Mitigation  Policy  designed  to  limit  the 
exposure  of  human  divers  and  marine  mammals  to  low-frequency  sound. 


Background 

The  following  procedures  have  been  formally  instigated  during  trials: 

•  Visual  and  passive  acoustic  (i.e.,  listening)  monitoring  to  ensure  that  marine  mammals 
are  distant  from  an  acoustic  source 

•  Slow  increase  of  acoustic  power  to  allow  marine  mammals  sufficient  opportunity  to 
leave  the  test  area  in  the  event  that  visual  and  passive  searches  are  unsuccessful 

•  Creation  of  a  database  of  marine  mammal  sightings  and  strandings  to  identify  the 
migratory  behavior  of  species  so  that  acoustic  testing  can  be  conducted  in  areas 
where  marine  mammals  are  unlikely  to  be  encountered 

•  Continuous  review  of  the  Marine  Mammal  Risk  Mitigation  Policy  to  incorporate  new 
techniques  and  information,  e.g.,  do  not  plan  experiments  near  known  diving  areas 
or  marine  mammal  breeding  grounds.  The  policy  requires  that  marine  mammals 
should  not  be  exposed  to  sound  levels  exceeding  160  dB  re  1  pPa,  which  in  practice 
amounts  to  maintaining  a  “mammal-free  zone”  around  the  source 

Advanced  passive  acoustic  sensors,  developed  originally  for  ASW  are  being  studied  as  a 
means  of  detection  and  classification  of  underwater  “biological”  sound,  in  order  to  reliably 
establish  and  monitor  that  an  experimental  area  is  free  of  marine  mammals. 

At-  sea  measurements  using  dual  use  technology 

Progress  in  basin  wide  oceanographic  modelling  and  remote  sensing  and  the  ability  to 
acoustically  monitor  the  environment  for  prolonged  periods,  allows  a  multidisciplinary, 
multinational  approach  to  predicting  the  population  density  of  cetaceans  on  a  seasonal 
basis.  A  multi-year  at-sea  effort  was  planned,  to  evaluate  risk  mitigation  techniques,  which 
use  existing  monitoring  technologies,  acquire  associated  environmental  information  and  to 
determine  the  correlation  of  these  parameters. 

Dual  Use  Technology  developed  at  the  Centre  and  satellite  remote  sensed  data  were  used 
during  the  first  sea  trial,  Sirena  ’99,  in  the  Ligurian-Corsica  basin,  which  was  recently 
designated  as  the  first  International  Marine  Sanctuary  in  the  Mediterranean  Sea. 

The  US  Office  of  Naval  Research  (ONR)  participated  by  providing  funds  to  map  the  lower 
trophic  level  distribution  and  to  tag  cetaceans  in  order  to  increase  knowledge  of  whale 
behaviour. 
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The  Italian  Navy  provided  the  research  vessel  for  the  oceanographic  and  lower  trophic 
level  measurements  and  a  Sea  King  helicopter  for  visual  surveys  prior  to  the  cruise. 

Although  it  is  inevitable  that  monitoring  and  observation  is  subject  to  significant  temporal 
and  spatial  discontinuities,  it  constitutes  a  valuable  supplement  to  existing  models. 


aT  -  -i  tv  •  JL— -1  V*  /} 4  jK  M  *>.  (  /  ■  Cuvier’s  beaked  whale  (Ziphius  cavirostris)  (Alliance) 

V  S  iv ^  Fin  Whale  (Balaenoptera  physalus)  (Alliance) 

|  ^  ®  Pi,ot  wha,e  (Globicephala  melas)  (Alliance) 

I  ®  Striped  dolphin  (Stenella  coeruloalba)  (Alliance) 

^  . ,/\  Jkf  - JR  ®  ^determined  Delphinidae  (Alliance) 

A  #  Indetermined  Cetacean  (Alliance) 

A fj  A  Fin  Whale  (Balaenoptera  physalus)  (Magnaghi) 

B  Striped  dolphin  (Stenella  coeruloalba)  (Magnaghi) 

.  .  .  ,.  ,;L  ;  ;  ®  Indetermined  Delphinidae  (Magnaghi) 

Nark  ThU itoJL  Sl?na  ???uise.  Alliance  conducted  acoustic  and  visual  measurements  along  the  track  shown  in 

ttvstracl^shown  f T'  Mf^naghi  conducted  oceanographic  and  prey  field  measurements  along 

the  track  shown  tn  green.  The  symbols  show  the  visual  sighting  of  the  cetaceans  by  both  vessels.  y 

Figure  06-C.1  summarizes  acoustic  and  visual  sightings  during  Sirena  99  (3394  nautical 
miles  of  navigation,  320  hours  of  visual  observations  and  over  340  hours  of  real  time 
acoustic  monitoring  in  an  1 1  day  period). 

Tagging  operations  from  a  small  vessel  were  conducted  prior  to  the  main  portion  of  the 
oceanographic  and  acoustic  cruise,  but  adverse  weather  conditions  prohibited  concurrent 
agging  operations.  On  the  Alliance,  a  dedicated  visual  survey  team  conducted  visual 
observations  from  the  flying  bridge, 

which  provided  a  platform  1 8m  above  , _ _ 

the  water  line.  The  acoustic  watch  i  ?  ^  f 

was  conducted  in  the  scientific  lab,  "  "  ” :  (  y:~L 

where  there  was  constant  monitoring  “  i  .. 

of  the  acoustic  display.  A  comparison  •:  Jr"  \ 

of  visual  and  passive  acoustic  ^  j  Jn  \ 

detection  methodologies, 

demonstrated  that  combining  these  J\  W  K*  $ 

techniques  is  most  effective.  Visual  L  »  "  f 

watch  is  restricted  during  high  seas  JjPf ...  \*r| J f  j 

and  low/no  light  conditions.  Passive  /f[, T-S  ' 


Figure  06-C.2  Representative  sea  surface  j  lE 
temperature  during  the  Sirena  99  cruise  period 
(5  August  1999).  Red  lines  outline  the  Ip 
boundary  of  the  Ligurian  Sea  Cetacean  I.N 
Sanctuary.  §&f:j 
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Figure  06-C.3  A  one  second,  25  Hz  - 18  Hz  low  frequency,  down  sweep  signal  from  a  fin  whale. 
Photographs  show  visual  sightings. 


acoustic  techniques  are  effective  for  detection,  localization  and  tracking  only  when  the 
animals  vocalize 

Remotely  sensed  data  (SeaWIFS  ocean  colour  and  AVHRR  sea  surface  temperature) 
were  acquired  in  real  time  at  SACLANTCEN.  Ancillary  physical  satellite  data  from  the 
Colorado  Center  for  Astrodynamics  Research  (CCAR)  were  used  to  interpret  mesocale 
physical  and  biological  oceanography  in  the  Ligurian  Sea  during  the  summer,  where  upwelling 
provides  a  nutrient  rich  cetacean  feeding  area.  (Fig.  06-C.2).The  concurrent  oceanographic 
and  lower  trophic  level  measurements  are  being  integrated  with  the  satellite  images  in  a 
geographical  information  system  (GIS).  Sirena  99  demonstrated  the  feasibility  of  the 
establishment  of  a  comprehensive  GIS-based  database  of  marine  life  in  conjunction  with 
relevant  physical/biological  parameters. 

Eight  species  of  cetaceans  are  commonly  found  in  the  Mediterranean.  These  include:  Fin 
Whale  (Balaenoptera  physalus),  Sperm  Whale  (Physeter  macrocephalus),  Cuvier’s  beaked 
whale  ( Ziphius  cavirostris),  Long  Finned  Pilot  whale  (Globicephala  melas),  Risso’s  dolphin 
(Grampus  griseus ),  Bottlenose  dolphin  (Tursiops  truncatus),  Common  dolphin  (Delphinus 
delphis)  and  Striped  dolphin  (Stenella  coeruleoalba).  Sightings  of  ten  other  species  have 
been  reported.  Cetaceans  emit  a  variety  of  sounds,  clicks,  grunts  and  whistles  with  widely 
differing  time  and  frequency  characteristics. 

During  Sirena  99,  acoustic  and/or  visual  detections  were  made  of  five  common 
Mediterranean  cetaceans,  including  fin  whales,  alone  or  in  small  groups  on  most  days  at 
ranges  of  up  to  9  km.  Figure  06-C.3  shows  a  very  low  frequency,  fin  whale  signal,  detection 
of  which  is  difficult,  due  to  high  background  noise  levels  in  this  frequency  band.  Long  finned 
pilot  whales  were  encountered  once,  when  two  pods  were  seen  in  close  proximity  to  the 
Alliance.  Figure  06-C.4  shows  wideband  whistles  recorded  from  a  single  individual. 
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Figure  06-C.4  Broadband  whistles  from  a  long  finned  pilot 
whale.  Photograph  shows  visual  sighting  of  two  adult 
juvenile 


Vocalizing  animals  can  be  localized.  Figure  06-C.5  shows  the  position  of  a  sperm  whale 
10.5  nautical  miles  from  the  Alliance,  exploiting  the  capacity  of  the  towed  array  to  distinguish 
azimuthal  direction  over  time  as  the  ship  moves.  Several  manoeuvres  were  executed  to 
cross  fix  the  animal  via  target  motion  analysis  (TMA).  As  the  whale  was  outside  the  authorized 
operating  area,  visual  identification  was  not  possible.  Figure  06-C.6  shows  acoustic 
localization  of  a  fin  whale  using  trilateration,  using  time  of  arrival  (TOA)  differences  from 
separate  sensors  to  cross  fix  the  source.  The  location  estimated  from  this  localization 
corresponds  well  to  the  estimated  range  and  bearing  of  a  visual  observation  minutes  earlier 
and  indicated  on  the  figure. 


Passive  localization  is  moderately  effective.  Marine  mammal  vocalizations  cover  a  very 
wide  frequency  range  and  vary  from  extremely  narrowband  to  extremely  wideband.  Bearing 
information  from  a  towed  line  array  is  useful,  e.g.,  as  a  cue  to  visual  observation,  however 
the  critical  parameter  for  mitigation  is  the  range.  TMA  is  effective  for  single  animals  that 
vocalize  regularly.  The  effectiveness  of  trilateration  of  time-of-arrival  information  for  any 
species,  using  simple  equipment,  is  limited  only  by  sensor  configuration  and  characteristics. 


Figure  06-C.5  Position  estimate  of  a  sperm  whale 
measured  as  a  distance  from  the  Alliance.  The 
lines  ending  in  crosses  are  the  beam  directions  at 
various  times  along  the  tracks.  The  beams  on  the 
left  side  of  the  image  converge  on  the  area 
indicated  by  the  circle.  The  beams  on  the  right  of 
the  image  diverge. 


Cetacean  database  (in  partnership  with  the  Centro  Studi 
Cetacei) 

The  Mediterranean  cetacean  database  was  conceived  to  determine  areas  of  high  and  low 
cetacean  population  density  in  the  planning  stages  of  trials  using  high  power  sound  sources. 
The  database  is  designed  to  accelerate  data  entry,  integration  and  dissemination,  based 
on  an  interactive  WEB  server.  The  server  will  incorporate  an  interactive  geographic 
visualization  system  allowing  data  to  be  searched  by  area.The  database  consists  of  several 
thousand  records  of  sighting  and  stranding  data  in  Italian  waters  and  on  Italian  coasts. 


Eleven  years  of  stranding  data  (2000  records  from  1986-1996)  have  been  digitized.  Figure 
06-C.7  shows  the  data,  sorted  by  family.  The  data  may  be  summarized  and  plotted  by 
species,  year,  season,  geographic  area  and/or  cause  of  death.  The  causes  of  these  events 
include  entanglement  in  fishing  nets,  animals  which  have  been  deliberately  killed  and  natural 


death.  More  than  4000  Italian  sighting 
records  have  been  added  to  the 
original  database. 

Sirena  sea  trials  will  continued 
systematic  monitoring  coupled  with 
environmental  measurements  to 
provide  the  understanding  of 
cetaceans  in  their  environment.  The 
database  will  be  extended  to  include 
other  regions  in  the  Mediterranean. 

For  public  outreach  and  education, 
the  Centre  has  established  working 
relationships  with  the  Acquario  di 
Genova,  the  University  of  Pavia,  the 


whale  using  triiateration. 

Istituto  Centrale  per  la  Ricerca 
Applicata  al  Mare,  (ICRAM)  World  Wide 
Federation  for  Nature  and  the  Centro 
Studi  Cetacei  to  disseminate  SOLMR 
products. 


i?  Delphinidae 
c  Ziphidae 

•  Phisiteridas 

■  BaJenopteridaa 

•  Unknown 

Figure  06-C.7  Location  of  cetaceans 
(see  key)  stranded  along  the  Italian 
coast  1986-1996. 
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Project  06-C  publications  and  presentations 

Bondaryk,  J.E.,  D’Amico,  A.,  Portunato,  N.  Techniques  for  Passive  Detection  and  Localization  of  Marine  Mammals 
using  Wide  Aperture  Arrays.  Journal  of  the  Acoustical  Society  of  America,  1 999, 106, 1 999:21 88:3aABA2. 

D’Amico,  A.,  Bondaryk,  J.E.,  Fossati,  C.,  Manghi,  M.,  Pavan,  G.,  Portunato,  N.,  Priano,  M.  Acoustic  detection  and 
localization  of  marine  mammals  in  the  Ligurian  Sea.  CSC  Annual  Conference,  Milan,  11-12  November  1 999. 

D’Amico,  A.,  Bondaryk,  J.E.,  Portunato,  N.  Passive  Detection  and  Localization  of  Marine  Mammals  in  the  Ligurian 
Sea.  Journal  of  the  Acoustical  Society  of  America,  1 06, 1 999:21 64:2pAB1 3. 

D’Amico,  A.,  Fossati,  C.,  Pavan,  G.,  Podesta,  M.,  Application  of  Graphical  Digital  Tools  to  Stranding  Information: 
Mapping  of  the  Strandings  off  the  Italian  Coast,  1986-1996,  Proceedings  of  the  European  Cetacean  Society  13,h 
Annual  Conference,  Valencia,  Spain,  April  1 999. 

D’Amico,  A.  Sirena  ’99  -  Quick  look  presentations  and  preliminary  data,  SACLANTCEN  CD-27,  October  1 999. 

Fossati,  C.,  D’Amico,  A.,  Manghi,  M.,  Pavan,  G.,  Podesta’,  M.,  Portunato,  N.,  Priano,  M.,  Teloni,  V.  Marine  mammals 
and  oceanographic  parameters: The  Geographic  Information  System  as  a  tool  for  their  organization  and  integration  in 
the  Mediterranean  Sea.  CSC  Annual  Conference,  Milan,  11-12  November  1 999. 


Angela  D’Amico  received  a  Master’s  degree  in  Marine  Science  from  the 
College  of  William  and  Mary,  Williamsburg,  Virginia.  She  has  worked  in 
underwater  acoustics  since  1977.  From  1985  until  joining  SACLANTCEN 
in  1997,  she  worked  at  the  Space  and  Naval  Warfare  Systems  Center,  San 
Diego,  CA,  as  an  Office  of  Naval  Research  program  manager  in  Multistatic 
Active  Surveillance.  During  her  tenure  at  SSC,  she  was  awarded  the  US 
Department  of  the  Navy  award  for  Meritorious  Civilian  Service  in  the  field 
of  Multistatics  in  1996.  She  also  received  two  awards  for  outstanding 
contributions  to  the  USN  Critical  Sea  Test/Low  Low  Frequency  Active 
program  (1995,  1996).  At  SACLANTCEN,  she  has  worked  on  the  Low 
Frequency  Active  program  and  currently  is  leader  of  the  Sound, 
Oceanography  and  Living  Marine  Resources  (SOLMR)  project.  The  focus 
of  her  research  is  to  understand  the  effects  of  anthropogenic  noise  in  the 
marine  environment  in  support  of  SACLANTCEN’s  Acoustic  Risk  Mitigation 
Policy.  Ms.  D’Amico  is  the  scientist  in  charge  of  a  multi-year,  multinational 
at-sea  measurement  program,  “Sirena”,  which  employs  Dual  Use 
Technology  to  acoustically  detect  and  localize  cetaceans.  She  is  responsible 
for  the  functional  design  of  the  SOLMR  data  base  which  contains  cetacean 
sighting  and  stranding  information  for  the  Mediterranean  Sea. 


i 
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Joseph  E.  Bondaryk  received  his  Ph.D.  in  Oceanic  Engineering  from 
MIT  and  WHOI  in  1994.  Subsequently,  he  was  a  Research  Engineer 
and  Lecturer  in  the  MIT  Department  of  Ocean  Engineering  and  a 
Principle  Scientist  at  Engineering  Technology  Center,  Inc.  He  has  over 
twelve  years  of  practical  government  and  corporate  experience  with 
radar  and  sonar  systems  in  the  areas  of  signal  and  array  processing, 
ocean  and  structural  acoustics  and  target  detection.  His  current  project 
is  marine  mammal  risk  mitigation  for  high  power  sonar  systems. 
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Ship  Management  Office, 

NRV  Alliance  and  T  Boat  Manning 


R.V  Alliance 


The  outstanding  operational  availability  of 
the  Alliance  allowed  participation  in  8 
major  scientific  experiments  during  1999, 
including  testing  for  the  Mammal  Risk 
Mitigation  Policy.  As  all  operations  took 
place  in  the  Mediterranean,  the  ship 
remained  on  heightened  security  alert 
status  throughout  the  year,  due  to  NATO 
operations  in  Kosovo. 

Notwithstanding  an  average  of  200 
experimental  sea  days  per  annum  since 
1989,  the  vessel  is  in  excellent  condition, 
with  no  signs  of  steel  plate  consumption 
at  recent  ultrasonic  tests.  Alliance,  with 
continued,  directed  investment,  will 
operate  successfully  for  twenty  years. 

The  Ship  Management  Office  has 
attempted,  within  the  limitations  of 
resources,  to  heighten  the  marketing 
profile  for  Alliance  by  attending  two 
defence  related  exhibitions:  Underwater 
Defence  Technology  (UDT),  Nice  (France) 
in  June  and  Defence  Systems  and 
Equipment  International  (DSEi),  Chertsey 
(UK),  in  September, 


The  success  of  the  marketing  strategy  may  be  judged  by  the  six-day  charter  to  the  Italian 
Navy  during  1999  and  the  34  charter  days  in  2000,  bringing  the  total  number  of  charters 
and  charter  days  since  1997,  when  the  chartering  programme  began,  to  10  and  123 
respectively. 

SMO  continues  to  participate  in  and  contribute  to  the  Research  Vessel  Operators’  Committee 
(RVOC);  this  group  comprises  a  significant  number  of  the  UNOLS  (USA)  and  other  national 
marine  managers  of  research  vessels  and  is  a  professionally  illuminating  and  knowledgeable 
forum  which  meets  at  least  once  a  year. 
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T-Boat  Manning  replacement 

The  pursuit  of  funding  to  design  and  construct  a  replacement  for  the  46  year  old  T-Boat 
Manning  has  met  with  success.  As  a  result,  a  significant  proportion  of  SMO  resources  have 
been  devoted  to  developing  Mission  Profile  Statements,  Requirements  Document  and 
Indicative  General  Description  for  the  proposed  vessel.  Following  an  international  competition 
the  Design  and  Project  Management  contract  was  placed  with  the  Consulting  Naval  Architect 
company  Burness,  Corlett  and  Partners.  A  period  of  six  months  dedicated  to  design  and 
specification  will  precede  a  further  International  Competition  for  the  construction  phase  in 
October  2000.  Delivery  is  scheduled  for  the  final  quarter  of  2001. 


CREW:  To  be  determined 

DIMENSIONS:  27.5m  /  8.5m  /  2.4m  /  260  tons 
RANGE:  1500nm  at  11  knots 

PROPULSION:  Diesel  electric  dynamic  positioning,  one  silent  condition  0-5  knots 
STATUS:  Italian  public  vessel  NATO  operated  classed  by  ABS  and  RINa 


Chris  Gobey  entered  the  Britannia  Royal  Naval  College,  Dartmouth,  in  1961. 
His  subsequent  career  included  a  period  of  loan  service  with  the  Royal 
New  Zealand  Navy  and  four  years  as  commanding  officer  of  HMS  Hecate 
in  the  Falklands,  South  Georgia  and  Antarctica.  He  was  deputy  leader  of  a 
joint  services  expedition  to  the  northern  ice  cap  of  Chile.  He  has  been 
Head  of  the  Ship  Management  Office  at  SACLANTCEN  since  1986. 
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The  site  of  the  new  pier. 
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B  Annual  Bibliography  of  SACLANTCEN 
reports  with  abstracts 


In  addition  to  8  journal  papers,  52  conference  presentations,  4  CD-ROMs  and  one  book 
(listed  under  the  appropriate  project),  the  following  SACLANTCEN  reports  were 
published 


SM-359  _ 

Hermand,J.-P. 

Broad-band  geoacoustic  inversion  in  shallow  water  from  waveguide  impulse  response  measurements  on  a  single 
hydrophone:  theory  and  experimental  results. 

Performance  models  such  as  the  Allied  Environmental  Support  System  (AESS)  require  environmental  input  parameters 
such  as  the  geoacoustic  properties  of  the  battlespace.  The  models  ingest  the  input  parameters  from  databases  to 
predict  the  AS  W  and  MCM  performance  in  the  operating  area.  These  databases  have  a  scarcity  of  data,  particularly 
in  littoral  environments.  A  Rapid  Environmental  Assessment  (REA)  technique  has  been  developed  which  allows  for 
the  determination  of  the  geoacoustic  properties  without  resorting  to  large  acoustic  receiving  apertures.  The  method  is 
based  on  reception  of  broadband  acoustic  signals  and  modelling  the  littoral  environment  between  the  acoustic  transmitter 
and  a  signal  receiving  hydrophone.  In  addition,  the  technique  has  been  extended  to  drifting  sonobuoys  such  that  the 
bottom  properties  of  the  battlespace  can  be  determined  remotely.  This  report  presents  the  measurement  and  modelling 
results  and  demonstrates  that  the  bottom  features  can  be  resolved. 


SM-360 _ _ _ 

Mozzone,  Lv  Berni,  A.,  Guerrini,  P. 

Long  range,  large  throughput  radio  data  link  for  DUSS  (Deployable  Underwater  Surveillance  Systems). 

Two  mono-directional  radio  telemetry  systems  are  described  operating  from  a  buoy  to  NRV  Alliance  at  frequencies  of 
0.4  and  2.28  GHz.  Digital  data  at  2  Mbps  were  transmitted  close  to  the  sea  surface,  collecting  information  on  error 
statistics  and  propagation  loss  versus  buoy  distance,  antenna  height  and  radio  parameters.  A  candidate  system  was 
configured  for  both  frequency  bands  and  the  goal  of  10  n. mi  range  was  achieved.  Field  tests  were  supported  by 
computer  simulation  for  validation  and  a  better  insight  into  the  results.  An  additional  test  assessed  the  performance  of 
a  low-power,  full  duplex,  spread-spectrum  radio  link,  operating  at  the  data  rate  of  128  kbps,  up  to  3.5  n.mi.  The 
experiments  and  conclusions  provide  useful  input  to  the  design  of  a  Deployable  Underwater  Surveillance  System 
( DUSS)  for  scientific  and  operational  purposes. 


SM-362  _ _ 

Tesei,  A.,  Maguer,  A.,  Fox,W.J.L.,  Schmidt,  H. 

Measurements  of  acoustic  scattering  from  partially  and  completely  buried  spherical  shells. 

Recent  work  at  NATO  SACLANT  Undersea  Research  Centre  has  concentrated  on  investigating  the  use  of  low  frequency 
sonars  (2-20  kHz)  in  order  to  better  exploit  scattering  features  of  buried  targets  that  can  contribute  to  their  detection 
and  classification. 

Part  of  the  recent  GO  AT 98  experiment  performed  off  the  island  of  Elba,  involved  controlled  monostatic  measurements 
of  scattering  by  spherical  shells  which  were  partially  and  completely  buried  in  sand,  as  well  as  suspended  in  the  water 
column.  Analysis  is  mainly  addressed  to  a  study  of  the  effect  of  burial  on  the  dynamics  of  backscattered  wave  families, 
which  can  be  clearly  identified  in  the  target  responses.  Data  interpretation  results  are  in  good  agreement  with  theory. 
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SM-363 


LePage,  K.D.,  Schmidt,  H. 

Spectral  integral  representations  of  scattering  from  volume  inhomogeneities 

In-situ  measurements  of  scattering  strength  can  be  obtained  by  analyzing  the  early-time,  high  angle  reverberation 
from  bottom  and  sub-bottom  features.  In  order  to  provide  insight  into  the  mechanisms  which  cause  bottom  reverberation 
and  to  their  distinguishing  characteristics,  it  is  necessary  to  have  a  capability  for  modelling  both  the  rough  surface  and 
the  volume  scattering  mechanisms.  For  high-angle,  early-time  backscatter ;  the  most  appropriate  approach  is  to  use  a 
spectral  integral  representation,  which  naturally  includes  the  continuous  spectrum  important  for  this  angular  regime. 
A  rough  surface  scattering  theory  developed  earlier  in  this  framework  has  provided  important  insights  into  wave 
scattering  and  penetration  physics  at  the  seafloor.  Here  a  consistent  representation  for  the  subbottom  scattering  is 
developed  and  implemented  in  the  same  spectral  integral  code,  OASES. 


SM-364 _ 

Holland,  C.W.,  Hollett,  R.,Troiano,  L. 

Bottom  scattering  measurements  in  shallow  water. 

Sonar  performance  predictions  of  reverberation  in  shallow  water  rely  upon  good  estimates  of  the  bottom  scattering 
strength.  However,  little  is  understood  about  bottom  scattering  in  shallow  water  in  the  frequency  range  400  -  4000  Hz, 
particularly  its  dependency  upon  frequency  and  its  relationship  to  the  physical  properties  of  the  seafloor.  In  order  to 
address  these  issues,  new  measurement  techniques  have  been  developed  to  probe  the  frequency  and  angular 
dependency  of  bottom  scattering  strength.  The  measurement  techniques  also  appear  to  be  capable  of  revealing  the 
physical  mechanisms  that  give  rise  to  the  scattering.  Several  experimental  techniques  will  be  described,  including  use 
of  coherent  and  incoherent  sources  (lightbulbs).  The  general  experimental  approach  is  also  described  which  includes 
auxiliary  acoustic  and  geoacoustic  measurements  designed  to  allow  exploration  of  the  relationship  between  bottom 
scattering  and  the  physical  properties  of  the  bottom.  Measurement  and  modelling  results  for  two  shallow  water  sites 
are  presented.  Atone  site,  the  scattering  appeared  to  arise  from  at  or  near  the  water-sediment  interface.  At  the  other 
site,  scattering  from  a  25  m  sub-bottom  horizon  is  clearly  apparent  in  the  data  at  and  below  1800  Hz. 


SM-365 _ 

Fioravanti,S. 

An  efficient  SAS  processor 

This  work  describes  the  design  and  the  implementation  of  an  efficient,  high-speed  SAS  processor  able  to  receive  the 
navigation  system  data  for  an  accurate  generation  of  the  synthetic  image.  The  SAS  beamformer  is  a  pixel-based  one 
working  in  the  time  domain  with  real  or  complex  data  (in-phase  and  in-quadrature).  The  adopted  approach  allows  the 
generation  of  physical  images  during  the  process.  The  report  explains  the  technique  developed  for  a  fast  phase 
reconstruction  with  complex  beam  forming. 

The  SAS  processor  makes  use  of  Taylor  approximation  for  an  improvement  in  the  computational  speed.  The  analysis 
of  the  region  of  validity  of  the  approximation  is  mathematically  computed  depending  on  the  order  of  the  Taylor  expansion. 

The  algorithm  has  been  implemented  in  standard  ANSI -C  language  and  intensively  tested  on  an  HP-C180  workstation 
with  simulated  and  real  data. 

An  accurate  analysis  of  the  computational  efficiency  is  reported. 


SM-366 _ 

Hermand,  J.-P.,  Nascetti,  P.,  Cinelli,  F. 

Inverse  acoustical  determination  of  photosynthetic  oxygen  productivity  of  posidonia  seagrass. 

As  part  of  geoacoustic  inversion  experiments  (Yellow  Shark)  in  the  Giglio  basin,  off  the  west  coast  of  Italy,  low  frequency, 
broad-band  propagation  measurements  were  performed  in  the  winter  and  spring  of  1995  over  a  dense  and  extensive 
prairie  of  Posidonia  seagrasses  which  surrounds  Scoglio  d’ Africa,  a  minor  island  of  the  Tuscan  Archipelago.  The 
purpose  of  the  measurements  was  to  determine  the  applicability  of  model-based,  geoacoustic  inversion  techniques 
developed  for  marine  sediments  to  the  monitoring  of  oxygen  synthesis  by  Posidonia. 

A  dual-flextensional  acoustic  projector  and  a  4-element  vertical  receiving  array  were  positioned  at  1541  m  distance  in 
an  isobath  (25  m  water  depth),  dense  and  homogeneous  part  of  the  Posidonia  prairie.  The  waveguide  impulse  response 
was  measured  during  one  day  by  1~min  repeated  and  alternated  transmissions  of  3  s  chirp  signals  with  frequency 
bands  0. 1-0.9  kHz  and  0.8-1. 6  kHz. 

The  water  sound  speed  profiles  calculated  from  repeated  CTD  measurements  were  slightly  downward  refracting  and 
exhibited  little  temporal  variability  except  for  mild  surface  heating  in  the  afternoon.  Contemporaneous  oxygen  and 
CTD  profiles  as  a  function  of  daytime  and  season  were  obtained  in  1997  to  support  the  present  study. 
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In  this  paper  acoustic,  solar  radiation,  oxygen  and  CTD  data  are  analyzed  and  discussed.  The  analysis  shows  strong 
correlation  between  photosynthesis  and  the  impulse  response  of  the  acoustic  waveguide.  The  most  evident  feature  is 
an  abrupt  and  marked  change  of  attenuation  and  time  dispersion  characteristics  at  the  onset  of  photosynthesis. 
Frequency-  and  depth-dependent  rapid  variations  of  received  energy  (2-5  dB)  and  time  spread  (3-10  ms)  are  observed. 
The  time  of  occurrence  and  rate  of  change  of  these  variations  are  consistent  with  solar  time  and  oxygen  concentration 
measured  in  situ.  The  phenomena  is  attributed  to  bubbles  of  photosynthetic  oxygen  formed  on  the  Posidonia  leaf 
blades.  The  bubble  layer  creates  an  absorbing,  dispersive,  low-speed,  thin  parallel  waveguide  which  modifies  interaction 
of  acoustic  energy  with  the  Posidonia  “matte”.  It  is  demonstrated  that,  after  sunrise,  low-order  modes  begin  to  travel  in 
the  seagrass  layer,  absorbing  a  portion  of  acoustic  energy  from  the  main  waveguide.  A  similar  effect  can  occur  in  a 
near-surface  waveguide  when  supersaturation  (undissolved  oxygen)  conditions  obtain. 

Modeling  results  indicate  that  the  inverse  problem  of  determining  gas  and  oxygen  void  fractions  in  the  seagrass  layer 
could  be  solved.  Parameters  such  as  surface  density  and  photosynthetic  efficiency  of  Posidonia  can  be  derived  from 
the  variations  of  inverted  void  fractions. 

These  results  may  be  applicable  to  the  monitoring  of  the  state  of  health  of  Posidonia  and  other  seagrasses  in  the 
Mediterranean  and  other  oceans 


SM-367 _ 

Berni,  A.,  Mozzone,  L. 

The  application  of  spread-spectrum  communications  to  REA  tactical  networks  and  deployable  underwater  surveillance 
systems. 

The  present  document  is  a  part  of  an  ongoing  study  on  advanced  communication  techniques  in  support  of  Rapid 
Environmental  Assessment  (REA)  and  Deployable  Undersea  Surveillance  Systems  (DUSS).  It  aims  at  the  definition 
of  data  transmission  architectures  for  both  scientific  data  acquisition  at  SACLANTCEN  and  for  operational  concept 
demonstration. 

Emerging  methodologies  for  Anti-Submarine  Warfare  (ASW)  and  Rapid  Environmental  Assessment  (REA)  increasingly 
rely  on  communication  technology,  in  order  to  exchange  information  with  other  naval  units  or  military  commands 
ashore. 

The  specific  requirements  (such  as  ranges  and  transmission  data  rates)  are  addressed.  They  vary  from  2.4  kbps  to  2 
Mbps  and  from  10  to  20  n.mi.  In  addition  to  that,  all  military  applications  have  common  requirements  in  terms  of 
reliability,  availability  and  security.  The  eligibility  of  classical  and  spread-spectrum  radio  communication  techniques  to 
the  fulfillment  of  such  requirements  is  discussed  and  shown.  Performance  is  estimated  and  compared. 

Spread-spectrum  techniques  offer  such  features  as  interference  rejection,  anti-jam  capability  and  low-density  power 
spectra  for  covert  operations:  field  tests  of  spread-spectrum  equipment  have  successfully  been  conducted  in  1998 
during  SACLANTCEN  experiments  GOATS  98  and  DUSS  98.  The  availability  of  commercial  off-the-shelf  (COTS) 
products  including  sophisticated  communication  protocols  also  represents  a  relevant  issue  for  scientific  applications. 
Such  features  as  the  creation  of  a  multi-point,  error-free,  packet  switching  Local  Area  Network  (LAN)  deserve  to  be 
carefully  investigated  in  their  impact  on  DUSS  and  REA  applications. 


SM-368  '  “ 

Tesei,  A.,  Pinto,  M. 

Application  of  aided  inertial  navigation  system  to  synthetic  aperture  sonar  micronavigation. 

In  order  to  improve  navigational  accuracy,  critical  to  successful  SAS  processing,  a  high-quality  Aided  Inertial  Navigation 
System  (AINS)  is  investigated,  which  consists  of  an  inertial  navigation  system  (INS),  Doppler  velocity  log  and  intermittent 
DGPS  fixes,  the  data  of  which  are  fused  by  an  extended  Kalman  filter.  Although  the  strapdown  INS  recently  acquired 
at  SACLANTCEN  is  one  of  the  best  available,  the  accuracy  of  positioning  remains  insufficient,  in  particular  for  high- 
frequency  SAS. 

Through  a  navigational  simulator/estimator  provided  by  NDRE  (Norwegian  Defence  Research  Establishment),  the 
analysis  of  expected  performances  of  the  AINS  has  been  performed.  In  particular,  the  short-term  residual  errors  of 
position  displacement  and  attitude  have  been  analyzed  in  terms  of  dominant  components  and  their  sources.  Preliminary 
results  show  that  the  attitude  error  is  estimated  to  be  negligible  due  to  the  very  low  values  of  gyro  bias  and  white  noise. 
The  position  error  is  shown  to  be  dominated  by  an  evident  quadratic  drift  due  to  accelerometer  biases. 

Aiding  the  AINS -based  navigation  with  SAS  motion  compensation  techniques  is  foreseen  to  significantly  reduce  this 
critical  error  component.  The  design,  development  and  test  of  the  combined  SAS- AINS  architecture  will  be  the  main 
objective  of  next  activities. 
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SM-369 


Akal  ,T.,  Di  iorio,  D.,  Guerrini,  P.,  Boni,  R,  Cavanna,  A.,  Stoner  R.,  Ferla,  C.,  Kuperman,W.A.,  Hodgkiss,  W.,  Song, 
H.C.,  Edelmann,  G.,  Kim,  S. 

Focused  acoustic  field  at  3.5  kHz:  data  report  for  FAF-99  sea  trial. 

Two  low  frequency  (-450  Hz)  Focused  Acoustic  Field  (FAF)  experiments  previously  were  carried  out  in  1996  and 
1997 in  - 125  m  water  adjacent  to  Formiche  di  Grosseto  (a  small  island  approximately  1 00  miles  SIV  of  SACLANTCEN). 
These  experiments  demonstrated  that  FAF  is  both  feasible  and  stable  at  low  frequencies  in  shallow  water  and  that 
focusing  of  the  retransmitted  energy  is  possible  at  ranges  of  at  least  30  km.  Also  demonstrated  was  the  ability  to  shift 
the  range  of  focus  to  ranges  other  than  that  of  the  probe  source  by  a  simple  method  involving  a  frequency  shift  of  the 
received  time  series  prior  to  retransmission. 

As  an  outgrowth  of  the  successful  low  frequency  phase  conjugation  experiments ,  a  high  frequency  (3.5  kHz)  FAF 
experiment  was  carried  out  within  the  Joint  Research  Project  between  SACLANTCEN  and  MPL  in  July  1999. 
Experiments  were  carried  out  adjacent  to  both  Formiche  di  Grosseto  and  N  of  Elba  Island.  The  former  provided  a  link 
to  our  previous  experiments  while  the  latter  provided  a  brief  opportunity  to  explore  a  new  environment.  During  this 
experiment  we  demonstrated  the  high  frequency  FAF  at  ranges  out  to  14  km  in  both  flat  and  sloping  coastal  regions, 
the  short-term  temporal  and  spatial  properties  of  the  focal  region  in  two  different  bottom  types  and  the  initial 
demonstration  of  phase-conjugation  processing  in  acoustic  communications. 

In  this  memorandum  we  present  the  experimental  procedures  and  acoustic  and  environmental  data  collected  during 
the  July  1999  experiment.  An  additional  CD  containing  the  data  is  provided. 
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Noseworthy,  D. 

Seismic  stratigraphy  of  the  Capraia  Basin,  Northern  Tyrrhenian  Sea. 

The  distribution  and  thickness  of  sediments  within  the  Capraia  Basin,  Northern  Tyrrhenian  Sea,  have  been  mapped 
using  high-resolution  seismic  reflection  data.  The  data  were  studied  within  the  geologic  framework  determined  by 
previous  works  conducted  in  the  Northern  Tyrrhenian  Sea.  This  study  was  especially  influenced  by  the  work  of  Brizzolari 
et  at  (1991).  Collectively  these  studies  show  that  sediment  above  the  basement  in  the  Capraia  Basin  consists  of  at 
least  two  seismic  units,  deposited  since  the  Mid-Upper  Pleistocene. 

Unit  1  is  correlated  with  Brizzolari’s  eta!.,  (1991)  unit  B  and  unit 2  with  his  unit  C.  Unit  1  lies  unconformably  below  unit 
2  and  covers  a  vast  area  within  the  Capraia  Basin.  As  determined  by  Brizzolari  et  al.,(1 991),  the  deposition  of  unit  B 
and  unitl  was  strongly  influenced  by  Mid-Upper  Pleistocene  glacio-eustatic  fluctuation.  Subsequent  rise  and  fall  of 
sea  level  within  shallower  regions  of  the  basin,  especially  between  1 15- 100  m  water  depth,  resulted  in  the  deposition 
of  stacked  sedimentary  successions  (ofunit2)  which  display  strong  erosional  character  and  discontinuity.  Basinward, 
unit  1  thickens  drastically,  from  10 -80  ms  twtt,  into  a  package  of  aggregated  to  slightly  progradational  reflectors. 
Conversely,  unit  2  thins  basinward.  It  is  limited  geographically  to  the  northwestern  region  of  the  Capraia  Basin.  As 
determined  by  Brizzolari  etal.,  (1991),  the  geometry  of  unit  2,  is  strongly  influenced  the  depositions  setting.  Unit2is 
composed  of  a  thin  (1- 15  ms  twtt)  package  of  stratified  reflectors,  which  taper  westward.  Studies  have  shown  that 
successions  correlating  with  unit 2  are  predominantly  composed  of  sediment  transported  into  the  basin  by  modern 
day  feeding  sources,  with  dominant  influence  coming  from  the  Cecina  River,  located  northeast  of  the  study  area. 
Brizzolari  etal.,  (1991)  interpreted  unit  2  to  be  the  result  of  present-day  highstand  sedimentation  in  the  Capraia  Basin, 
deposited  since  the  Holocene. 


SM-372 _ _ _ 

Acunto,  S.,  Lyons,  A.R,  Pouliquen,  E, 

Characteristics  of  the  Mediterranean  seagrass  Posidonia  oceanica  of  contributing  to  high  frequency  acoustic  scattering. 

Posidonia  oceanica  meadows  are  the  most  important  ecosystem  for  the  life  cycle  of  coastal  Mediterranean  benthos 
having  a  fundamental  role  in  the  primary  production  of  the  neritic  system  and  a  decisive  influence  on  other  vegetation 
and  animal  communities.  The  P.  oceanica  meadows  are  presently  undergoing  a  slow,  but  progressive  regression 
throughout  the  Mediterranean  and  generally,  the  most  common  cause  of  regression  is  enhanced  turbidity  and  the 
consequent  reduction  of  water  transparency.  The  first  step  towards  the  safeguarding  of  such  important  coastal 
environments  is  that  of  defining  their  extent  and  state  of  health.  Echographic  survey  allows  general  maps  to  be 
obtained,  but  accurate  seafloor  characterization  requires  knowledge  of  the  acoustic  properties  of  seabed  vegetation. 

In  this  report  we  present  some  of  the  microscopic  characteristics  of  P.  oceanica  leaves.  The  presence  of  air  channels 
running  the  length  of  each  leaf  blade  and  the  presence  of  cellulose  that  is  the  main  constituent  of  cell  walls  are 
assumed  to  be  extremely  important  for  explaining  the  strong  high-frequency  scattering  of  Posidonia  meadows.  A 
preliminary  study  has  been  conducted  that  represents  a  first  attempt  to  quantify  the  air  within  the  leaves  of  P.  oceanica 
and  to  consider  it  as  a  function  of  the  plant  life  cycle.  An  order  of  magnitude  of  the  volumes  of  the  different  elements 
considered  and  their  relative  importance  expressed  as  percentages  has  been  given.  The  values  obtained  have  been 
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considered  with  reference  to  the  seasonal  changes  of  the  phenological  variables  from  which  they  depend  and  also  to 
the  main  structural  variable:  the  density  of  leaves. 

A  future  step  from  the  present  work  will  be  the  development  of  a  model  for  the  acoustic  scattering  of  P.oceanica 
derived  from  one  developed  for  gassy  sediments.  The  final  result  could  be  a  fundamental  tool  to  improve  the  present 
capabilities  of  acoustic  instruments  and  make  it  possible  from  an  ecological  point  of  view  to  obtain  fast  and  large- 
scale  information  not  only  on  the  extension  of  the  meadows,  but  also  on  the  density  and  height  of  the  plants.  On  the 
other  hand ,  from  an  operational  point  of  view  it  could  improve  the  performances  of  mine-hunting  sonars. 


SR-310  - 

Zerr,  B.,Tesei,  A.,  Maguer,  A.,  Fox,  W.L.J.,  Fawcett,  J.A. 

A  classification  technique  combining  aspect  dependence  and  elastic  properties  of  target  scattering. 

Discrimination  between  man-made  and  natural  underwater  objects  and  between  man-made  objects  of  different 
characteristics  are  the  key  objectives  of  target  classification.  The  current  approach  is  mainly  based  on  the  analysis  of 
the  target  signature  imaged  by  high  resolution  (typically  <20  cm)  high  frequency  (typically  <100  kHz)  sonars.  To 
estimate  the  potential  of  alternative  classification  schemes  based  on  more  detailed  acquisition  of  target  echo  features, 
SACLANTCEN  has  investigated  low  frequency  ( 1 0  kHz)  sonar  systems  with  high  fractional  (or  relative)  bandwidth  (2- 

The  classification  method  relies  on  two  features  of  the  target  echo:  aspect  dependence  and  elastic  scattering.  The  two 
dimensional  reflection  map,  reconstructed  from  multiple  aspects,  serves  as  the  basis  for  pre-classification.  For  objects 
of  external  shape  recognized  as  man-made  at  the  pre-classification  stage,  resonance  information  is  extracted  by 
autoregressive  spectral  estimation  techniques  and  further  processed  for  particular  aspects  (e.g.  broadside).  The 
analysis  of  the  resonance  scattering  provides  an  estimate  of  some  geometrical  and  physical  target  parameters  (i.e. 
shell  wall  thickness  and  material  and  inner  fluid  properties)  and  can  be  improved  by  introducing  the  object  dimensions, 
estimated  from  the  reflection  map.  This  report,  which  describes  the  classification  methodology  and  the  results  obtained 
with  steel  cylinders,  demonstrates  the  potential  of  the  method. 
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Osier,  J.,  Algan,  O. 

A  high  resolution  seismic  sequence  analysis  of  the  Malta  Plateau. 

A  high  resolution  seismic  study  of  the  Malta  Plateau  was  conducted  to  characterize  the  Quaternary  sediment  layers, 
mapping  their  regional  extent  and  determining  the  environmental  conditions  during  their  deposition.  Seismic  stratigraphic 
analysis  indicates  that  there  are  six  seismic  units  (or  depositional  sequences)  bounded  by  erosional  unconformities. 
The  basement  surface,  that  is  an  erosional  truncation  in  some  places,  indicates  that  a  hiatus  probably  occurred  during 
Messinian  to  Early  Pliocene.  The  terminations  of  all  the  sequence  boundaries  onlap  the  acoustic  basement  in  a 
landward  direction,  revealing  an  overall  relative  sea  level  rise.  The  termination  of  the  internal  strata  (reflectors)  in  Units 
2,3  and  4  have  an  oblique  parallel  pattern  and  toplap  at  their  upper  surface,  characteristic  of  periods  of  sea  level 
stillstand  and/ora  decrease  in  the  rate  of  sea  level  rise.  It  is  anticipated  that  this  depositional  history  would  cause  the 
physical  properties  of  the  sediment  in  the  immediate  vicinity  of  the  erosional  unconformities  to  differ  from  that  in  the 
adjacent  sequences.  Within  a  given  sequence  and  with  respect  to  the  other  sequences,  it  is  anticipated  that  the 
physical  properties  of  the  sediment  are  relatively  homogeneous,  with  the  exception  of  the  uppermost  sequence, 
Seismic  Unit  1. 

Carbonate  buildup  structures  are  observed  at  the  boundary  between  Units  2  and  3  in  the  northwest  part  of  the  study 
area  and  imply  a  very  shallow  water  depositional  condition  during  their  formation.  The  uppermost  layer,  Unit  1,  is  the 
most  extensive  geographically  and  gets  thicker  towards  the  coast  of  Sicily.  The  parallel  configuration  of  its  internal 
reflectors  and  the  progradational  character  of  the  unit  suggest  that  it  is  comprised  of  coastal  sediments  deposited 
during  the  Holocene  highstand  of  sea  level.  In  the  shelf-slope  transition,  sequences  are  deformed  by  mass  movements 
that  might  be  attributable  to  tectonic  activity  in  the  region.  Several  basement  outcrops  were  observed  in  the  study 
area.  Some  have  northwest-southeast  trending  lineations,  as  observed  on  Sicily,  while  others  are  isolated  structures, 
possible  dikes  or  salt  domes. 


SR-312  _ 

Watermann,  J. 

The  magnetic  dipole  moment  of  RV  Alliance. 

A  SACLANTCEN  magnetometry  sea  trial  was  performed  in  September  1998  in  the  northern  Tyrrhenian  Sea  in  the 
vicinity  of  the  Formiche  di  Grosseto  islands.  One  of  the  trial  objectives  concerned  the  measurement  of  the  magnetic 
dipole  moment  of  the  NATO  research  vessel  Alliance  with  an  array  of  fixed  Ocean  Bottom  Magnetometers.  The 
Alliance  sailed  with  constant  velocity  along  three  specified  tracks  in  reciprocating  directions.  The  magnetic  signature 
of  the  Alliance  was  subsequently  extracted  from  the  magnetic  field  data  recorded  during  the  experiment.  From  the 
magnetic  signatures  taken  from  six  successive  runs,  quantitative  magnetic  dipole  models  were  derived.  A  comparison 
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between  the  magnetic  signatures  obtained  from  northwestward  and  southeastward  runs  yields  estimates  of  the 
permanent  and  induced  magnetic  dipole  moments  in  the  horizontal  plane.  In  the  vertical  direction ,  only  the  total 
magnetic  dipole  moment  can  be  determined.  The  positions  of  the  OBMs,  uncertain  to  some  extent,  were  fictitiously 
varied  such  that  the  horizontal  components  of  the  permanent  magnetic  dipole  moment  and  the  vertical  component  of 
the  total  dipole  moment  remained  constant  during  the  experiment.  The  induced  horizontal  dipole  moment  tended  to 
increase  steadily  and  monotonically  in  a  northward  direction.  Explanations  of  this  phenomenon  are  discussed. 


SR-315 _ _ _ _  - _ 

Maguer,  A.,  Fox, W.L.J.,  Zerr,  B.  Jesei,  A.,  Bovio,  E.,  Fioravanti,  S. 

Buried  mine  detection  and  classification  (Research  Summary  1 996-1 999) 

The  applicability  of  low-frequency  sonar  (2-16  kHz)  to  buried  mine  detection  has  been  investigated.  Experiments 
were  performed  on  sound  penetration  into  sediment,  buried  target  detection  and  broadband  multiple  aspect  classification. 
The  results  of  the  experiments  are  given  in  this  report  and  compared  with  modelled  results. 

One  of  the  main  results  is  our  success  in  understanding  the  physical  mechanisms  contributing  to  subcritical penetration 
into  sediment.  It  has  been  demonstrated  that  the  evanescent  wave  was  dominant  in  the  lower  frequencies  [2-5  kHz]  of 
our  bandwidth  of  interest  [2-1 6  kHz].  Roughness  scattering  dominates  at  higher  frequencies  (above  5  kHz)  for  our 
bottom  type  (RMS  roughness  1.5  cm,  cross-ripple  correlation  length  25  cm).  Although  roughness  scattering  has  been 
shown  to  be  one  mechanism  for  explaining  “anomalous" penetration  into  sediment,  its  potential  for  detection  and 
classification  of  buried  objects  is  unclear  due  to  the  low  level  and  the  lack  of  coherence  of  the  received  signals .  It  is 
demonstrated  that  sound  speed  variation  with  frequency,  could  exist  for  permeable  sandy  bottoms,  which  could 
influence  the  design  of  a  buried  mine  sonar. 

The  detection  of  buried  targets  is  shown  to  be  very  effective  at  above  the  critical  angle.  At  subcritical  angles,  detection 
becomes  difficult.  Significant  gains  in  signal-to-reverberation  radio  below  critical  angle  were  obtained  either  by 
emphasizing  a  relatively  narrow  band  of  frequencies  at  the  lower  end  of  the  transmitted  bandwidth  (below  3  kHz)  or  by 
using  a  larger  physical  array  or  synthetic  array  processing  which  improve  the  sonar  resolution.  Simulations  have 
shown  that  lower  frequencies  (of  the  order  of  0. 5- 1  kHz)  are  essential  to  the  detection  of  buried  targets  at  low  grazing 
angles  and  that  the  de3tection  at  those  frequencies  will  only  be  effective  for  shallow  buried  targets. 

A  method  based  on  multiple-aspect  target  echo  analysis  in  time  and  frequency  domains,  which  considers  the  rigid 
and  resonance  responses  is  presented.  Its  potential  was  demonstrated  in  simulation  and  on  real  data  (exercise  mine, 
cylinders  and  rocks)  for  proud  targets  and  for  buried  spheres 


SR-316 _ _ 

Siderius,  M.,  Snellen,  M.,  Simons,  D.,  Onken,  R. 

Environmental  assessment  in  the  Strait  of  Sicily:  measurements  and  analysis  techniques  for  determining  bottom  and 
oceanographic  properties. 

In  October  1997,  the  EnVerse  97  shallow  water  acoustic  experiments  were  jointly  conducted  by  SACLANT  Centre, 
TNO-FEL  and  DERA  off  the  coast  of  Sicily,  Italy.  The  primary  goal  of  the  experiments  was  to  determine  the  seabed 
properties  through  inversion  of  acoustic  data.  Using  a  towed  source,  the  inversion  method  is  tested  at  different  source- 
receiver  separations  in  an  area  with  a  range-dependent  bottom.  The  sources  transmitted  over  a  broad-band  of 
frequencies  (90-600  Hz)  and  the  signals  were  measured  on  a  vertical  array  of  hydrophones.  The  acoustic  data  were 
continuously  collected  as  the  range  between  source  and  receiving  array  varied  from  0.5-6km.  An  extensive  seismic 
survey  was  conducted  along  the  track  providing  supporting  information  about  the  layered  structure  of  the  bottom  as 
well  as  layer  sound  speeds.  The  oceanic  conditions  were  assessed  using  current  meters,  satellite  remote  sensing, 
wave  height  measurements  and  casts  for  determining  conductivity  and  temperature  as  a  function  of  water  depth. 
Geoacoustic  inversion  results  taken  at  different  source/receiver  ranges  show  seabed  properties  consistent  with  the 
range  dependent  features  observed  in  the  seismic  survey  data.  These  results  indicate  that  shallow  water  bottom 
properties  may  be  estimated  over  large  areas  using  a  towed  source  fixed  receiver  configuration. 


SR-317 _ _ 

Mozzone,  L.,  Lorenzelli,  P. 

Deployable  underwater  surveillance  systems. Target  localization  with  multiple  sonar  receivers. 

Deployable  Underwater  Surveillance  Systems  (DUSS)  are  a  network  of  small  multistatic  transmitter /receiver  sonar 
nodes.  This  study  analyzes  the  contact  localization  capabilities  of  DUSS  in  term  of  range,  time  and  bearing  error.  This 
information  is  used  in  Monte  Carlo  simulations  to  estimate  the  accuracy  of  multistatic  localization  methods  suing  2  or 
3  receivers.  Simulations  are  validated  by  real  data.  Time  -  only  localization  of  active  sonar  echoes  with  2  receivers 
produces  error  estimates  of  150  mat  10  km.  Active  pinger  localization  with  2  receivers  produces  average  errors  of 
83  mat  10  km.  Bearings  -  only  passive  localization  with  2  receivers  produces  average  errors  of 250  mat  10  km.  Buoy 
separation,  buoy  localization  accuracy,  acoustic  travel  time  estimation,  beam  width  and  compass  accuracy  are  the 
most  critical  system  parameters.  The  use  of  three  receivers  further  improves  accuracy. 
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Mozzone,  L.,  Bongi,  S.,  Ziegenbein,  J.  ~~ 

Diversity  in  deployable  underwater  surveillance  systems. 

A  multistatic  network  of  autonomous  active  sonar  transmitters  and  receivers,  such  as  the  DUSS  (Deployable  Underwater 
Surveillance  Systems)  currently  being  studied  at  the  NATO  S AC  LANT  Center,  often  provides  multiple  contacts  of  the 
same  target.  The  present  paper  examines  receiver  diversity  (“spatial  diversity")  on  a  multistatic  system,  where  three 
independent,  spatially  separa ted  sonar  receivers  detect  the  same  target,  insonified  by  the  same  FM  or  CW pulse  (in 
a  symmetrical  configuration)  and  the  contacts  are  merged  together  after  detection. 

Real  data  collected  South-East  of  the  island  of  Elba  (Italy),  during  “DUSS  ’97”  tests  with  SACLANTCEN’s  experimental 
system  are  analyzed.  The  statistical  characteristics  of  measured  background  noise  and  reverberation  are  estimated. 
Signal  to  Noise  Ratio  (SNR)  series  of  echoes  from  a  towed  Echo  Repeater  are  analyzed  and  decomposed  into 
slowly-varying  and  rapidly-varying  components.  The  degree  of  inter-receiver  cross-correlation  is  estimated  and  related 
to  the  corresponding  overall  detection  performance  of  the  multistatic  system.  Experimental  distributions  of  detections 
and  false  contacts  are  computed  and  experimental  Pd  -  Pfa  (Probability  of  detection  -  Probability  of  false  alarm) 
curves  are  extracted.  Receiver  Operating  Characteristic  (ROC)  curves  are  produced  with  Monte  Carlo  simulations  of 
a  classical  statistical  model  and  adopted  as  a  reference. 

The  analysis  is  extended  to  time  diversity,  comparing  single-ping  to  multi-ping  detection.  Frequency  diversity  is  also 
addressed  by  comparing  the  classical  single  -  frequency  CW  pulse  with  a  “ stepped  CW”,  i.e.  a  sequence  of  three 
short  CW  sub-pulses  at  different  frequencies,  combined  together  after  detection. 

A  previous  paper  demonstrated  the  feasibility  of  post-detection  fusion  of  contacts  received  at  different  sensors. 
[  Localization  and  fusion  of  echoes  with  deployable  multistatic  active  sonar:  evaluation  of  feasibility  using  experimental 
data",  L.  Mozzone,  S.  Bongi,  4th  ECU  A  ”98,  Roma]. 

The  performance  enhancement  produced  by  the  three  types  of  diversity  studied  here  is  quantified  in  the  experimental 
data  and  confirmed  by  simulations.  An  increase  of  Pd  varying  between  35%  and  50%  is  found  (with  Pfa  of  10-5.  ..10- 
6)  when  contacts  from  three  independent  receivers  (or  frequencies,  or  pings)  are  combined  together.  The  same 
improvement  is  expressed  in  terms  of  an  equivalent  increase  of  echo  SNR  ranging  from  7  to  8  dB. 


SR-319 _  " 

LePage,  K.D. 

Time  series  variability  in  fluctuating  ocean  waveguides. 

The  variability  of  signals  propagating  through  an  uncertain  sound  speed  structure  is  addressed.  Signals  are  assumed 
to  travel  in  a  narrow  band  adiabatically  in  modes  and  to  experience  fluctuations  in  sound  speed  which  are  characterized 
according  to  the  vertical  and  horizontal  distributions  of  these  fluctuations.  The  sound  speed  fluctuations  are  assumed 
to  affect  only  the  phase  speed  and  the  group  speed  of  the  modes  in  a  perturbative  way.  The  changes  in  the  local  phase 
and  group  speeds  are  expanded  for  small  perturbations  to  the  sound  speed.  Sound  speed  perturbations  are  described 
in  terms  of  their  statistical  characteristics.  Vertically,  the  sound  speed  fluctuations  are  decomposed  into  empirical 
orthogonal  functions  (EOFs),  while  horizontally  they  are  assumed  to  be  correlated  on  some  horizontal  length  scale 
much  smaller  than  the  propagation  ranges  of  interest.  Thus  the  cumulative  phase  and  group  speed  fluctuations  over 
the  propagation  path  are  assumed  to  be  distributed  Gaussian  according  to  the  central  limit  theorem. 

The  framework  outlined  above  is  used  to  derive  the  first  and  second  moments  of  the  signal  envelope  received  over  an 
ensemble  of  ocean  realizations  following  the  distribution  properties  outlined  above.  Since  the  mean  and  variance  of 
the  expected  signal  are  obtained  in  the  time  domain,  the  stability  of  model  arrivals  in  time  can  be  predicted  for  a 
variety  of  different  sound  speed  fluctuation  distributions.  Since  the  phase  and  group  speed  fluctuations  are  linear  in 
identical  terms  involving  the  inner  product  of  the  mode  shape  functions  with  the  EOFs,  the  fluctuations  of  these 
quantities  are  entirely  correlated.  However,  as  the  different  EOFs  express  themselves  differently  on  each  set  of 
propagating  modes,  the  modal  interference  structure  becomes  less  certain  due  to  the  fluctuations.  The  theory  estimates 
this  degree  of  decorrelation  as  a  function  of  the  signal,  waveguide  and  fluctuation  parameters. 

In  order  to  benchmark  the  theory,  the  first  moment  of  the  short  time  average  of  the  signal  intensity  is  also  predicted 
using  realizations  of  propagation  through  an  ensemble  of  sound  speed  fluctuations  consistent  with  the  statistical 
description.  Excellent  agreement  is  found  between  these  self-consistent  Monte-Carlo  estimates  of  the  signal  variability 
and  the  closed  form  expressions. 
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Haralabus,  G.,Capriulo,  E., Zimmer,  W.M.X. 

SWAC  4:  Broadband  data  analysis  using  sub-band  processing. 

The  frequency  dependence  of  broadband  active  detection/localization  is  examined.  The  analysis  is  based  on  1200  Hz 
(2300-3500  Hz)  LFM  signals  acquired  during  the  SWAC  4  sea  trial.  A  sub-band  matched  filter  scheme  is  devised 
according  to  which  a  replica  of  the  transmitted  pulse  is  segmented  into  ten  120  Hz  sub-bands  and  processed 
independently  through  a  matched  filter  detector.  Comparison  of  target  detection  and  ranging  results  indicate  comparable 
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performance  for  all  sub-bands.  However ;  ping-to-ping  variability  of  the  ten  correlator  outputs  suggest  that  the  detection 
performance  may  be  improved  by  employing  incoherent  processing  schemes.  Signal-to-noise  ratio  is  proved  to  be 
controlled  mainly  by  noise  (reverberation  is  the  predominant  noise  source)  rather  than  signal  variations.  The  signal 
intensity  remains  proportional  to  the  distance  between  source  and  receiver  due  to  favorable  propagation  conditions. 
Doppler  effects  and  sub-band  detection  synchronization  problems  which  may  lead  to  performance  degradation  in 
large  time-bandwidth  signal  processing  are  addressed.  A  method  to  estimate  range  rate  (relative  velocity  between 
source  and  receiver)  based  on  single  ping  differential  time  delay  between  sub-band  MF  outputs  is  developed.  This 
intra-ping  technique  is  an  alternative  to  the  standard  inter-ping  method  which  requires  multiping  detection  history. 
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Osier,  J.C.,  Gualdesi,  L.,  Michelozzi,  E.,Tonarelli,  B. 

Piston  coring  capabilities  at  SACLANTCEN:  minimizing  and  assessing  core  disturbance. 

Piston  and  gravity  coring  and  techniques  used  to  collect  samples  of  seabed  material.  There  are  many  variables  that 
can  be  adjusted  when  operating  a  piston  corer,  such  as  free-fall  height,  core  liner  material,  weight,  piston  friction, 
rigging  and  winch  speed.  In  order  to  develop  a  capability  at  SACLATNCEN  to  collect  longer  cores  with  minimal 
disturbance,  the  aforementioned  variables  were  adjusted  in  a  systematic  manner  in  order  to  determine  their  relative 
effects.  During  the  Coring  Engineering  Trial  in  March  1999,  multiple  cores  were  collected  at  an  experimental  site  in  the 
Capraia  basin,  north  of  Elba  Island,  where  a  geoacoustic  model  has  been  developed  based  on  a  time  domain  inversion 
of  wide-angle  reflection  data  and  frequency  domain  modelling  of  bottom  loss  data.  The  disturbance  of  the  seabed 
material  during  the  coring  process  may  have  an  adverse  effect  on  its  physical  properties,  such  as  sound  speed, 
magnetic  susceptibility  and  bulk  density.  Accordingly  the  amount  of  compression  in  each  core  has  been  estimated  by 
an  analysis  of  magnetic  susceptibility  data,  correlating  and  aligning  nulls  in  the  response.  Laboratory  measurements 
of  sound  speed  and  bulk  density  on  three  of  the  cores  have  been  compared  with  each  other  and  with  the  geoacoustic 
model. 

From  a  seabed  dominated  by  silt  and  clay  material  with  some  thin  layers  containing  shells,  the  properly  configured 
piston  corer  was  able  to  recover  cores  of  5  to  6  min  length  with  a  compression  of  approximately  10%.  When  it  is  not 
properly  configured  or  does  not  function  properly,  the  compression  may  be  considerably  higher,  25  to  35%.  Critical 
factors  in  the  operation  of  the  piston  corer  included  the  piston  design,  the  piston  friction  against  the  liner,  the  strength 
of  the  shear  pins,  the  cable  lengths  and  the  winch  speed.  Factors  that  are  less  significant  include  the  liner  material  and 
free-fall  height.  The  laboratory  measurements  of  bulk  density  are  higher  that  those  determined  for  the  geoacoustic 
model,  but  may  be  explained,  in  part,  by  the  compression  of  the  material  during  coring. 
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Tesei,  A.,  Maguer,  A.,  Zerr,  B. 

.Multiple-aspect  acoustic  scattering  analysis  of  fluid-filled  cylindrical  targets  in  water. 

Research  on  mine  classification  has  been  focusing  for  some  years  at  SACLANTCEN  on  new  methods  based  on  the 
broadband  analysis  of  target  echo  at  low  to  intermediate  frequencies  in  time  and  frequency  domains. 

The  proposed  methodology  has  been  presented  in  previous  works  being  limited  to  multiple  aspect  target  scattering 
analysis  of  rigid  features  in  the  time  domain,  integrated  with  single-aspect  (broadside)  elastic  feature  analysis  in  the 
frequency  domain.  This  approach  allows  first  the  classification  between  man  made  and  natural  objects  based  on 
target  external  shape  and  rough  estimation  of  its  size  and  then,  if  the  target  is  classified  as  man-made,  its  eventual 
identification  by  estimating  internal  target  parameters,  such  as  filling,  shell  wall  material  and  thickness.  This  low- 
frequency  approach  allows  therefore  discrimination  between  objects  of  the  same  shape  and  size,  but  with  different 
internal  structures. 

The  present  work  proposes  an  extension  of  that  approach  to  the  multiple-aspect  analysis  of  both  rigid  and  elastic 
features  also  in  the  frequency  domain.  This  extension  is  addressed  to  make  the  basic  method  more  robust  and 
accurate.  Its  potential  is  demonstrated  here  on  simulated  and  real  data  for  targets  in  the  free  space  (i.e.,  suspended  in 
the  water  column).  The  low-frequency  approach  should  allow  the  method  applicability  also  to  partially  and  completely 
buried  targets,  as  will  be  demonstrated  in  future  work. 

Hence,  future  modelling  and  experiments  will  be  dedicated  to  the  extension  and  validation  of  the  methodology  for  the 
classification/identification  of  targets  either  lying  proud  on  the  sea  bottom  or  buried  in  the  sediments.  Further,  methods 
for  the  automatic  extraction  of  the  selected  scattering  features  will  be  investigated. 
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Zerr,  B., Tesei,  A.,  Maguer,  A.,  Houston,  B.,  Sletner,  P.A. 

Proud  target  classification  based  on  multiple  aspect  broadband  low  frequency  response. 

The  aspect  dependence  of  the  acoustic  signature  has  been  demonstrated  to  be  an  essential  indicator  to  the 
discrimination  between  man  made  and  natural  underwater  objects.  A  classification  method  has  been  define  using  the 
variation  with  incidence  angle  of  the  acoustic  waves  scattered  by  an  elastic  object.  As  the  experiment  conducted  in  a 
basin  on  free  field  cylinders  has  given  encouraging  results,  more  realistic  acoustic  measurements  have  been  conducted 
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on  natural  and  manufactured  objects  positioned  on  the  seabed.  The  external  shape ,  extracted  from  its  reflection  map 
reconstructed  by  tomography,  allows  selection  of  candidate  objects  for  detailed  analysis  of  their  scattering  properties. 
The  resonance  scattering  analysis,  limited  to  selected  aspects  in  its  original  version  (e.g.  broadside  fora  cylindrical 
shape)  has  been  extended  to  incorporate  aspect  varying  features.  The  variation  with  incidence  of  the  acoustic  wave 
diffracted  by  object  discontinuities  has  also  been  introduced.  This  paper  reviews  the  classification  technique  and 
describes  the  SACLANT CEN-NRL  TASCOE  (TArget  Scattering  in  Controlled  Environment)  experiment  conducted  in 
water  depth  of  15  mat  Marciana  Marina  (Elba,  Italy)  in  October  98.  The  results  obtained  from  the  data  analyzed  show 
that  the  aspect  dependence  of  the  acoustic  waves  scattered  by  elastic  objects  (ka  2-20)  allows  clear  discrimination 
between  manufactured  and  natural  objects. 
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Levesque,  I.,  Bondaryk,  J. 

Performance  issues  concerning  Doppler-only  localization. 

Target  localization  can  be  achieved  using  only  frequency  information  obtained  from  source-receiver  platforms.  The 
numerical  method  proposed  achieves  machine  precision  accuracy  in  computing  target  position  and  velocity  from 
multiple  simultaneous  Doppler-shift  measurements.  The  area  over  which  localization  can  be  performed  depends  very 
strongly  on  the  configuration  in  which  the  sensors  are  placed  and  the  drop-off  in  performance  is  step-like.  A  negligible 
decrease  in  performance  is  encountered  when  switching  from  “N-source,  N-sensors"  to  a  “N-sensor,  single-source " 
configuration.  The  number  of  sensors  to  be  used  is  an  important  criterion  and  it  is  found  that  increasing  this  beyond  4 
is  only  worthwhile  with  certain  layout  configurations.  Uncertainties  in  the  sensor  locations  decrease  the  area  over 
which  accuracy  is  acceptable;  a  10-meter  uncertainty  translates  into  a  5%  loss  in  coverage.  Noise  in  the  Doppler 
information  decreases  accuracy  gradually;  one  should  expect  an  error  3  times  greater  with  the  addition  of  10%  noise. 
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Holland,  C.W ,  McDonald,  B.E. 

Shallow  water  reverberation  from  a  time  reversed  mirror. 

The  ability  to  re-focus  energy  using  a  time  reversed  mirror  in  the  ocean  has  been  demonstrated  by  Kuperman  et  al. 
( 1 998),  Song  ()  and  others.  The  metric  in  these  experiments  has  been  the  focal  spot  itself.  Another  potentially  useful 
metric  of  the  time-reversed  mirror  concept  is  the  reverberation.  Exploring  the  reverberation  has  at  least  two  important 
aspects.  First,  the  reverberation  is  an  important  consideration  for  potential  tactical  applications  inasmuch  as  it  is  the 
target  illumination-to-reverberation  ratio  that  is  important  rather  than  just  the  characteristics  of  the  spot  size.  Second, 
the  time-reversed  mirror  has  the  potential  to  be  useful  as  a  scientific  tool  to  study  shallow  water  bottom  scattering. 
This  report  presents  results  from  shallow  water  reverberation  measurements  at  3500  Hz  designed  to  meet  these  two 
goals.  Reverberation  reduction  was  not  observed  when  the  focus  was  in  the  water  column  (far  from  the  boundaries). 
In  fact,  an  increase  in  the  reverberation  was  observed.  One  conclusion  from  these  results  is  that  in  order  to  exploit  the 
FAF concept  for  systems  concepts,  (perhaps  substantial)  azimuthal  resolution  is  required.  The  data  do  show  a  successful 
demonstration  of  a  focused  field  at  the  seafloor  boundary.  Thus  it  appears  that  FAF  can  be  exploited  as  a  scientific  tool 
for  studying  long  range  scattering  in  shallow  water. 
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Akal,  T.,  Di  lorio,  D.,  Guerrini,  P.,  Boni,  P.,  Ferla,  C.,  Hodgkiss,  W.S.,  Song,  H.C.,  Kuperman,  W.A., 
Jackson,  D.R. 

Range  and  time  dependence  of  focused  acoustic  field  and  iterative  focusing  at  450  Hz. 

A  phase  conjugate  “mirror"  time  reverses  the  incident  signal  precisely  returning  it  to  its  original  source  location.  This 
phenomenon  occurs  independent  of  the  complexity  of  the  medium.  The  time  reversal  process  can  be  accomplished 
by  the  implementation  of  a  retransmission  procedure.  A  signal  received  at  an  array  is  time  reversed  and  retransmitted. 
A  full  water  column  source  array  excited  by  the  phase  conjugated  (time-reversed)  signal  received  at  the  array  position 
will  focus  at  the  position  of  the  radiating  target.  The  medium  fluctuations  are  embedded  in  the  received  signal  so  that 
if  retransmission  can  occur  on  a  time  scale  less  than  the  dominant  fluctuations,  the  medium  variability  will  be  eliminated 
since  one  back  propagates  and  “undoes"  the  variability. 

SACLANTCEN  and  MPL  under  a  Joint  Research  Program  are  investigating  phase  conjugation  process  on  Focused 
Acoustic  Fields  (FAF)  since  1996.  Focusing  at  450  Hz  was  demonstrated  in  a  1996  experiment.  1996  experiment  was 
the  first  demonstration  of  focusing  with  a  phase  conjugate  array  in  an  ocean  environment.  The  following  year  a  second 
experiment  has  been  carried  out  where,  we  extended  the  range  of  focus  from  the  earlier  result  of  6  km  out  to  30  km, 
verified  anew  technique  to  refocus  at  ranges  other  than  that  of  the  probe  source  and  demonstrated  the  time  dependence 
of  the  FAF  successfully  with  extensive  environmental  data. 

In  this  report  we  present  the  results  of  the  final  low  frequency  (445  Hz)  Focused  Acoustic  Field  (FAF)  experiment  and 
describe  the  iterative  time-reversal  process  developed  with  low  frequency  FAF  data. 
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Sellschopp,  J. 

Rapid  environmental  assessment  for  naval  operations. 

Rapid  Environmental  Assessment  (REA)  is  a  methodology  that  is  being  implemented  in  order  to  close  knowledge 
gaps  and  to  provide  useful  environmental  information  in  a  tactically  relevant  time  frame.  REA  surveys  are  set  up  for 
operational  needs  rather  than  to  give  a  full  scientific  picture.  Emphasis  is  on  an  optimal  organizational  structure ,  fast 
data  processing  tools  and  modern  data  communication  channels.  Distributed  data  processing  and  product  generation 
is  complemented  by  a  particular  data  fusion  centre  in  an  Internet-like  network. 

Specific  instrumentation  for  REA  is  developed  under  the  premise  of  immediate  data  availability.  New  scientific  methods 
will  reduce  the  time  spent  for  data  collection.  Techniques,  concepts  and  procedures  must  be  tested,  trained  and 
developed  in  REA  exercises  that  may  or  may  not  be  connected  with  military  exercises.  The  Rapid  Response  surveys 
are  described  as  examples  for  extensive  REA  effort  prior  to  naval  exercises. 
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Hughes,  D. 

Aspects  of  card io id  processing. 

It  is  well  known  that  conventional  active  towed  arrays  suffer  from  ‘The  Ambiguity  Problem’ in  which  it  is  impossible  to 
distinguish  returns  from  port  or  starboard.  Although  operational  methods  for  overcoming  this  problem  exist  they  are 
time-consuming  and  cannot  be  used  for  a  single  ping.  One  proposed  method  of  coping  with  these  difficulties  is  to  use 
a  triplet  array  in  which  the  direction  of  arrival  of  a  signal  can  be  ascertained  on  a  single  ping. 

This  report  details  a  broadband  algorithm  which  has  been  developed  and  is  currently  being  used  at  the  SACLANT 
Undersea  Research  Centre.  The  algorithm  was  specifically  designed  to  be  used  on  a  delay-type  beamformerin  which 
shading  cannot  be  applied  at  the  beamlevei.  A  second  algorithm  employing  the  full-flexibility  of  a  3d  beamformer  is 
also  described  and  the  two  algorithms  are  compared  for  their  performance  in  terms  of  left-right  suppression  and  SNR 
for  which  analytical  expressions  are  computed.  It  is  shown  that  both  of  the  algorithms  have  desirable,  but  different 
features.  The  algorithms  are  applied  to  data  taken  from  the  trial  BACCHUS’98  and  the  results  obtained  are  used  to 
help  validate  the  theoretical  results  and  to  give  an  indication  of  the  performance  which  is  achievable  for  the  triplet 
system.  We  show  that  for  the  Centre  algorithm  a  typical  achievable  left-right  suppression  is  greater  than  15dB. 
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Laterveer,  R. 

MRF  segmentation  for  low  frequency  active  sonar:  further  results. 

The  use  of  low  frequency  active  sonar  in  shallow  water  leads  to  large  numbers  of  clutter  detections. 

This  high  false  alarm  rate  can  overload  automatic  tracking  and  classification  algorithms.  Traditional  detection  algorithms 
operate  on  each  beam  output  individually  searching  for  targets  at  all  ranges.  However,  the  target  echo  and  bottom 
features  may  extend  over  several  beams,  either  because  a  reflector  is  extended  over  space  or  because  of  the  sidefobe 
structure  of  the  beamformer.  This  suggests  the  association  of  detections  over  bearing,  e.g.  apply  image  processing  to 
the  range-bearing  sonar  data. 

A  previous  study  described  an  automatic  method  of  image  segmentation  based  on  a  Markov  random  field  (MRF) 
model  to  reduce  clutter.  The  segmentation  is  treated  as  a  labelling  problem,  assigning  to  each  range-bearing  cell 
either  a  target  or  background  label,  removing  small  objects  which  do  not  exhibit  the  correct  signature  over  beams . 
Separate  detections  corresponding  to  one  large  reflector  are  combined  and  removed  if  they  are  too  large  to  be  a 
submarine . 

In  this  report,  the  algorithm  is  evaluated  on  additional  data.  On  a  single  ping  basis  the  MRF  segmentation  shows 
slightly  lower  detection  performance  at  the  same  number  of  false  alarm  objects  as  the  Page  test  detector.  MRF 
segmentation  still  shows  potential  for  fixed  feature  removal  over  multiple  pings  and  will  be  studied  in  the  feature. 
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van  Velzen,  M.R.,  Laterveer,  R. 

Performance  measurement  in  active  sonar  using  object  counting. 

At  present  two  measures  of  performance  (MOP)  are  frequently  being  used  at  the  Centre  for  the  evaluation  of  detection 
algorithms  in  Low  Frequency  Active  Sonar  (LFAS)  data  processing.  The  performance  of  the  basic  processing  elements 
of  every  active  sonar,  consisting  of  beamforming,  matched  filtering  and  normalization,  is  expressed  by  the  standard 
Receiver  Operating  Characteristic  (ROC)  curve,  which  relates  the  probability  of  detection  to  the  probability  of  false 
alarm  per  range-bearing  cell.  However,  for  highly  time  varying  data  or  after  some  clustering  has  been  applied  to  the 
data,  performance  is  usually  expressed  as  the  probability  of  detection  versus  the  number  of  false  target  like  objects 
per  ping.  The  existents  of  two  MOPs  depending  on  the  type  of  data  or  the  level  of  processing  is  unsatisfactory  as  it 
prevents  quantifying  the  improvement  achieved  at  every  successive  processing  stage . 
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To  obtain  the  ROC  curve  for  experimental  data ,  target  echo  and  background  are  often  separated,  using  a  priory 
knowledge  of  the  target  position,  to  measure  the  probability  of  detection  and  probability  of  false  alarm  separately  For 
an  objective  measure  of  detection  performance  target  and  background  should  not  be  separated  as  this  may  hide 
some  cross  effects  like  the  increased  false  alarm  probability  due  to  side  lobes  in  the  beam  direction  of  the  target  echo. 

A  final  measure  of  detection  performance  should  also  be  a  function  of  the  accuracy  which  the  actual  position  of  the 
target  can  be  found.  This  could  show,  for  example ,  the  improved  detection  performance  due  to  fine  bearing  estimation. 
Positioning  errors  also  plays  an  important  role  in  multi-statics  where  detections  from  different  platforms  will  have  to  be 
compared. 

In  this  report  we  propose  a  method  that  counts  all  connected  objects  in  the  data  and  defines  a  unique  point  in  each 
object  that  defines  the  exact  location  of  the  detection.  The  method  can  be  applied  at  every  stage  of  processing,  is  a 
function  of  the  accuracy  with  which  the  target  location  can  be  found  and  no  target  location  information  is  being  used 
in  counting  detections  and  false  alarms.  The  resulting  MOP  gives  the  probability  of  detecting  the  target  within  a  box 
around  the  actual  target  position  versus  the  number  of  false  alarms  per  ping. 
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Lyons,  A.P.,  Fox,  W.L.J.,  Hasiotis,T.,  Pouliquen,  E. 

Characterization  of  the  two-dimensional  roughness  of  shallow-water  sandy  seafloors. 

Surface  roughness  is  a  fundamental  seafloor  property  affecting  a  variety  of  physical  phenomena  including  sediment 
transport  and  the  interaction  of  acoustic  energy  with  the  seafloor.  Characterization  of  seafloor  surface  roughness  and 
its  dynamics  is  therefore  essential  for  understanding  and  quantifying  the  influence  of  sediment  microtopography.  Field 
measurements  have  been  taken  recently  with  an  end-to-end  digital photogrammetry  system  providing  quantitative, 
two-dimensional  seafloor  surface  roughness  measurements  on  spatial  scales  of  approximately  a  millimeter  to  a 
meter.  Results  of  these  measurements  have  shown  that  sediment  surfaces  in  shallow  water  are  often  anisotropic  and/ 
or  exhibit  non-Gaussian  height  distributions,  both  of  which  have  the  potential  to  strongly  affect  high-frequency  seafloor 
acoustic  scatter.  For  these  kinds  of  surfaces,  simple  roughness  parameters  such  as  rms  height  or  the  slope  and  offset 
of  a  power  law  representation  of  the  power  spectra  will  not  give  a  sufficiently  complete  description.  Two-dimensional 
statistical  models  are  needed  to  capture  the  anisotropic  nature  of  sediments  with  oriented  features,  while  for  seafloors 
with  peaked  forms,  it  is  the  phase  information  in  the  frequency  domain  that  is  required,  as  this  controls  the  shape 
characteristics  of  a  surface.  Characterization  of  seafloor  roughness  based  on  these  ideas  will  be  presented  using 
results  from  the  digital  photogrammetry  system. 
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Cdr  B.  Williams,  USN  (from  July) 

Scientific  Division  Coordinator 

LCdr  J.R.  Staveley  (UKN)  from  July) 

Security  Officer 

Cdr  C.  Gabrini,  ITN  (until  August) 

CdrG.  Piaggio,  ITN  (from  September) 

Budget  Officer 

LCdr  R.  Reilly  (USN) 

ASW  Programme  Officer 

Cdr  B.  Egenberg  (NON) 

Mine  Warfare  Officer 

LCdr  J.R.  Staveley  (UKN)  (until  June) 
LCdrO.  Molino,  ITN  (from  July) 

ADP  Security  Officer 

LtJ.G.  M. Shannon  (USN) 

MILOC  Survey  Officer 

D.  Schaa  (Civilian)  (GE) 

FINANCE  OFFICE 


Controller: 

M.  Baldassari  (IT)  (until  January) 

L.  Parise  (IT)  (from  February) 

Fiscal  Officer: 

J,  Jenkins  (US)  (until  August) 

M.  Newman  (US)  (from  August) 

Internal  Auditor: 

G.  Lambie  (UK) 

PERSONNEL  &  ADMINISTRATION  OFFICE 


Head  of  PAD 

R.  Albini  (IT) 

Logistics  &  Supply 

U.  Varlese  (IT) 

Information  Services 

A.W  Green  (UK) 

INFORMATION  TECHNOLOGY  OFFICE 


Head  of  ITO 

R.  Seynaeve  (BE) 

Senior  ADP  Systems  Spec 

P.  Colucci  (IT) 

SHIP  MANAGEMENT  OFFICE 


Head  of  SMO 

C.  Gobey  (UK) 

Senior  Engineer 

G.  Soncini  (IT)  (until  Feb) 

A,  Spairani  (IT)  (from  March) 

T-Boat  Captain 

A.  Spairani  (IT)  (until  August) 

A.  lacono  (IT)  (from  September) 
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CHIEF  SCIENTIFIC  DIVISION 


Deputy  Director 


ACOUSTICS  DEPARTMENT 


Head  of  ACD 

W.  Roderick  (US) 

Senior  Principal  Scientist 

F.  Jensen  (DA) 

Senior  Principal  Scientist 

N.  Pace  (UK) 

Principal  Scientist 

T.  Akal  (TU) 

Principal  Scientist 

D.  Burnett  (US) 

Principal  Scientist 

J.P.  Hermand  (BE) 

Principal  Scientist 

J.  Watermann  (GE)  (until  April) 

Principal  Scientist 

C.  Harrison  (UK)  (from  March) 

Senior  Scientist 

O.  Bergem  (NO)  (until  February) 

Senior  Scientist 

M.  Ferla  (IT) 

Senior  Scientist 

C.  Holland  (US) 

Senior  Scientist 

K.  LePage  (US) 

Senior  Scientist 

R.  Tyce  (US)  (from  August) 

Scientist 

A.  Lyons  (US) 

Scientist 

P.  Nielsen  (DA) 

Scientist 

E.  Pouliquen  (FR) 

Scientist 

M.  Siderius  (US) 

Head  of  ETD 

A.  Barbagelata  (IT)  (until  February) 

0.  Bergem  (NO)  (from  March) 

Principal  Engineer  (Head,STO) 

F.  de  Strobel  (IT) 

Principal  Engineer  (Head, ENG) 

P.  Guerrini  (IT) 

Senior  Engineer 

L.  Gualdesi  (IT) 

Senior  Engineer 

E.  Michelozzi  (IT) 

Engineer 

V.  Grandi  (IT) 

Engineer 

B.  Miaschi  (IT) 

Engineer 

P.  Sletner  (NO) 

Engineer 

R.  Stoner  (UK) 

Engineer 

L.  Troiano  (UK) 

Head  of  OCD 

J,  Sellschopp  (GE) 

Principal  Scientist 

F.  Askari  (US) 

Principal  Scientist 

B.  McDonald  (US) 

Senior  Scientist 

R.  Onken  (GE) 

Scientist 

A,  Alvarez  (SP)  (from  June) 

Scientist 

D.  Di  lorio  (CA)  (until  April) 

Scientist 

J.  Osier  (CA)  (until  September) 

Sr  Admin  (Scientific  Specialist) 

F.  Spina  (IT) 

- ■ 

OPERATIONAL  RESEARCH 

DEPARTMENT 

Head  of  ORD 

K.  Pye  (UK) 

Principal  Scientist 
Principal  Scientist 
Principal  Scientist 

J.  Redmayne  (UK) 

D.  Ruskin  (US)  (until  August) 

P.  Simcock  (UK)  (from  September) 

Senior  Scientist 
Senior  Scientist 

P.  Simcock  (UK)  (until  August) 

E.  Verhoeff  (NL) 

Scientist 

G.  Arcieri  (IT) 

Naval  Scientist 

P.  Golmayo  (SP) 

3 

Head  of  SSD 

J.  Ziegenbein  (GE) 

Senior  Principal  Scientist 

A.  D'Amico  (US) 

Senior  Principal  Scientist 

M.  Pinto  (FR) 

Principal  Scientist 

G.  Field  (UK) 

Senior  Scientist 

G.  Davies  (UK)  (from  February) 

Senior  Scientist 

D,  Grimmett  (US) 

Senior  Scientist 

D.  Hughes  (UK) 

Senior  Scientist 

A.  Maguer  (FR)  (until  October) 

Senior  Scientist 

B.  Zerr  (FR) 

Senior  Scientist 

M.  Van  Velzen  (NL) 

Scientist 

A.  Bellettini  (IT)  (from  December) 

Scientist 

J.  Bondaryk  (US) 

Scientist 

S.  Fioravanti  (IT) 

Scientist 

G.  Haralabus  (GR) 

Scientist 

R.  Laterveer  (NL) 

Scientist 

L.  Mozzone  (IT) 

Scientist 

A.  Tesei  (IT) 

Sr  Admin  (Scientific  Specialist) 

E.  Bovio  (IT) 

Sr  Admin  (Scientific  Specialist) 

W.  Zimmer  (GE) 
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Visitors 


January 

Professor  Roberto  Schmid 

Professor  Mario  Pavan 

Dr  Gianni  Pavan 

University  of  Pavia 

March 

Group  of  Italian  Navy  officers  (ASW 
and  MW  specialists) 

ITN  Academy,  Livorno 

Officers  from  Italian  Navy 

Hydrographic  Institute  of  Genova 

Hydrographic  Institute  of  Genova 

Jack  F.J.  Venneulen 

TNO-FEL  Netherlands 

DrTsihYang 

NRL  Washington,  D.C. 

Dirk  G.  Simons  and  Mirjam  Snellen 

TNO-FEL  Netherlands 

Richard  J.  Brind  and  Alex  S.  Penfold 

DERA,  U.K. 

April 

RADM  L.C.  Baucom,  USN 

SACLANT  ACOS  Plans  and  Policy 

VADM  Manlio  Galliccia,  ITNA 

Commander  in  Chief,  Upper  Tyrrhenian 
District,  La  Spezia 

May 

Group  of  Norwegian  Navy  officers 

Norwegian  Tactical  Training 

(ASW  and  MW  specialists) 

Establishment 

Dr  M.l.  Yarymovych 

Chairman,  NATO  Research  and 
Technology  Board,  Paris 

Dr  E.A.  van  Hoek, 

Director,  NATO  Research  and 

Technology  Board,  Paris 

June 

Mr  J.  Johnsen 

SCNR  Norway 

Mr  Thorsen 

NDRE  Oslo 

Dr  H.  Schneider  and  Dr  Anton  Homm 

FWG  Kiel 

July 

Dr  J.  Keil 

Executive  Director,  NUWC  Newport 

Dr  B.  Kittredge 

Deputy  Director,  NUWC  Newport 

August 

Dr.  D.  Wyn  Williams 

Operations  Director,  Sea  Systems,  DERA 

September 

Dr  Richard  W.  Spinrad 

Technical  Director,  Office  of 
Oceanographer  of  Navy,  U.S. 

Thomson  Marconi  Representatives 

France 

Prof.  F.  Abail  and  CDR  P.K.  Veldhuis 

Netherlands  Naval  College 

Mr.  F.  Asbeck 

Western  European  Satellite  Centre 

Dr  G.  Notarbartolo  di  Sciara 

President  ICRAM 

October 

RADM  Dino  Franco  Vene,  ITNA 

President  of  MARIPERMAN,  La  Spezia 

RADM  Paul  Gaffney,  II,  USNA 

French  Officers  from  FS  DE  GRASSE 

Chief  of  Naval  Research 

November 

RADM  Gino  Bizzari,  ITNA 

COMFORDRAG,  La  Spezia 

ADM  E.  Bengtsen 

GDRE  J.  Parker 

CDR  G.  Stamp 

SACLANT  HQ 

VADM  Jose  Romero  Iglesias,  SPNA 

Visit  of  STANAVFORMED  ASW 

Officers 

Deputy  COMNAVSOUTH 

December 

Dr  Wyn  Williams 

UK  member  SCNR 

Mr  Michael  White 

DERA 

Mr  N.  Fisher  and  Mr.  J.  Downing 

MOD  UK 
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Meetings 


Scientific  Committee  of  National  Representatives 


and  National  Liaison  Officers 


BELGIUM 

National  Representative 

LCDR  Claude  L.  Renard,  BENA 

Section  Systemes  d’Armes 

Lutte  anti-sous-marine 
de  PEtat-Major  de  la  Marine,  Brussels 

CANADA 

Acting  National  Representative 

Mr  Warren  C.E.  Nethercote 

Deputy  Director-General,  Defence  Research 
Establishment  Atlantic,  Dartmouth,  Nova  Scotia 

National  Liaison  Officer 

Dr  Dr.  Dale  Ellis 

Defence  Research  Establishment  Atlantic, 
Dartmouth,  Nova  Scotia 

DENMARK 

National  Representative 

Rear  Admiral  Niels  Mejdal,  DENA 

Defence  Command  Denmark,  Vedbaek 

National  Liaison  Officer 

Mr  Bjarne  Damsgaard 

Danish  Defence  Research  Establishment, 
Copenhagen 

FRANCE 

National  Delegate 

ICA  Daniel  Long 

Service  des  Programmes  Navals, 
chef  du  Departement  Lutte  Sous  la  Mer 

Paris 

GERMANY 

National  Representative 

Mr  Siegfried  Tympe 

Bundesministerium  der  Verteidigung,  Bonn 

National  Liaison  Officer 

Dr  Dirk  Tielbuerger 

Forschungsanstalt  der  Bundeswehr  fur 
Wasserschall-und  Geophysik  (FWG),  Kiel 

GREECE 

National  Representative 

Captain  Anastasios  Sklavidis,  HENA 

Hellenic  Navy  Hydrographic  Office,  Holargos, 
Athens 

ITALY 

National  Representative 

Rear  Admiral  Dino  Nascetti,  ITNA 

Direttore,  Arsenale  Militare  Marittimo,  La  Spezia 

National  Liaison  Officer 

Commander  Carmelo  Di  Natale,  ITNA 

MARISTAT,  Rome 

NETHERLANDS 

National  Representative 

Mr  Coenraad  M.  Ort 

Physics  and  Electronics  Laboratory, 

FEL-TNO,  The  Hague 

National  Liaison  Officer 

Captain  Gijs  J.  Van  Beeck  Calkoen,  NENA 

Director  of  Naval  Research  and  Development, 
Royal  Netherlands  Navy,  The  Hague 

NORWAY 

National  Representative 

Mr  Jarl  Johnsen 

Norwegian  Defence  Research,  Establishment, 
Horten 

PORTUGAL 

National  Representative 

LCDR  E.  Ferreira  Coelho,  PONA  Instituto 
Hidrografico,  Lisbon 
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SPAIN 

National  Representative 

Captain  Jorge  Juan  Moreno  Sanmartin, 

SPNA 

JEFE  Spanish  MCM  Flotilla,  Cartagena 

National  Liaison  Officer 

LCDR  Juan  A.  Rico  Palma,  SPNA 

Instituto  Hidrografico  de  la  Marina,  Cadiz 

TURKEY 

National  Representative 

Eng.  Captain  Nazim  Qubukgu,  TUNA 

Head  of  the  Department  of  Navigation, 
Hydrography  and  Oceanography,  Istanbul 

National  Liaison  Officer 

Eng.  LCDR  H.  Basaran,  TUNA 

AKAG  Bsk,  Gdlcuck/Kocaeli 

UNITED  KINGDOM 

National  Representative 

Dr.  D.  Wyn  Williams 

Operations  Director,  Sea  Systems  Sector, 
Defence  Evaluation  and  Research  Agency, 
Portsdown  West 

National  Liaison  Officer 

Mr  David  Lewis 

DSc(SEA)3.3.,  Ministry  of  Defence,  London 

UNITED  STATES 

National  Representative 

Dr  Eric  0.  Hartwig  (NRL  7000) 

Associate  Director  of  Research  for  Ocean  and 
Atmospheric  Science  and  Technology,  Naval 
Research  Laboratory,  Washington  D.C. 

Alterna  te  Na  tiona  1  Represen  ta  ti ve 

Dr  Steven  E.  Ramberg  (ONR-32) 

Head,  Ocean,  Atmosphere  and  Space  Science 
and  Technology  Department,  Office  of  Naval 
Research,  Arlington,  VA 

Na  tiona  i  L  ia  ison  Officer 

Commander  Scott  M.  Tilden,  (ONR  321 B) 
Littoral  Warfare  Advanced  Development, 

Project  Manager  and  Military  Deputy,  Sensing 
and  Systems  Division,  Office  of  Naval 

Research,  Arlington,  VA 

SECGEN  NATO 

Representative 

RADM  Guillermo  Leira,  SPNA 

Deputy  Assistant  Secretary  General,  Defence 
Support,  NATO  Headquarters,  Brussels 

NAMILCOM 

Representative 

CoL  Pat  Nutz,  USAF 

Logistics,  Armaments  &  Standardization 

Division,  NATO  Headquarters,  Brussels 

SACLANT 

Representative 

Rear  Admiral  Dieter  George  Leder,  GENA 

Assistant  Chief  of  Staff  (Operations), 

Representative 

Rear  Admiral  Eivind  Bengtsen,  DANA 

Assistant  Chief  of  Staff  (Resources), 

Acting  Representative 

Commodore  Jeremy  Parker,  UKNA 

Deputy  Assistant  Chief  of  Staff  (Operations) 

Acting  Representative 

Captain  Eugene  Alleman,  BENA 

Deputy  Assistant  Chief  of  Staff  (Resources) 

Liaison  Officer 

Commander  Gordon  Stamp,  UKNA 

Personnel  by  category  and  nationality 
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Personnel  Strength  by  category  and  nationality  at  31  December  1999  (With  aggregate  Scientist  Numbers  1959  -  1 
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Including  Director  and  Deputy  Director 
Since  1  May  1959 


5»l  o\ 


Underwater  Warfare  Workshop,  1  -  3  September,  1999 


l 


[■  SG-21  and  NG-2  meetings,  20-21  and  22  -  24  September,  1999 


NATO  SACLANT  Undersea  Research  Centre 
Viale  San  Bartolomeo,  400 
19138  La  Spezia,  Italy 
Address  for  mail  from  USA  and  Canada: 

SACLANTCEN 
APO  AE  09613-5000 

telephone  +39  0187  5271  or  527  +  known  extension 
facsimile  +39  0187  527  700 
http://www.saclantc.nato.int 
email:  library® saclantc.nato.int 

Applying  science  to  NATO  maritime  operational  requirements  since  1959 


Monthly  infrared  images  of  sea  surface  temperature,  during  1999,  from  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  satellite  sensor  depict  an  anticyclonic 
gyre  in  the  eastern  Mediterranean  Sea,  southeast  of  Crete.  The  images  commence  in  the  upper  left  with  January  and 
end  in  the  lower  right  with  December.  Dark  blue  designates  cooler  surface  temperatures,  black  represents  cloud 
cover.  Shown  is  the  lerapetra  Gyre  with  clear  seasonal  variations.  Maximum  temperature  differences  occur  in  the 
summer  when  cooler  water  flowing  southeast  from  the  Aegean  Sea  mixes  with  the  warmer  Eastern  Mediterranean 
water.  There  are  gyres  throughout  the  western  and  eastern  Mediterranean,  which  are  a  result  of  atmospheric  forcing 
and  complex  dynamics  of  water  masses  entering  and  leaving  the  Mediterranean  Sea.  Satellite  remote  sensing  is  one 
component  of  Rapid  Environmental  Assessment  (REA),  which  provides  valuable  information  on  sea  surface  parameters 
in  littoral  and  deeper  water.  Progress  is  being  made  in  applying  these  data  to  the  prediction  of  oceanographic 
variables  which  affect  undersea  warfare  operational  capabilities,  (see  pages  7  and  16  to  21). 


